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Abstract 

In this paper, we introduce a class of spacetimes {Ai,g) which sat- 
isfy the vacuum Einstein equations and dynamically approach a Schwarz- 
schild solution of mass M, a class we shall call ultimately Schwarzschildean 
spacetimes. The approach is captured in terms of boundedness and decay 
assumptions on appropriate spacetime-norms of the Ricci-coefhcients and 
spacetime curvature. Given such assumptions at the level of k deriva- 
tives of the Ricci-coefBcients (and hence fc — 1 derivatives of curvature), 
we prove boundedness and decay estimates for k derivatives of curvature. 
The proof employs the framework of vectorfield multipliers and commuta- 
tors for the Bel-Robinson tensor, pioneered by Christodoulou-Klainerman 
in the context of the stability of the Minkowski space. We provide mul- 
tiplier analogues capturing the essential decay mechanisms (which have 
been identified previously for the scalar wave equation on black hole back- 
grounds) for the Bianchi equations. In particular, a formulation of the 
redshift-effect near the horizon is obtained. Morever, we identify a cer- 
tain hierarchy in the Bianchi equations, which leads to the control of 
strongly r-weighted spacetime curvature-norms near infinity. This allows 
to avoid the use the classical conformal Morawetz multiplier K, therby 
generalizing recent work of Dafermos and Rodnianski in the context of 
the wave equation. Finally, the proof requires a detailed understanding of 
the structure of the error-terms in the interior. This is particularly intri- 
cate in view of both the phenomenon of trapped orbits and the fact that, 
unlike in the stability of Minkowski space, not all curvature components 
decay to zero. 
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1 Introduction 

A major open problem in general relativity is to establish the non-linear stability 
of the Kerr family of solutions. Considerable mathematical progress towards 
this goal has been achieved in recent years, mainly by studying the linear wave 
equation, 

□ gV. = 0, (1) 

on fixed Schwarzschild [1, 2, 3] or Kerr-black hole backgrounds [4, 5, 6]. 

In particular, a precise and robust understanding of the role of the horizon 
for the decay mechanism of linear waves is now available within the framework 
of vectorficld multipliers and commutators, which is the common ground on 
which almost all of the currently available results stand. In the same context, 
there is also a good understanding of the obstruction to decay associated with 
the phenomenon of trapped orbits and the associated loss of derivatives in the 
estimates, at least for Schwarzschild- and Kerr black holes. We refer the reader 
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to [7] for a detailed discussion of these characteristic properties of black hole 
backgrounds. 

While the study of (1) provides important insights regarding the role of the 
geometry for the decay problem, its connection with the black hole stability 
problem is a-priori rather remote: (1) is certainly not the linearization of the 
Einstein equations with respect to a fixed black hole background g, but merely a 
"poor man's" linearization, which forgets entirely about the tensorial character 
of the original equations. 

Encouraged by the rapid progress regarding (1), we initiate in this paper 
a new approach to address the linear stability problem for black holes. This 
novel point of view will turn out to be much more intimately connected to the 
non-linear stability problem, as it is based entirely on a study of the Einstein 
and the Bianchi equations governing the metric evolution. 

The setting we suggest is the following. Consider a non-stationary black hole 
spacetime {A4,g), which satisfies the vacuum Einstein equations 



and in addition settles down to a stationary solution for late times. In the 
familiar geometrical framework, the latter assumption of approach will mani- 
fest itself in appropriate decay assumptions on deformation tensors and Ricci 
rotation coefficients, once an appropriate coordinate system in which decay is 
measured, has been fixed. 

Clearly, in view of (2) being satisfied, the Weyl-curvature tensor W of such 
a spacetime is equal to the Riemann-tensor and obeys the Bianchi equations 



Exploiting the rich structure of (3) one can then try to prove estimates for the 
curvature tensor from appropriate assumptions on the rates of approach for the 
metric. 

Now by the very definition of curvature, decay assumptions on the future 
asymptotic behavior of the metric and its derivatives will imply certain decay 
of the curvature tensor.^ Hence in order for the estimates arising from (3) to 
be useful, one should be able to establish more decay on the curvature than 
what immediately follows from the assumptions on the metric coefficients and 
deformation tensors. More precisely, in order to be valuable in a potential 
non-linear application, the decay proven for W should be sufficiently strong 
to stand a chance to eventually improve, via the null-structure equations, the 
assumptions initially made on the spacetime metric. This type of argument 
forms an integral part in the context of the bootstrap setting in which non- 
linear stability results are typically proven. 

In the present paper, we address the simplest example, for which the ideas 
outlined above can be successfully carried out: A spacetime approaching a 
Schwarzschild solution for late times. While this example is non-generic in 
some sense, a large class of such spacetimes is expected to exist. Moreover, as 
we will see, the problem introduced here already exhibits many of the difficulties 
that the analogous problem for an ultimately Kerr spacetime would do. 

^When we talk about decay of the curvature tensor here and in the future, we always 
understand this as decay to the curvature tensor of the spacetime we are approaching. 



R^.. (.9) = 0, 



(2) 



D'^Wai^jS (5) = . 



(3) 
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Roughly speaking, the statement we are going to prove in this context is 
the following. For sufficiently large k, houndedness of the -energy and the 
-spacetime-norm of k-derivatives of the Ricci rotation coefficients (and ap- 
propriate decay for lower order norms) implies houndedness and a degenerate 
version of integrated decay of the L^-norm of k derivatives of the curvature 
tensor. While in this paper we do not concern ourselves with the problem of 
improving the bounds on the Ricci coefficients from the curvature bounds estab- 
lished (such estimates should follow closely the ones in [8]), we will in addition 
prove decay estimates for lower order curvature norms in terms of norms on 
the Ricci-coefhcients. These estimates will reveal consistency of the approach 
in the sense that one does not have to assume stronger decay-rates on the Ricci- 
cocfficients than the ones one expects to eventually derive from the curvature 
bounds obtained via multiplier estimates from the Bianchi equations.^ 

Before we turn to a more precise statement of the theorem and the main 
difficulties associated with its proof, let us briefly outline the basic technique of 
obtaining estimates for W from (3). 

There exist natural energy currents for (3) or, more generally, (3) with an 
inhomogcneity Sp-ys on the right hand side. These currents arise from the so- 
called Bel-Robinson tensor 

Q m^p,8 = w«P7-w// + -^w^p^.^w// . (4) 

The tensor Q is symmetric and traceless and satisfies the divergence identity 
{Qc.p^sXPyz') = Af [W] + K^^^ [W] , (5) 

for any spacetime vectorfields X,y,Z. The term Ki^^^ [W] vanishes in the 
case of the homogeneous Bianchi equations. 

Moreover, Q {X, y, Z,n) > holds for any future directed causal vectorfields 
X, y, Z, n. If they are timelike, then Q {X, y,Z,n) in fact controls the sum of 
squares of all components of W. Note also that both the right hand side and 
the left hand side of (5)depend only on W and not derivatives thcrof. 

The basic strategy to derive estimates for (3) involves choosing appropriate 
vectorfields and integrating (5) over certain spacetime regions. For instance, if 
the vectorfields X,y, Z are Kilhng and J = (the homogeneous Bianchi equa- 
tions), the right hand side of (5) vanishes, and we in fact obtain a conservation 
law relating currents on future spacelike slices to those in the past. 

For our spacetime however, there are no exact Killing vectorfields but only 
what we are going to call "ultimately Killing fields". These are vectorfields. 



^As discussed below, the houndedness result for k derivatives of curvature would be much 
easier to prove if one were to assume strong decay in t for k derivatives of the Ricci-coefHcients, 
say -p^-decay of the L^-norm. However, this would be inconsistent in the above sense: 
Because of the trapping, multiplier estimates from the Bianchi equations can only produce 
■i-decay for k — 1 derivatives of curvature, which fails to match the decay-rate assumed on 
fc-derivativcs of the Ricci-coefBcients. Cf. the discussion in section 1.1.3. 
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whose deformation tensor does not vanish but decays to zero in time. The 
non-zero deformation tensor requires a careful analysis of the non-vanishing 
spacetime terms on the right hand side of (5).'^ 

We also note, schematically at least, the formula for commuting the Bianchi 
equations with vectorfields: 

D"(CxW) = ^■^K ■ DW + D^^\ ■ W (6) 

Here Cx & modified Lie-derivative which gives CxW the algebraic proper- 
ties of a Weyl-tensor (cf. section 3.2). Clearly, deriving energies for CtW will 
necessitate an analysis of the non- vanishing term Ki^'^'^ \W]- 

We are now in a position to explain at a heuristic level, some of the difficul- 
ties and challenges one encounters, when one applies the vectorfield techniques 
outlined in the previous paragraph to the problem we wish to study: 

1. Not all components of W decay for the spacetime under consideration 
(otherwise we would approach the Minkowski space). This causes difficul- 
ties both at the level of energies associated to (3) (as they do not decay), 
and at the level of the error-terms arising from commuting (3) with the 
approximate symmetries of the spacetime. As is already apparent from 
formula (6), in sharp contrast to the stability of Minkowski space, not all 
error-terms will exhibit a quadratic structure in which both components 
decay. 

2. What is the analogue of the rcdshift effect [1] for (3)? 

3. How is the obstruction for decay associated with trapped orbits captured 
for (3)? 

4. How docs one obtain decay (in the interior, at null infinity, ...) for Wl 

5. What are the appropriate (minimal?) assumptions on the Ricci rotation 
coefficients? 

As perhaps anticipated by the reader, the issues listed above are coupled to 
one another. 

1.1 Summary of the argument 

We proceed with a summary of the estimates carried out in the paper explaining 
how the above issues are addressed. In view of the frequently extensive and 
tedious formulae, we will regularly refer to the bulk of the paper. We hope that 
by this means, the introduction can also serve as a guide to the paper. 

1.1.1 Ultimately Schwarzschildean spacetimes 

Fix a regular coordinate system {t*,r,6,(j>) on the black hole exterior of a 
Schwarzschild metric gM with mass M. Consider perturbations of qm, such 
that the resulting metric g on the black hole exterior is what wc call ultimately 

^We remark at this point that it would not make sense to study (3) on a fixed Schwarz- 
schild background: Here the dynamics is trivial due to an algebraic constraint (known as the 
Buchdahl constraint in the literature [9]) on the Riemann- tensor. 
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Schwarzschildean to order k + 1 (cf. Definition 2.1). At the core of this definition 
are assumptions at the level of k derivatives of the Ricci-coefficients (i.e. k + 1- 
derivatives of the metric): In particular, both an energy on spacelike slices with 
null-ears (Er, cf. the figure in section 2.2) and a spacetime energy for the re- 
gion enclosed by two such slices are defined for the Ricci-rotation coefficients *H: 
These energies measure the approach of the U\ to their Schwarzschild values in 
the coordinate system {t*,r,9,4'). Roughly speaking, we assume boundedness 
of both of these energies for k derivatives, decay for fc — 1 derivatives, 
for fc — 2 derivatives and t"^ for fc — 3 or less derivatives (cf. (42) and (44)). 
The reason for this hierarchy is related to what we expect to be able to show 
for the energies involving curvature and will become more clear later. 

In addition, we assume the existence of a vectorfield T, satisfying g (T, T) = 
on the horizon and approaching the timelike (null on the horizon) Killing field 
of Schwarzschild in the sense that the deformation tensor associated with T 
decays. We will refer to the vectorfield T as being an ultimately Killing field. 
See Definition 2.3. 

1.1.2 Null-decomposition of the curvature tensor 

We will also pick a null-frame ei, 62, 63, 64 for g (5(63,64) — —2)), which is 
assumed to be appropriately close to a previously fixed Schwarzschild frame. We 
null decompose the curvature tensor with respect to that frame. In the standard 
notation of the subject, we obtain the components a, /?, p, cr, /3, a (cf. section 3.3). 
The equations (3) can then be written as equations for the null-components 
along null-directions (cf. (85)-(94)). We note that in Schwarzschild all these 
components except p = — are zero. 

1.1.3 The T-energy 

The first estimate that one may wish to obtain from (3) is, of course, the ana- 
logue of the energy estimate associated with the ultimately Killing field T. As 
p itself does not decay, and also in view of the fact that higher derivative es- 
timates need to be derived, we are going to commute the equation sufficiently 
many times with the vectorfield T. However, as is apparent from formula (6), 
irrespective of how many commutations are performed, there will always be a 
highest order error-term of the form p- C^~^DTr, which unlike most other terms, 
does not decay quadratically.'* As a consequence, in the energy estimate arising 
from (5) we will have to control an error-term of the form 

f p- (^^-^Dtt) -f^M/. (7) 

spacetime 

In view of the trapping, we will not be able to control the spacetime-integral 
of the highest order derivative term globally, which means that we will have 
to put it "L^ in space". It follows that to establish boundedness of the error- 
term at the highest level of derivatives (and hence boundedness of the fc*'' order 
Weyl energy), we would need very strong decay (at least -p-) of the L^-cnergy 
of fc derivatives of the deformation tensor. However, such a decay has only a 
chance of being obtained from the curvature components via elliptic estimates, 

^ There are also lower order terms which do not decay quadratically, in view of the fact 
that some derivatives of p do not decay. We will focus on the highest order term here. 
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if we can show the same decay for k — 1-derivatives of curvature. Unfortunately, 
with current vectorfield techniques, we can only expect to prove |--decay for the 
L^-energy of fc — 1 derivatives of curvature. 

The resolution of this difficulty relies strongly on the coupled character of 
the problem, i.e. the fact that the background is not fixed but related to the 
curvature components via the null-structure equations (cf. section 3.9). Working 
out the precise contractions of curvature components and components of the 
deformation tensor in (7) and inserting the structure equations relating Dtt 
to the curvature components, one realizes that all terms in (7) have a special 
structure. To give an example, 

^ p( (£^-1 [Daf + a + l.o.t.]) + (^C'f^ [divf + (3 + l.o.t.]) (8) 

is such a term, where f is an expression involving the Ricci rotation coefhcients. 
Remarkably, the derivatives of the deformation tensor have brought in precisely 
the right curvature components from the structure equations, so as to allow an 
integration by parts in T. This lets the T-derivative fall on p and makes this term 
a cubic error-term with all components decaying. The derivative term (which 
has Dsf and divf respectively) can be integrated by parts as well. Clearly, if 
the derivative falls on p, we obtain a lower order term. If it falls on the highest 
order derivative term, however, we can use the Bianchi equation 

aa = -2-0113 + l.o.t. (9) 

and obtain also a lower order term.^ At the end of the day, we can show that 
the worst possible error-terms are effectively of the form 

f p - (c!f^DnyC^W (10) 

J spacetime 

The gain of a derivative compared to (7) is essential to close the estimates. The 
analysis of this structure is the topic of section 7. 

1.1.4 The redshift 

As familiar from the wave equation, the vectorfield T does not provide control 
over all curvature components near the horizon, even if the metric is exactly 
Schwarzschild (in our case, T may not even be causal everywhere!). For the wave 
equation it is well-known how to stabilize the T estimate using a redshift vector- 
field [1]. In complete analogy to this case, we construct a timelike vectorfield N 
(agreeing with T far away from the horizon), which generates boundary-terms 
controlling all components globally, and a spacetime term Kf^^ [W] which 
has a good sign (and controls all curvature components) near the horizon. The 
error-terms introduced can be dealt with, once an integrated decay estimate in 
the interior (away from the horizon) is available. Cf. Proposition 6.4. 

Since until now we have only commuted with T, this technique will merely 
allow us to obtain non-degenerate control over arbitrary T-derivatives of the 
curvature components. This is sufficient away from the horizon in view of the 

^Noto that commutation of derivatives with Ct derivatives only introduces lower order 
terms and additional decay. 
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elliptic estimates (section 9) available. Close to the horizon we invoke a com- 
mutator version of the redshift estimate (developed in the context of the wave 
equation in [4]), for which one commutes with the redshift vectorfield to obtain 
estimates for the transversal derivatives on the horizon. For the Bianchi equa- 
tions, there are two (equivalent) possibilities to adapt this estimate: Either at 
the level of the tensorial Bianchi equations for the Weyl-tensor commuted with 
A^,^ or at the level of the null-Bianchi equations. In this paper, we follow the 
second approach, as it has the advantage that, at the lowest level, we can work 
with the renormalized null-Bianchi equations for p and a. Cf. section 8. 

1.1.5 Decay at infinity 

Very recently, Dafcrmos and Rodnianski introduced a novel approach to the 
decay of the wave equation 'Oi-ip = on black hole spacctimcs [10]. Using a 
new multiplier, they were able to generate a certain hierarchy of estimates for 
r-weighted energies near infinity provided an integrated decay estimate in the 
interior is available. A version of the pigeonhole principle finally allowed them 
to exploit that hierarchy to translate boundcdncss of r-weighted energies into 
interior decay in t of the natural L^-energy. 

This new method of obtaining decay at infinity and exporting it to the 
interior has an advantage over the traditional method of using the conformal 
Morawetz multiplier K = {t^ r^^ dt + Atrdr (which is conformally Killing in 
Minkowski space). Namely, as the latter carries i- weights, the error-terms it 
generates in the interior arc typically difficult to control. The multiplier of the 
new method, on the other hand, carries only r-weights and is only applied in 
a region close to null-infinity. This emphasizes the power of using multipliers 
locally, as one root of the difficulties with the /C- vectorfield is precisely that it 
has to be applied globally. 

It turns out that the approach of [10] has a very natural generalization to the 
Bianchi equations. The basic idea is easy to understand and can be explained 
at a heuristic level. Suppose we use a Minkowskian coordinate system {u,v,u;) 
near infinity, i.e. such that in particular 9„r « —di,r ~ — 5 there. In fact, for the 
purpose of the following heuristic, the reader can in fact consider the Bianchi 
equations on exact Minkowski space. In this region, the Bianchi equations for 
the components a and /3 take the form (cf. section 3.3) 

2dua--a = -2t)lp + lo.t. , 2dyP + -P = diva + lo.t. (11) 
r r 

Multiplying the first by r^a and using that ^2 is the adjoint of div on the 
2-spheres, we obtain 

du {rP\\a\\^) + 2a„ (rf ||/3||2) + {-p ■ 9„r - 1) 

+ \\PfrP^'^ (8 - 2pdyr) = error + tot. div on Sl^ . (12) 

Upon integration in a characteristic region (the shaded region in the figure 

^provided one also commuted with T first to eliminate the non-decaying p component 
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below) using the measure dudv sin 9d9d(f> 




the first two terms in (12) will generate positive future boundary terms. The 
remaining (spacetime) terms on the left hand side will both be positive as long 
as 2 < p < 8. In other words, provided we can control the error-term arising on 
the timclike boundary of the region (which can be done, provided an integrated 
decay estimate is available in the interior) and the error-terms on the right hand 
side (which are absent if the background was exactly Minkowskian and which 
cause a considerable amount of work in our case), we obtain an estimate for 
strongly r-weighted boundary and spacetime terms for the components a and 
/3. Considering the next pair of Bianchi equations, we will obtain a similar 
estimate for the pair (/3, p) with the condition 4 < p < 6. In this fashion we 
can estimate all curvature components with appropriate r- weights. From the 
weighted spacetime-terms we can generate decay in t in the interior using the 
pigeonhole principle as in [10]. Exploiting the hierarchy in the equation one 
easily obtains the well-known pointwise decay rates for the null-components in 
Minkowski space [11], for instance. 

As mentioned above, a significant amount of work goes into estimating the 
error-terms on the right hand side. They impose additional non-linear con- 
straints on the admissible r^- weights. In spirit, these estimates are of course 
similar to the ones in [8] , since only the asymptotic region and weights in r are 
concerned. We emphasize two important differences, however: On the one hand, 
we are not trying to minimize the number of derivatives. The cubic error-terms 
are always estimated by putting one term in and the others in (unlike 
estimates in [8] ) . Secondly, by using the null-structure equations directly on the 
error-terms, we exploit an additional cancellation of the terms which have the 
worst decay in r. This structure goes beyond the null form of the error-terms 
and simplifies considerably some of the estimates. The analysis at infinity is the 
content of section 10. 

1.1.6 Integrated decay 

We have seen that both the redshift and the new method of capturing the decay 
at infinity require an integrated decay estimate in the interior to control their 
error-terms.'' How does one obtain such an estimate? 

From the experience with the wave equation, one may try to use a vectorfield 
of the form X = f (r) dr in the identity (5). Carrying out the computation one 
observes that one cannot obtain a globally (positive, say) spacetime term, which 
controls all curvature components but that instead the quantities p and a will 
always enter with the opposite sign. More precisely, the main (assuming the 

'^A boundedness statement can in fact be obtained without the construction of an integrated 
decay estimate. See section 5. 
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metric is exactly Schwarzschild) spacetime term reads 



XTT 




(13) 



We immediately recognize the familiar trapping factor of (r — 3M) in the second 
line. Choosing / as a bounded function which changes sign at r = 3M we 
can control all components but p and cr in a region of bounded r around the 
degenerating set r = 3M.^ The boundary-terms arc in turn controlled by the 
T-energy. We remark that all these arguments need to be stabilized near the 
horizon by the redshift vectorfield (cf. section 11). 

The above observation regarding the role of p and a requires us to obtain 
an integrated decay estimate for p and a via different means. The important 
insight here, which goes back to Price (see [12] and [13]), is that the compo- 
nents p and a satisfy a wave equation (the so-called Regge- Wheeler equation), 
which in the case of exact Schwarzschild would be homogeneous and in our 
case involves inhomogeneous error-terms on the right hand side. Those terms 
exhibit a quadratic structure of (derivatives of) "Ricci coefficients • curvature 
components" . 

The natural energies for these wave equations do not involve p and a itself 
but rather the rescaled quantities r^p, r^a. On the one hand, this is good 
news because the quantity p ■ approaches the mass (a constant) for late times 
and hence we can expect decay for the L^-energy of its derivatives.^ On the 
other hand, this complicates the estimates for the error-terms appearing on the 
right hand side of the equation, as they now get multiplied with the weight of 
the rescaling factor, r^. In other words, one cannot prove an integrated decay 
estimate in the interior without understanding the r-weighted energies for (all) 
the curvature components! 

Here it turns out that the we can borrow an e of the strongly r-weighted 
spacetime integrals that we have available at infinity. Coupling the estimate at 
infinity with the anticipated integrated decay decay estimate in the interior will 
allow us to absorb the e-contribution and close the estimate. 

Besides controlling the inhomogcneity, an additional problem in the deriva- 
tion of energy estimates for the renormalized p and a components arises from 
the Regge- Wheeler operator itself. As is well-known for the Rcgge- Wheeler 
equation on a fixed Schwarzschild background, an everywhere positive space- 
time integral can only be derived if the angular momentum mode of the field 
satisfies I > 2. For us, however, the situation is more favorable because we 
are assuming that the metric approaches Schwarzschild and that we are esti- 
mating the curvature components of that particular spacetime. More precisely, 
the ultimately Schwarzschildean property implies that the zeroth order terms 
in the anticipated integrated decay estimate arc a-priori bounded, while the 

^Note that near infinity and near the horizon we obtain control over all curvature compo- 
nents with the given choice of /, since the terms in the second line eventually dominate the 
ones in the first. 

®We remark that the r-weight which is being introduced in this way corresponds precisely 
to the maximal r-weight that the quantity p gets associated via the r-weighted multiplier 
described in section 1.1.5 
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derivative terms can easily be made non-negative with the help of an appro- 
priate multiplier. In this way the integrated decay assumed in the ultimately 
Schwarzschildean assumption eliminates the I = and I = 1 modes (see section 
12). 

With these difficulties resolved one finally obtains an integrated decay esti- 
mate for the renornializcd p and a components, which is coupled to the inte- 
grated decay estimate for the entire curvature tensor and then, finally, to the 
r-weighted estimates at infinity. 

Having understood the basic mechanisms, we recommend the reader at this 
point to turn to section 5 immediately for a precise statement of the theorems. 
A brief glance at section 4 to familiarize with the energies which are being used 
may be helpful. 

1.2 Some remarks on the Kerr case 

It is natural to enquire which aspects of the mechanisms described stand a 
chance to (in principle) carry over to the case of a spacetime approaching a 
Kerr solution of, say, small angular momentum. In this case both the curva- 
ture components p and a will be non-zero, at least if a null-frame based on the 
principal null-directions of Kerr is being used (otherwise, other null-components 
will be non-zero, too). We do not expect any further difficulty entering from 
the non-decaying tr-component. The analysis of the intricate error-terms pro- 
portional to p in the commuted Bianchi equation carried out in section 7 will 
extend to the dual tj-component. The fact that T is no-longer (not even asymp- 
totically) timelike near the horizon, on the other hand, can ~ at least for the 
case of small angular momentum - be dealt with the redshift, cf. [4]. This leaves 
the trapping as the main difficulty. In view of the trapping for Kerr itself be- 
ing well-understood, a generalized approach of Dafcrmos and Rodnianski [7] or 
Andersson-Blue [6] (using commutation with the second order Carter-operator, 
i.e. an ultimately Killing tensor in our language) is highly promising to yield 
similar results to the ones in the paper. We postpone the analysis of this case 
to a future paper. 

2 Ultimately Schwarzschildean spacetimes 
2.1 The Schwarzschild geometry 

In the perhaps most familiar (t, r, 9, (j)) coordinates, the Schwarzschild metric is 
given by 

g = -{l- ^) df +1^- ' dr^ + {dO^ + sin^ Odcf") , (14) 

which defines a smooth metric on the manifold Aii = (—00,00) x (2Af, 00) x 
. Using the the transformations u = t — r* , v = t + r* , where r* = r — 
2 A/ log (r — 2M) — M + logSM is the tortoise coordinate, we may express the 
metric on Aie in double null-coordinates (u, v, 0, (p) G (—00, 00) x (—00, 00) x 5'^ 
as 

g^-{l- p) dudv + {d9^ + sin^ edc/}^) , (15) 
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where we set fi = It is well known that the singularity at r = 2M in (14) is 
a coordinate singularity and that the metric can be extended through r = 2M 
to a larger manifold. Indeed, let x i''^) be an interpolating function which is 
equal to 1 for r < 6M and equal to for r > 7M. Setting 



t*=t + f{r) with /'(r)=x(r) 



2M 
r - 2M 



we obtain 



g = -fc_ {dt* f + —xdt*dr'- 
r 



dr 



2 I 2 



{dO 



2 ' alv? ed<i)'^) , 



where we have introduced the shorthand notation k± = 1 ± /i. Clearly, in 
the coordinates {t* ,r,9, (j)), g can be defined on the larger manifold Aiji ~ 
(— oo, oo) X (0, oo) X S'^. We easily see that r = 2M is a regular nuU-hypcrsurface, 
which we denote "H"*" and call the future event horizon. 

Let us fix on Mji = (— oo, oo) x (0, oo) x S"^ a shce of constant t* = tq > 1 
(note that this crosses "H^) and define the submanifold 



7^ = [to,oo) X [2M,oo) X S"^ 



For convenience, let us also define 



2M 
1± X 

r 



!3rr 



x^k+ + ^ (1 - x') 



(16) 



(17) 



Note that fc^ is bounded above and below by positive constants depending 
on M only. 

Next we define various frames on M.11. The vectors 



R 



ei_2 — frame on S'j* ^ , 



z9r and n = sfkZdt* 



2Mx 



dr 



(18) 
(19) 



satisfy 

g {R, R) = —g {n, n) = +1 as well as g {R, n) — g (T, ca) = g {R, ba) — , 

and hence constitute an orthonormal frame. In addition, we may choose a null 
frame, i.e. vectors ei, 62, 63, 64 such that 



3(63,64) = -2 and 5 (63, 63) = g (64, 64) = g (63, e^) = 5 (64, 6a) = , 
as well as 



9 (6A, 6b) = Sab , 
for A, B = 1, 2. Such a frame may be defined via ei, 62 

2Af . 



63 = ^= (n - i?) = - — 



64 



dr , 



^ (n + i?) = fc^a** + h - ^x ] dr 



(20) 

(21) 
(22) 
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For Schwarzschild, the Ricci coefficients with respect to the null-frame (ei , 62, 63, 
are - using the conventions of [8] summarized in appendix A: 



i±!l^ ^-^ 



5ai 



Hab =1 X -Oab 

r r 



(23) 

(24) 
(25) 

Ya^^a = Za = Zj, = Va=^ . (26) 

Finally, we collect some useful formulae for further reference: The timelike 
Killing field expressed in the chosen null-frame reads 



O = 



M (1 - x) + x'^fc- 



M 



7 (x - fx) 



2M \ 1 / 2M , „ 
2T = 2dt. = 1 X 63 + ^ 1 + X 64 . 



The normal of constant t* slices Et* is computed to be 



n 

and hence the components 



/fc^a*. = T^Vh^a + 9^7=64 



n 



(27) 



(28) 



(29) 



2.2 The class of ultimately Schwarzschildean spacetimes 

Fix a Schwarzschild spacetime of mass M, g^^) , equipped with the regular 
coordinate atlas {t*,r,Oi,(j)i) as defined in section 2.1. The Penrose diagram of 
(TZ,gM) is depicted below, where we use the standard notation 1+ to denote 
future null-infinity. 




r = ry 



R 



Let us write St* for slices of constant t* . Slices with "null-ears" are denoted 
St*: More precisely, if S^t ^ denotes the sphere of intersection of St* and a 
fixed^'' r = R hypersurface, we define 

Nout {St r) ■= future outgoing part of the lightcone emanating from r 



^"The exact R will be fixed later in the paper, with figuring as a potential source of 
smallness in the argument. 
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St* := (St* U {r < i?}) n Nout (5^1, a) ■ (30) 

The spacetime slab enclosed by E^-i , and the horizon is denoted by M (ti , r2), 
while the spacetime slab enclosed by St-j , St-^ and the horizon is denoted Ai (ji , T2 ) 
We also define the asymptotic region V (ri,T2) = M (ti,T2) n {r > i?}. Close 
to the horizon we also fix a hypersurface r = ry, whose precise location will be 
given in Lemma 6.3. 

We will now fix the differentiable structure of TZ and consider small pertur- 
bations of the metric on the black hole exterior to the future of Eq. As the 
differentiable structure is fixed, 5^* ^, S,-, S,-, (ri,r2), J^{ti,T2) ,etc. will 
remain well-defined with respect to the perturbed metric, and we will hence 
use the same notation for them, as long as no confusion arises. We denote the 
covariant derivative of the perturbed metric g hy D. 

Definition 2.1. Fix the manifold TZ and its differentiable structure. A metric 
g on TZ is called ultimately Schwarzschildean (of mass M) to order k -|- 1 

(k > 2) if it has the following properties: 

1. The boundary is null with respect to g. 

2. The metric g is C^~^ -close to the Schwarzschild metric in that 

l5.-(/0,,l + l.9'^-(./0"l<;i (31) 

\^l9^, - (5^0 I + 1^5'^^' - dl, (/O"' I < 4 (32) 
holds for n — 1, fc — 1 and m G {t*, r, 6, 0}. 

3. There exists an admissible null-frame ei, 64 for g, i.e. a frame satisfying 

(a) closesness to the Schwarzschild frame^^ e^.' 

\g (eg, 64) + 2| + |,g (64, eg) + 2| + \g (e^, cb) - Sab\ < e (33) 

while all other combinations are e-small in the above sense. 

(b) The outgoing null-hypersurfaces generated by 64 foliate TZ. More pre- 
cisely, there is an optical function u whose level surfaces are the 
outgoing null-hypersurfaces generated by 64. This optical function 
is normalized such that u = 00 is associated with the horizon . 
Moreover, the ei, 62 are tangent to the spheres Sf* „ arising from the 
intersection of the timelike t* = const slices and the outgoing null- 
hypersurfaces u ~ const. 

(c) the gauge condition 2Y = g {D^ca, 64) = 

(d) the Ricci coefficients are ultimately Schwarzschildean to order k [see 
Definition 2.2] in this frame 

^^In view of the fixed differentiable structure, the (defined in terms of coordinate vector- 
fields in (21)) remain well-defined with respect to g. 
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I. There exist functions p : 7?. — > M and g : 7?. — > M such that p vanishes on 
the horizon and that the associated vectorfield 2T = pe^ + qe^ is ultimately 
Killing to order k + 1 [see Definition 2.3]. Moreover, the deformation 
tensor satisfies tr^'^^i: = and 



\T{p) I + \eA ip) \ + \Tiq)\ + \eA (q) \ + |03 [p-k. 



X - 



-104 [p -K]\ + \iP 











1+04 





|<e(0-^ (34) 



holds in the interior region r < t* . 
Remarks. 

• Our assumptions are certainly not minimal and leave room for improve- 
ment. In particular, the decay properties on the deformation tensor in 
assumption (4) may be derivable from (3d). 

• The assumption that the frame satisfies Y — (which corresponds to a 
partial choice of gauge) will be of importance in the analysis at infinity, 
cf. [8]. It can easily be achieved by Fermi-propagating ei,e2 in the 64 
direction. 

• The assumption that T can be chosen to be null on the horizon (i.e. the 
condition — 0) is not essential for the results of the paper. However, 
it allows an elegant proof of uniform boundedness of the highest order 
curvature norm. We will show later how this assumption can be removed 
in the context of integrated decay estimates. 

• The assumption that tr'^'^^ir = is a technical assumption which allows 
one to reduce the number of error-terms in the proof at some point. The 
additional errors could be controlled without further difficulty but would 
make the paper longer (cf. section 12.5, Remark 12.19, where this is dis- 
cussed). 

We will write P = P (t*, u) for the projection from the tangent space of [TZ, g) 
to the tangent space of S"^* ,^. If D denotes the covariant derivative on {TZ,g), 
then we denote by 03, 04 the projections to Sf* „ of D3 and D4. The covariant 
derivative intrinsic to 5*^^* „ will be abbreviated V. We will write I? shorthand 
for any derivative 2? £ {03i 04 7 ^1- 



2.2.1 The Ricci coefficients 

We denote the Ricci coefficients in an admissible null-frame collectively by 
While 9\ contains all components, we generally distinguish 

$H™"" = {H,Y,Z,Z,V,n,n,trH,trH} (35) 

and _ff , which has to be treated separately, as it contains information from the 
radiation field in its ^-term. Also, 5Hg™"™'' denotes the corresponding Schwarz- 
schild values. Defining the norm of the null-decomposed Ricci coefficients in 

^^That is the values of the components of the Ricci coefficients of (TZ,g^^) evaluated with 
respect to the frame e^. 
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the obvious way, i.e. 

||r^B||= E \T{eA.eB)\ and ||[/^|| - ^ |[/ (e^) |, (36) 



A,B=OA 



A=0,l 



we can state 



Definition 2.2. The Ricci coefficients of g are ultimately Schwarzschildean to 
order k with respect to the null-frame (61,62,63,64) if the following estimates 
hold: 



globally on the black hole exterior 



2 / 2M 
trH - - 1 X 



trH 



2 / 1- 



2M 



< 



\\Hab\\<— and \\Hj^b\\ < TTT* T 

r^ r max {\,t* — r) 

\^~^ss\ + U-Q.ss\ < 4 



||y^|| + ||ZA|| + ||z^||<- 



(37) 

(38) 

(39) 
(40) 



• in the interior region, r < t* , the pointwise estimate 

p-SHssll < e(t*)"^ (41) 
i.e. (37)-(40) hold in r < t* with all right hand sides replaced by e (i*)~* . 

• The estimate 



[Dl] (r) B" [<n] (t, 00) < 



holds for n = 1, ...,k and 
l{n) 



if n = k 

1 if n = k — \ 

2 ifn = k-2 
I ifn<k-2, 



(42) 



(43) 



where 



r [m] (ri, T2) = sup E [5H] ( ) + E [m] {H (n, ra)) 

-re(Ti,T2) 



+f [5n] (7W(Ti,r2)) (44) 



with the weighted energies for the Ricci rotation coefficients (5 ~ -^) 



n „ 



dt*drdujr^-\'^'''+^''-' 



'l(n) 
,3=0 



rnm(2fci-i+l,0)^i-'(n) 



(45) 
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E 



dvdLor^-'r''^-+^'''\\$l'$l'i/'''Hf (46) 



as well as 



n „ 

r [91] (7^ (n, T2)) = ^ / IIP' (9^ - 



(47) 



i=0 "''"(^I'-^a) 

• /n addition, for the boundary term of the highest derivative 

sup l'' [9^] (Sr) + e''' [9^] CH (ri, ra)) < e • (n)"^ (48) 

and in the asymptotic region 

f [9^](P(ri,T2)) <e.(Ti)-^ . (49) 

Remarks: 

• Clearly, taking a 4- or an angular derivative improves the r- weight in the 
L^-energy by a power of 2: this is the factor j-'^^^+'ik-i ^ 

• By convention, the j in X]j=o runs 0, 1, 2, |. This sum incorporates that 
in general, one loses an r-weight for every power in t gained. However, if 
a three-derivative is applied (fci > 0), decay in r is gained allowing one to 
extract decay in t without losing powers of r. This is the reason for the 
factor min {2ki — j + 1,0). The component H_ is special as it decays less 
in r near infinity (cf. (38)). 

• The additional assumptions (48) and (49) are "justified" by the fact that 
we are going to show decay for the boundary terms and the degenerate 
spacetime energy norm squared of curvature at this level of derivatives, 
which "improves" this assumption. We remark that assumption (49) is 
technical and can be eliminated with some further work. It simplifies the 
error-estimates near infinity. 

• the precise reason for the various weights will become much more clear in 
conjunction with the curvature-energies and our method to prove decay 
for the latter. Note however at this point that the assumptions on k 
derivatives of the Ricci-coefRcients immediately imply estimates on fc — 1 
derivatives of curvature, which are easily inferred from the null-structure 
equations of section (3.9). 



^We write r'^dvduj for the induced (with respeet to 64) volume form on Nout {s'^ 
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Before we proceed, we introduce two more energies which will be useful. 
First an energy on the untilded slices 



n „ 



(50) 



and second, a A-weighted energy on the Ricci-coefficients: 

[m] (ti,T2) = supl" [$H] + Bx ■ D"-i [m] (ri,T2) 

+X-r[y\](^M{n,T2)) +r[m]{V{n,r2)) (51) 

and B" [9^] (r) [91] (r, od). Note that this energy only requires a A-small 

contribution of non-dcgcnerate integrated decay for n derivatives of the Ricci 
coefficients near the photon sphere, since the last term regards only the asymp- 
totic region. In particular, we see that if {A4,g) is ultimately Schwarzschildean 
to order k + 1, then [$H] (t) is e- A-small for sufficiently large r. This definition 
is again tailored to what we are going to show for the curvature, cf. the remarks 
at the end of section 5. 

2.2.2 The timelike ultimately Killing field 

We recall the definition of the deformation tensor of a vectorfield X. It is the 
covariant two-tensor 

(^V {y, Z)^^{g (Z, DyX) + g (y, DzX)) (52) 

whose norm is computed in the orthonormal frame: 

||7r||:= k(e-e,)|. (53) 

i,i=l,..4 

We also introduce the null-decomposition of its traceless part 

(^)m^ = W7f4^ , Wm^ = W^3B, 

Wn=W^44 , (^)n=W^33. (54) 

Definition 2.3. The vectorfield T is called ultimately Killing to order k -I- 1 

if its deformation tensor satisfies 

• on the entire black hole exterior the pointwise estimate 
II '•"^V except ^'^'^iAsW < -j- as well as W^'^^'iAsW < 



r ■ max (1, — r) 
• in the interior region, r < t* , the pointwise estimate 

\\^^K\\<e{t*)-i (55) 
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• For any tq < ti < T2 



mm(2fei-i+2,l)^J-'(") ^2k2+2k3 



J2 I dt*drdLur^-^ i^r ,^ 

fci+fc2+A;3=»-^^(^i'^^) \j=0 

'''1 

where i ~ 1, fc and I (n) as in (43). 

3 The Bianchi equations and the Bel-Robinson 
identity 

In this section we collect the main energy identities and commutation formulae 
that we are going to need in the paper. We use the conventions of [8] and refer 
the reader to this reference for detailed derivations. 

3.1 The energy identity 

As we are going to work with higher order energies arising from commuting the 
Bianchi equation with various vectorfields, it is convenient to study directly the 
more general inhomogeneous Bianchi equations 

D'^Wcp-tS = Jm& , (56) 

for any tensor W which has the algebraic properties of a Weyl-tensor. Defining 
the left and right duals 

*Wa/37A- = ^ eafif^u W^^^^ and W*aP',5 = ^^0/3'"' ^^^^5 

with (zap-yS denoting the volume form of g, we can also write (56) as 

D"*Wa.PyS = Jp^s ■ (58) 
We define the symmetric, traceless Bel-Robinson tensor 

Q Ma^PyS = W.p^.W// + *Wap7.*W/,'^ . (59) 

It satisfies the identity (5), from which we derive 



/ Q[W]iX,y,Z,n^)d^l^+ f Q[W]iX,y,Z,nn)d,^n 
- [ {K^y^ [W] + [W]) d^,= f Q [W] {X, y, Z, ns) dfi^ , 



'M{Ti,r2) 

(60) 

for vectorfields X,y,Z, where 

K^y^ [w] = g"^^^ [^''^TTc.py^Zs + ^^\^pz^xs + ^^\c.px^ys) m 
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K^y^ [W] = {divQ)p^s X^yz' (62) 

where 

{divQ)f,^, = W// J-^,. + W// J-^^. + J-,:,,, + J';^, . (63) 

An analogous identity holds for the tildcd slices S and M.. In that case, there 
is an additional boundary term at null-infinity. 

Regarding the boundary term, we note [8] that Q [W] {X,y,Z,n) is non- 
negative for any set of causal vectorfields X^y,Z,n. It will sometimes be 
convenient to evaluate this boundary-term for the Schwarzschildean normal 
ns = ^fc^e3 + 5-^e4, which is C^-close to ns by the ultimately Schwarzschildean 
assumption and use the estimate 

IQ [W] (A-, Z, ns) - m [X, y, Z, ns) I < m {N, N, N, N) , (64) 

where is a timelike vectorfield with b < g {N, N) < B for constants 6, B 
depending on M only. 

For future applications, we finally collect the null-components of the defor- 



mation tensor of a vectorfield of the form 

X = «8e3 + 1)64 with <B and D C^-bounded functions . (65) 
They read (writing the short-hand tt for ^'^^tt) 

TT^"* = i [-63 (») - 64 (D) + -f 22)17] , (66) 
7r33 = 1 [-264 (*B) - 4Q5f7] , tt^'' = J [-2e3 (D) - 4Dfl] , (67) 

TT^^ = -i [-6A (2)) + QSr^ + 2)Z^ + DV^] , (69) 

TT^^ = - i [-6A (*B) - "BV^ + SZ^ + Sy^] . (70) 



3.2 The (tensorial) Bianchi equations under commutation 

We now state the commutation formulae that arise from commuting the Bianchi 
equations with vectorfields. For this we need the notion of a modified Lie- 
derivative (cf. [8]), which ensures that CxWap-yS is again a Weyl field. 

Definition 3.1. The modified Lie derivative of a Weyl-tensor W is defined as 

CxW CxW -l^^^[W] + ltA^^nW (71) 
2 8 

with 
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Proposition 3.2. Let 

D''W^0^s=jp^s (73) 
and X an arbitrary vectorfield. Then 



+ [Ds^c^x - Dx^c^s) ^ j = CxJM& + Jm& {X, W) . (74) 

Here the last equality is to be understood as a definition of Jp^s . 
Proof. This is Proposition 7.1.2 of [8]. □ 
Corollary 3.3. 

Ti-l 

(c^CtW) = d^Jp^s (T, W) + y] &tJpiS (t, CI^^'-^CtW) . (75) 

' fc=0 

3.3 Bianchi equations for the null-curvature components 

For W a Weyl tensor we define its null-eomponents 

= KK^py-s {esr (63)' , (76) 
a^, = nP^F^Wp^aS (e^y {aY , (77) 
= ^P^WP7-^ (^3)" {esr (e4)' , (78) 
= iK^py^^ (eif ie^r (64)' , (79) 
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P = ^Wp^.S (63)" (64)" (63)^ (64)' , (80) 

a - ^e^^Wp^-i (63)" (64)' , (81) 

with the projection onto the 2-spheres ^. A computation reveals 

Lemma 3.4. For Schwarzschild and its null-frame 60,61,62,63 all curvature 
components are zero except 

2M , , 

p = -—=:-p. (82) 

We will consider certain derivatives of these null-components in the null di- 
rections. The following definitions are useful for a (differentiated) null-component 
u of signature s and rank k (cf. [8]): 

u, = lP,u^^-^^±^trEu, 

^,^0^^+^±l(4±AMi,,/,. (83) 
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Here, the signature is defined recursively as follows. For a p-tupel {p > 0) Up of 
3's and 4's (in any order) 

s (ckrip) = 2 + (number of 4's - number of 3's in rip) , 
s (/3rip) = 1 + (number of 4's - number of 3's in Up) , 
s (prip, Crip) = + (number of 4's - number of 3's in Up) , 

s (^P^ ^ = —1 + (number of 4's - number of 3's in Up) , 
s ^f^rip) = ^2 + (number of 4's - number of 3's in Tip) . 

The rank k is simply the length of the tuple, e.g. k [ctup) ~ P- Note that 
s (q;4„p) = s (ctrip) + 1 and s (a3„p) = s {oiup) — 1- We also define the short- 
hand notation 

^± M = ^^ "^ . (84) 

Following [8], we use the following notation for operators on the spheres Sf* „: 

• takes any S'^-tangent 1-form ^ into the pair of functions (d/w^, ci/rl^) 

• 'p*i, the L^-adjoint of 'Pi, takes any pair of scalars p, a into the ^^-tangent 
1-form -^aP + f^Asf cr- 

• ^2 takes any 2-covariant symmetric traceless tensor ^ into the S'^-tangent 
1-form d/u^. 

• ^2 1 the adjoint of ^2 takes any 5'^-tangcnt 1-form ^ into the 2-form 

-5 {"^bU + fA^B - [dH) $ASj 

Finally, we decompose the Bianchi equations as equations for the null-curvature 
components (sec Proposition 7.1 for the definition of (J) ,6 (J)): 

ag = + \4.rH_OL = -2Tt)\p, + E-^ (a) (85) 
E-i (a) = ma - 3 [up + *Ra \ {V ^ AZ) ®f3 - 26 (J) 



/34 = $4^ + 2trHI3 = diva + E4 (/?) (86) 
{13) = -217/3 + (2F + Z) • a + 3rp + 3 {*Ya) - J4A4 



"4 = 04" + ^^rif a = +2'0ll3_ + Ei (a) (87) 
E4 {a) = 4r2a - 3Hp - 3 f-*^cr) + (V' - 4Z) §/3 - 26 (J) 



^3 = -^3^ + '^t'^K^ = -diva + E3 (P) 
E3 {p) = -2^1 +{2V - Z)-a-iYj) + Z {*Y_(j) + J3A3 
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/33 = + trHj3 = +Tt)^ {-p, a) + E^^ (/?) (89) 
E^iiP) ^ + 2H ■ P + a + iZ p + 'i {*Z(j) + .Haa 



Ea 



= IPip + trHp_ = +1^; [p, a) + Ei (/?) 
ip) = +2n/3 + 2H-P-Y-a-3{Zp - *Za) - J4A3 



(90) 



P4 = 04/9 + -trF/j = d/w/? + E4 {p) 



Ei (p) = -^1 •a + V^-/3 + 2(Z-/3-y-^)- ^ J434 



(91) 



(T4 = 040- + -trHa = -ci/rip + E4 (cr) 

Ei (ct) = +^1 • - y • -^z? - 2 (z • */3 - r • *^) - i j:34 



(92) 



P3 = 03/0 + l^trHp = -divP_ + E3 (p) 
Esip) = -^H ■ a + V ■ P + 2{Y- p ~ Z ■ p) + ^J334 



(93) 



(T3 = 03(7 + -trHa = -ci/rip_+ E3 (a) 
E3 {<J) = -\h -^a + V ■*P_+2{Y-*P + Z-^P) + ij3*34 



(94) 



Note that at the lowest order we are going to renormahze p = p + 2M and 
a = r^a OT ahernatively p = for which we expect decay. Hence we also 
collect the equations 



P3 



2M 



-divP + S3 {p) - 



6M flp 



Pi 



3 

2M 



1 



d^vP + E.ip)-'^^^^' 



-trH 



1 

-trH 
2 



-divP + Esip) , 

+ i?4 (p) . 
(95) 



3.4 Commuting the null-Bianchi equations 

Besides the important commutation formula of section 3.2 we arc also going to 
commute the Bianchi equations for the nuU-curvaturc components directly. 

From the first order Bianchi equations of the previous section we derive the 
higher order Bianchi equations inductively using the following Lemma, which 
we quote from [8] (Lemma 7.3.3). 

Lemma 3.5. Let Uai...a^. be an S'^ - tangent k-covariant tensor on (JZ^g). Then 

Jj)^^ bUa^...A^ - 1/BlPiUA,...A^ + HBcfcUA,...A, = FiBA,...A, , (96) 
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I1)^11bUa,...A^ - 1/b1P3Ua^...A, +HscfcUA,...A, = F3BA,...A, , (97) 



ll)^U)iUA,...A, - J1)SzUa,...a, = , (98) 

where 



F3BA^...A, = YB0iUA^...A, + [Z B - Vb) 1/)3Ua,...A, + ^ (Ha,bY. 



-C 



i=l 



-HbcILa, + Ka^B^C - KBC^Ai- &A,C * P_iR) b jUAi...A,-iCAi+iAk 7 

k 

FiBM...A, = Yb$3Um...a, + (Ab + Vb) $iUM...A, + ^ (Sa.b>^c 

i=l 

-KBC'YAi + Ha,bZ_q - HbcZ_a^+ ^AiC *P {R)b ) UAi...A,.iCA,+iAk , 



F3iA,...A, = -2n$3UA,...A, + 2n$^UA,...A, + {Zb - Zb) fBUA,...A,+ 

k 

2 II (YaXc - YcY^. +Za^Zc - Za,Zc + eA,c<J {R))Ua,...a^^,ca,^,a, • 



Remark 3.6. There is a typo in [8] regarding the sign of the first term of the 
expression for ^34^^...^^. 

We recall the definition (83) and declare that for fti a basis of angular mo- 
mentum vectorfields in Schwarzschild, the signature of jt-n-Un^, is defined to be 
identical to the signature of Un^, ■ 

Proposition 3.7. Let ni be an l-tuple of 3's and 4 's, I > 0. The Bianchi 
equations imply the following equations for higher derivatives: For \ a \ 

a3n, - -2-plPn, + S3"' (a) (99) 
with the error defined recursively as 
El^"" {a) = 03£;3"- (a) + (a3„J trHE^^ (a) 

- d- (/?„,) {ftrH) §/3„^ - Mvpn, + *Hc4riPn^ + 2s {F3 (/3„J) 
E;;^" (a) = 04£;3"'' (a) + i?+ (a3„J trHE^^'' (a) 

~~ ^+ (/3„J {I'trH) §/3„^ - HdivPn, + *HcilriPn, + 2s (F4 (/3„J) 

and also 

/^>3„, = -21)1 [iZAn] + («) (100) 

£;3-"" (a) = (Ef'^'^ (a)) + 2 [fl /I^'/^n, 



For 



/33«, - ^Pn, + + ^3 ' (/3) (101) 

/34n, = d^a„, + El' (/?) (102) 
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i?3"^' (/3) - $:,E;^ (/?) + (/33nJ trHE"^^ + F3 (p„J + 

- ??- (p„J PnJtrH - d- (p„J cTn^*ftrH - H ■ ^p,,^ - • fa^^ 
El"^ ifi) = U)^E;^ iP) + d+ (/33„J trHE';^ (/3) + F4 (p„J + ^^^^4 (^„J 

- 7?+ (p„J Pn,ftrH - 79+ ((T„J a^/ftrH - H ■ ^ p^^ - *H ■ fa„^ 
El-^ (13) = $^E2- + d- (/34„J tvHEl- (/3) 

- (an J • a„^ - H_ ■ fa,,^ + trF^ (a„J 
E^-" {P) = Ip^E'l- + d+ (/34„J trHEl^ (/3) 

- 79+ (a„J 1/trH ■ a„„ - H ■ fa,,^ + trFi (a„J (103) 
and in addition 

InAn, = nZPr. +*flZ-n, +E-,'''^ 

/n. An, = dMn, an, + £;r''''" (104) 

i^a'''^" (/?) = /o. (/?)) - /or' (P.'-)n, 

i?r'"" (/3) = In, {eT""^'' (/?)) - [d^, In,] (105) 

for 1^ 

P3„, = + (p) (106) 

P4n, = divf3^, + E2' (p) (107) 

ip) = IP^E'^' (P) + (P3nJ irHi?^'- ip) - trFs (^„J 
+ (^ )(5 ■ ftrH + H-'^I3 

\!—np J —rip —rip 

Es"' (p) = IPiE^- ip) + ^+ (p3nj trHE'i" (p) - trF^ (^^^J 

El-^ {p) = Ip^El- ip) + (p4„J trHE^- ip) + trFs (/3„J 

- ^" (/3nJ /3np • ftrH - H ■ fprip 

E2''^ [p) = IP.Er ip) + ^+ (P4nj trHE2^ ip) - trF, (/3„J 

- 79+ (/?„J /3„^ • yirH - H ■ fPr^p (108) 

and in addition 

/™P3n,--«0.^„,+C''^ (P) 

/™ P4n, = divlnAu + E^'""^' ip) (109) 

i?3-"^ (p) = tn. {e-^''^'' ip)) - tz:"L, 

Ef'^ ip) = /o. {Ef'^'' ip)) [div, ^n,] (110) 
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For [a] 



crsm = -ct/rlB + E^' (ct) 



(111) 

(112) 



E;^^' (a) = Ip,E^^ (a) + z?- (ct3„J trHE^^ (a) - trF^ (*^^J 
-19^(0 V/3 • ytrH - H ■ f*B 

\— jip / —rip —Tip 

E';-" {a) = Ip^E'^-^ (a) + {a^np) trHE'^" {a) - trF^ 

\^ lop J 1 tp ftp 

i;^"-^' (a) = 03^2" (a) + z?- (a4„p) tr^S? (a) - trFg ('^/3„J 

+ t?+ (/3„p) '^/3„p • ^<rS + K ■ ^*/3„p 
ii;^-^ (a) = 04i?4"* (a) + (a4„p) trffi?:'' (a) - trF^ (*/3„J 

+ ,9+ (/3„J */3„^ • ytrif + 7? • T^up 



(113) 



^n.cTsn, = -cilrlQfi + (cr) 



c4rlt:i{B.n,^El'''^ (a) 



(114) 



for 



3 



(a) = (i?,;""'"" (a)) - [-c^W, 

(a)) - [-c^/W,^J/:,7'a 



Til £1^ 

■4 



^3„, =-d/««„,+£^3"' (^) 



(115) 



(116) 
(117) 



iSg"''^ {P) - ^giSg"" (^) + ^" (£^^^J tr^i?3"'' (^) - trFa (a„p) 

+ (a„p) i^trH ■ + 1 • fa^^ 
El-' (^) = Ip.El- (^) + ^+ (^3^^J trHEl- (^) - irF4 («„^ 



.n„3 



trHE' 



F3 (Pnp) - (c^np) 



iJ^"* (^) = 04i?4"'' (^) + (^„J ^^^^^4"" {§) + {Pnp) Pupl/trH 

- {(Jnp) 0.ap*1/trH - H ■ yp„^ - *H ■ l/Unp 
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-3ni 



E 



(118) 
(119) 



For 



E. 



rii 



f2i \ 



71 1 n 



m-1 , s 



(120) 
(121) 

(122) 



Et' (a) = 03i?4"' (a) + 1?- (a4„,) tri7£;4"' {a) 
Et' («) = IPiEl^ (a) + (a3„J tvHEl^ (a) 



Ei 



-21^* (/a,^„,) + E. 
'n. {eT""^^'' (a) 



4 



(123) 



3.5 More commutation formulae 



The commuted Bianchi equations in the last section are not yet of the form 
most useful for the multiplier estimates. The point is that we would like to 
write the left hand side of (99) as (am)^, so that the equation becomes an 
honest inhomogeneous Bianchi equation for the components a„, and /3„, . Of 
course, this simply means to pushing the 3 trough the ni. 

Let u be any curvature component, ni a tuple of 3's and 4's and a, 6 e {3, 4}. 
We define the commutator 

Cab [Uni] = Uabni " Uan,6 ■ (124) 

Note that wc can write both 

U3n,4: = Usini - C34 [u„J and lt3„,4 = M43,i, - C34 [m„J + [U34 - M43]„, ■ 

Let us agree on the following convention: If the Bianchi equation for U3„,4 
is derived using U34„, of the previous section, the resulting error-term will be 
denoted with a tilde. On the other hand, if the equation for U3„,4 is derived 
using the Bianchi equation for U43„, , the error-term will acquire a hat. For the 
Bianchi equations for W3n,3 and W4„j4 there is no ambiguity, as they arc derived 
using W33„, and M44„, respectively. In this case, the error-term is also denoted 
with a tilde. For instance, for /3: 



■3ni 



Ef^^ {P)=El-^ (^)-C, 
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—rii 



(125) 
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Ef 



--34 



, (126) 



K iP,<^)su,+Er' 



--43 



(127) 



E. 



On the other hand, we could write 

t\{P^^),n,+ET' 



4"! I R\ — J^^ni 



--44 



/3 

—ni 



(128) 



-3n,4 



El""' 



-c. 
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/3 



(129) 



^„,3 = -d^a4„,+4'"' (^) 



4ni 



4ni 



^3-^34 



— C43 



/3. 



(130) 



The formulae for the other components are easily derived. We also collect the 
following formulae 



a43ni — aurii + 



a [i?" (a) $itrH - d+ (a) Ip^trH] + F34 (a) (131) 



/343n, — Pmui + 



P [d+ (/3) $^trH - {)- ip) $itrH] + F34 (/3) (132) 



P43n! — P34n, 



p [d+ (p) Ip^trH - d- {p) Ip^trH] + i^34 (p) (133) 



CT43rii — Cr34,i, 



a (a) Ip^itrH - [a) IPitrH] + F34 (a) (134) 



^34n =^43„ + ^H^(^)^3irH + ^9-(^)^4t'-S] +i^43(^) • (135) 



3.6 Commutation with angular momentum 

We can commute the Bianchi equations for the quantities u^mi etc. derived in 
the previous section with the fij. In this way we obtain, for m > (using the 
convention that E'.^'^^ (^u) = (u) and (y) = i?;"' (u)) the formulae 



(136) 



^3 



^2' *-J2i 



[03, /oj (/n' 'a4„, ) - (a4„J (^i.tr^) /^"^ 'a^^, (137) 



-Oi P4ti, 
m — 1 
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For /3: 



E, 



{tZ^Pin), = (/o.a4„,) + ^r'""^ (/?) 

(/^r V4„,) - ^+ (/34„, ) (f^^tr-T?) /o7'/34„, (138) 



^3 



(139) 



?n— 1 -iTTi— 1 



(/3) = (4^-""" (/?)) + 

- [03, /nj (/nr'/34n,) - ^" (/34„,) (f^»tr^)/I?,"V4„, 



For p and a 



(/r.>4,.)4='^^-(^'^Z/34«,)+4'"'''^ (p) 

" (p) = (Sf '"""^ (p)) + (ro-V4„,) 



(/n.cr4„,)4 = -cilrl (/^./34„J + i:^"'"' (cr) 



E, 



' (a) = (^^"'""'^ (a)) + [-c^rUoJ (/rV4„, 
[04, /oj (/JZ" - d+ (a4„, ) (fi.iri/) V4„, 



-^^^n,l -rEp''' (p) 



[03, AiJ (/orV4n, ) - (P4„,) (r!,trH) /,7-V4. 



4«,73 

(p) = (4^-"'" ' (p)j + [-d/.,/oJ (/rV^n, 



,)3 = -cVw(/:?AJ + 



i^g'"'''^ (a) 



4ni S7 



" (a) = /a. (4'"''^" (-)) + [-c^^UnJ (/^"'^„, 

[03, (/rr'^4r J - (a4„J(17,tr^)/^fV4. 



For /3 



p4n,n™ 



- [04, /n.] V4„, - ^+ ^4,. (^»*'^^) ^r>4„, (146) 



(140) 



(141) 
(142) 

(143) 
(144) 

(145) 
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(147) 



n.i^inj + E3 



(148) 
(149) 



E3 



(150) 
(151) 



[03, /oj (/r^^3,J - ^- (^-.^^i^) ^r^^3„, (152) 



For a 



{^n.a,J,^2TPl + i^f'^^'" (a) 



Ei (a) 



A^Tl— 1 



[04, /oj Un, - («4,J mrH) 'a^^^ (153) 



-ni^4nJ3 



(^n,^3. 



^4 (^) 



(154) 



" («) = /^i, (^f '"-""^ (a)) + 2 (/r^s 

[03, (/fr'£^4„,) («4n,) (^^.i'^^)/rr'^i4n, 



(155) 



n:«3,J4 = 2^2(/IZ^3,J+^4"'''^'fe) 



(156) 



^4'"'"'' («) = /o. ( i^I'"""' (a) ) + 2 1^*, (/^;r V, 



, rn— 1 



■3ri, 



^+ (a3„,) (^o.iri/ 



-3n, 



^rn— 1 



(157) 
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3.7 Commutation properties for null decomposition and 
Lie-derivatives 

We also need to relate estimates for the null components of the Weyl-field j^j^W 
to "slashed" derivatives of the null-curvature components. This relation is easily 
inferred from Proposition 7.3.1 of [8]: 

Lemma 3.8. Let u he any null- component of a Weyl-tensor W . Then 

u [CtW^ = 2tU 
where = denotes equality up to lower order terms of the form 
(a, a, /3, /3, p, (t) • (^decaying RRC, . 

Moreover, iftr^^^n = (as is assumed in the context of the ultimately Schwarz- 
schildean assumption), the p-term is absent and the lower order terms are 
quadratically decaying. 

3.8 Null components of the Bel-Robinson tensor 

We collect the components of the Bel-Robinson tensor in the null-frame: 

<33333 = 2|a|^ , (54444 = 2|a|^ , 

Q3334 = 4|^p , Q3444 = 4|/3p , 

(33344 = 4 -f fj2) , 
Qa44A = '^aAB^B , Qa333 = -'^^AB§_b ' 

Oa344 = ^pPa - 4.(T*Pa , (3a433 = -4p^^ - 4cr*^^ , 

QABiA = '2'5ab\P\^ + 2paAB - 2a*aAB , 
Qab33 = 2^As|^P + 2pa^g + 2a*aAB , 

QAB3i = 2 [SabP ■ P - Pa§_^ - PbP_^ + 25 AB (p' + d^) ■ 

3.9 The null-structure equations 

We collect the null structure equations (Proposition 7.4.1 of [8]). 

IP3H. = +trHH = -2'PIy_ - 2 VLH +{{Z + Z_- 2V) ®r) - a (158) 

p3 {trH) + i {trHf = 2divY - 2ntrH + 2Y- {Z + Z- 2V) - H- H (159) 

ci/rlY^YAiZ + Z-2V) (160) 

^iH + trHH = -2'PlZ + 2nH_ - ^trHH + {Y®Y) + {Z®Z) (161) 
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^4 (trH) + ^trH (trH) 
= 2divZ + 2ntrH - H ■ H + 2{Y ■¥ + Z- Z) + 2p (162) 

ci/rlZ^^HAH-YAY-a (163) 
{divH)^ - VbH^b - I'AtrH - VAtrH + (164) 
{divH)^ + VbHab = I^AtrH + F^^riJ + Pa (165) 

K = -\trHtrH + Ih -H - p (166) 

03 I/a = -2yAfl - (Vb + Zb) + 2Q {Va - Za) + HabY^ + 2r!r^ - 

(167) 

04"^^A = 21Ia^ - Hab {-Vb + Z^) + 2n {Va + Z^) - R^bYb - 2nYA - Pa 

(168) 

+ ^trHH = -2t>lZ + 2nH - ^trHH + {Y®Y) + {Z®Z) (169) 

03 {trH) + ^ {trH) trH = 
2divZ_ + 20. {trH) - H ■ H + 2 {Y ■ Y + Z ■ Z) + 2p (170) 

]piH + trHH = -2t>\Y - 2VlH + {{Z + ^+ 2V) ®Y) - a (171) 

04 {trH) + i {trH)^ = 2(i/t;y - 2ntrH + 2Y ■ {Z + Z + 2V) - H ■ H (172) 

^411^ - 03^A = 4f^i:A + Rab {Zb - Zb) - (173) 
IP3YA - 04^A - ^[IYa + Hab {Zb - Zg) + Pa (174) 

04fl + 03^ = F- y + y-(z-z)-z-z + Ann + p (175) 
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4 The energies 

We define the following energies for W: 

\\Wr:=\\a\\' + \\ar + \\Pr + \\f3r + \p\' + \a\\ 



(176) 



\\Wr := Waf + \\af + + + \p-p\' + \a\' 
Recall 2? = {^3, 04, y}. We define the higher order norms 
WV'^Wf := IIP'^Q'IP 



\V''af + \\V''P\\^ + \\V"P\\^ 
+ {p-p)\^ + \V''{a)\\ 



(177) 



(178) 



where I?'^ denotes any fc-tupel of operators from V. We avoid ambiguities re- 
garding the ordering of derivatives by understanding the definition as including 
all permutations of the k derivatives applied. 

4.1 The master energies 

We define the zeroth order energies 



^(0,P2,P3,P4,P5,0,0) 



,^,^[w]{j:r]= i www 

E^n{r<H} 



dvcLu 



rP''\\a\\^ +rP'\\P\\^ +rP^\\p,a\\^ +rP'-\\l3\\ 



(179) 



dudio 



WW' 



[W] {III) 



2 + rP«||a||2 



and weighted spacetime energy 



1° 



(0,P2,P3,P4,P5,P6,0) 



[W] LM(Ti,r2) = 



dt*drduj 



rP^WaW' 



M{ti,T2) 

+^P3 ||^||2 ^ ^P4||^^ ^j|2 ^ ^p, j|^||2 ^ ^p, ||^||2 



(180) 



(181) 



These energies should be understood as capturing the characteristic decay in 
r of the curvature components familiar from the stability of Minkowski space. 
Which tuples of p are admissible will be discussed in section 4.2. 

Recalling the notation (83) we define the following n*''-order weighted ener- 
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gies on the shces Er (n > 1): 



^ ' Js^n{r<i?} 
+ / dvduj 



2k2+2k3 



fel _ifc2_J,fc3 



, (182) 



I I ; ^ 1 "^^T? 



(0,pi,p2,P3,P4,P5,P6) 

fcl+fe2+fe3=n-l ■'■^^1 

^-'l /W, ,n \ ||9 . II -nlfel rkk^^^k^ 



,2fe2 +2fe3 



rP^||0^-^0:^y^^ (ya,/34) |P + r^'3||03^04=r^ (a,, ^/J, P4, ^4) ||^ 

+rf^||03'^0^^y'-^(/33,^,yp,ycT) f 

+r^-^||0^^04'^y'-^(p3,^3,y^,«4) f 



(183) 



and the 71 order weighted spacetime energies 



I 



■n,{non)deg 



[W] (m ( 



/ dt*drduj 

fei+fc2 + fc3=n-l"'^(''i'^2) 



rPMi03^04^y'-M«4)| 



«^(„o„)des (r) • r^'^^+^fc.^ 
+rP^||0,^^04'^y'^' (ya,/34) |p + rP3||03^^0,^^y'=-^ (ag, ^Z?, P4, ^4) f 

+r'^=||03"04^y'^ (^3,^3,^^,04) IP 



(184) 



with the weight w„ondeg {r) = 1 in the non-degenerate case and Wdeg = ^2^'' 



in the degenerate case. 

Remark 4.1. It is apparent from these energies that in the asymptotic region 
differentiating by 04 or ^ gains two powers of r. To avoid ambiguities in the 
ordering of the derivatives, the definition is understood as including all possible 
permutations of the derivatives applied. Such permutations only differ by lower 
order terms, as we shall see in detail later. 
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Remark 4.2. Note that because we are dealing with the characteristic energies, 
a4 is missing on I, while is missing on Nout {S"^ j^) • 

In addition, we define the summed energies 

n 

K [W] (S.) = E ,0) m (S.) (185) 

n 
i=0 

where P = {p\,P2tP\,p\,p\,p\,p\) denotes an (n + 1) x 7 matrix encoding the 
decay at each order. These matrices are defined in the foUowing section. For 
convenience, we also define the unweighted energies 

E" [M^] (e,) = ^ / llP^T^f EXCEPT Wlp^af on Nout (186) 
i=0 "^^^ 

n „ 

j„,(no„).e, ^^^^^ ^^(_).e,||I?^l^f (187) 

Finahy, for the horizon we define 

E"[W]in{Ti,T2)) ^ [ dt*duj\\V"Wf EXCEPT \\lp'^a\\^ (188) 

n „ 

l"[VF](H(ri,T2)) / deduj\\V'Wf EXCEPT li^galj^ (189) 

4.2 The decay matrices 

We need some definitions: 

Definition 4.3. We call a tuple (0,|)2,P3,P4,P5j 0, 0) boundary 0- admissible if 
all Pi > 2 and in addition 

2<p2<7-S , P3<min(6,p2) , P4<min(4,p3) , < 2 . 

To a boundary 0-admissible tuple p we associate the bulk 0-admissible tuple 

p= (0,P2-1,P2-1,P3-1,P4-1,P5-1,0) (190) 

in case thatp^ < 6, p4 < 4 andp2 < 2. In case thatps = 6 we replace p^ — 1 = 5 
by b — 5. In case that pi = 4 we replace p4 — 1 = 3 by 3 — S and in case p^ = 2 
we replace P5 — 1 = 1 by \ ~ 5. 

Definition 4.4. We call a tuple (pi,P2,P3,P4,P5,P6i 0) boundary admissible if 

2 < Pi < 10 - (5 , p2 < niin (9 - (5,pi) , < min (8 - (5,p2) 
P4 < min (6, Pa) , ps < min (4,p4) , Pe < 2 , 
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or if it is equal to the tuple (10 — (5, 9 — (5, 8, 6, 4, 2, 0). To a boundary admissible 
tuple we associate the bulk-admissible tuple 

P= {Pi- l,Pi - 1,P2 - l,P3 - 1,P4 - 1,P5 - 1,P6 - 1) 

m case i/iaf p3 < 8, p4 < 6, p,5 < 4 and pe < 2 hold. In case that p^ ^ S we 
replace p^ — 1 = 7 by 7 — 6. In case that p4 = 6 we replace p^ — 1 = 5 by 5 — S. 
In case that p^ = 4 we replace ^5 — 1 = 3 by 3 — S and in case Pq = 2 we replace 
PQ-l = l by 1 - S. 

Finally, an {n + l)x7 matrix P is called boundary admissible, if its first line 
is a boundary 0-admissible tuple, and the remaining n lines consist of identical 
boundary admissible tuples. A boundary admissible matrix P has an associated 
bulk- admissible matrix P: The n*'' line of P is defined to be the bulk admissible 
tuple associated with the n*'' line of P. 

We define the following (n + 1) x 7 boundary admissible decay-matrices: 

7-S 6 4 2 
10 -S 9-S 8 6 4 2 



Po = 



(191) 



Pi 



6-6 6 

9-5 9-6 8 



(5 4 2 
(5 6 4 2 



(192) 



5-6 5 

8-6 8-6 8 



(5 4 2 
(5 6 4 2 



(193) 



A-6 A-6 A- 6 2 
7-6 7-6 7-6 6 4 2 



(194) 



Wc will also need the auxiliary decay matrix 

/0444200 
6 6 6 6 4 2 



V 



(195) 



Note that the natural L^-based energy for the curvature components would 
correspond to a matrix P in (185) with all non-zero entries equal to 2 and for 
which the weight r'^^'^^'^^^ is dropped in (182). 

We also define an ordering relation on these matrices. We will say P < Q if 
all entries of P are smaller or equal than that of Q. 

Remark 4.5. Note that the energy induced by a boundary admissible decay ma- 
trix with second line (pi, ...,Pq,Q) controls in particular the lowest order energy 
induced by a boundary 0-admissible tuple with first line (0,p2 — 2, ...,P6 ~ 2,0) 

in view of the Poincare inequality /^a {f — T) r^dAs2 < cr^ Jg2 \^ f\'^'r^dAs2 . 
Cf. Lemma 10.5. 

Remark 4.6. Applying Sobolev inequalities one can obtain the familiar charac- 
teristic pointwise decay in r of the curvature components from the boundedness 

2 —2 

of the Ep^ [W] and Ip^ [W] energies. 
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4.3 Auxiliary norms 

Recall that we fixed 6 ~ . For W a Wcyl tensor and U a spacetime region 



hh\'+^m'+^\§\'+'^{\p\'+ wn I (196) 



E [W] (±r) = / r^dvdio {\\af + + mf + \pf + \af) 

r'drdio (||af + \\af + \\f5f + + \pf + \af) (197) 



i?[W](H(Ti,r2))= / dt*r'du; {\\a\\' + m\' + m\^ + \p\' + \a\^) 

J-H{ti,T2) 

Note that these energies are only expected to decay if W arose from commuta- 
tion with an approximate Killing field. 

In the context of the spin-reduction, it is convenient to introduce also certain 
weighted energies: 

Esp^ni m (S.) = / r^dvdco (||a|p + \\pf + + \p\' + \af)) 

T^drduj (||a||2 + + ||^j|2 + ||/j||2 + |^|2 + |^|2) 

s^n{r<i?,} 

which gives the components /3, p, a an additional weight. This additional 
weight has already been incorporated into the spacetime energy (196). For the 
spin 0-curvature components p and cr, more precisely the renormalized quantities 
^„ = p {0:^W) = r^p {C^W) - 2M(5^, V« = ^ (£^1^) = (£^1^) + 2M5'^ we 
define the energy 

E [V<t>n] itr) = [ drdij\V(j)nf + [ dvduj [1^4 0« I' + 

(198) 

and 

Ideg[V^n]iU)^ f dt* drduj Il-I-^JL ('(r0„)2 + (^0^) " 

Ju r \ 



(199) 



where T-*- = —pe^ + qe^ is perpendicular to T = pe^ -\- 564. Note that there 
is globally non-degenerate control of the T-'^-derivative only. The analogous 
definitions are made for ip^ (replace 0„ by V'n everywhere in the above). 



^''The reason that (3 docs not also have it is that /3 will not appear weighted on the charac- 
teristic hypcrsurface Nout- Cf. the spin reduction performed in section 11.1. 
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4.4 Constants and Conventions 

Recall that S = has already been fixed. Unless noted otherwise, B (b) 
will denote a large (small) constant depending only on the mass M. If the 
constant depends on other values, it will typically be indexed by them, e.g. Bh 
denotes a constant which also depends on h. The quantity e is the ultimate 
source of smallness and arises from the closeness assumptions in the definition 
of ultimately Schwarzschildean. We use the usual algebra of constants Be = e, 
B B = B etc. 

5 The main theorems 

Recall that the definition of an ultimately Schwarzschildean spacetime (Defi- 
nition 2.1) depends on the choice of an e > 0. Recall also the energies of the 
Ricci-coefficients defined in section 2.2.1. We have the following statement of 
uniform boundedncss: 

Theorem 5.1. [Boundedncss] There is an e > such that the following state- 
ment is true: If (TZ, g) is a spacetime which is ultimately Schwarzschildean to 
order k + 1 (for some k > 7) in the sense of Definition 2.1, then the estimate 



C" (n, T2) = sup e' [m] (SO + e' [m] (n (n, T2)) + 1 [5H] (M (n, T2)) 



and C'^ [91] (tq) = C*^ [91] (to,oo), which is bounded by the ultimately Schwarz- 
schildean property. 

We have formulated the statement in terms of the untilded A4,T, as it is 
easier to prove in this form. Of course, the estimate (200) then also holds for 
the tilded slices E and the D)" -energy. 

In other words, the energy involving fc-derivatives of curvature is uniformly 
bounded. Note that the assumptions regard only k-\-l derivatives of the metric, 
while the assertion provides bounds for k derivatives of curvature, i.e. k -\- 2 
derivatives of the metric. To prove Theorem 5.1 it will be sufficient to exploit 
multiplier estimates based on the ultimately Killing vectorfield T and a version 
of the rcdshift vectorfield. 

To formulate any type of decay statement one has to be aware of two im- 
portant issues, which can be heuristically understood as follows. 

1. As discussed in the introduction, the components p and a (and Lie-T- 
derivatives thereof) have to be renormalized and treated separately, if one 
wants to prove integrated decay estimates. These renormalized compo- 
nents, 4>, tpi may be thought of as satisfying an inhomogeneous Regge- 
Wheeler equation on a Schwarzschild background of the type 




(200) 



holds, where 



re(ri,T2) 



VRegge-Wheeler (</>) = r^ D {W ■ m) , 



(201) 
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the r^-weight arising from the renormalization. Clearly, to prove energy 
estimates for (jj one has to understand r-weighted curvature energies to es- 
timate the inhomogeneity. This means that one cannot understand local 
integrated decay without at the same time understanding the character- 
istic decay in r of the null-curvature components! 

2. It is well known that for the homogeneous Regge- Wheeler equation on 
Schwarzschild only the angular modes I > 2 oi (j> decay, while the I = 
and I = 1 modes describe perturbations to nearby Kerr solutions. Conse- 
quently, any decay statement will need to eliminate these modes. 

The first problem will be resolved by coupling the local integrated decay esti- 
mate - which will require an e of an r-weighted spacetime integral of curvature 
at infinity by the above heuristics - with the r-weighted multiplier estimates 
near infinity, as described in the introduction. 

The second problem will be resolved by the ultimately Schwarzschildean 
assumption, which will allow us to a-priori estimate the troublesome lowest 
order terms in an integrated decay estimate (these lowest order terms are the 
manifestation of the I = and I = 1 modes). We emphasize that in the case of 
the Regge- Wheeler equation on a fixed Schwarzschild background, these terms 
can only be controlled from the angular term by applying a Poincare inequality 
which excludes the I = and I = 1 modes. Here the situation is more favorable 
as we are assuming approach to Schwarzschild at a certain order, which takes 
care of the zeroth order terms. 

Theorem 5.2. [Integrated weighted energy decay] With the assumptions of The- 
orem 5.1 we have for any boundary- admissible decay matrix with P > Pp and 
for any n = 0, fc — 1 the estimate 

[W] (s.,) +T;''''' [W] [m (n,r,)) < B ■ [W] (S.,) 

+S-B"+M9^](ri,T2) (202) 

and, for any A > 0, 

Ip [W] [K) + 1^-°"^^^ [W] (m (n,r,)) < B, ■ [W] (e.,) 

+B-ni+^m{Ti,T2) , (203) 

where D| [% (ti,T2) was defined in (51). Here the constants B also depend on 
lower order energies of curvature and Ricci- coefficients, which are bounded by 
the ultimately Schwarzschildean assumption. 

Remark 5.3. Applying the estimate (202) for n ~ k — 1 from the data up to 
any slice in the future with decay matrix Pq provides a stronger ( than Theorem 
5.1) statement of boundedness, which moreover, as we will see, does not hinge 
on T being null on the horizon. It also yields boundedness of a degenerate (at 
the photon sphere) r-weighted spacetime integral (integrated decay). 

Remark 5.4. In fact, one can write Bx ■ Ip [W] [tr^^ + B ■ Wp"^ [W] {tr^^ 
for the first term on the right hand side of (203). 
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Note that the non-degenerate estimate (203) needs only the U""'"^ [5H]-energy 
on the right hand side, thus requiring only a small amount of the non-degenerate 
integrated decay estimate for n-derivatives of the Ricci-coefficients, cf. (51). 

Iterating Theorem 5.2 one can obtain 

Theorem 5.5. [Decay] With the assumptions as in Theorem 5.2, we also have, 
for n = 1, fc — 1 and any A > the estimates (t > tq) 

%[W] (tr)+T[W] {n{T,oo)) 

<^-K% [W] (So) +3^2 ^K-"' m ir) (204) 

i=i 

(.M(r,oo)) 

< ^ • Ep„ [W] (So) + S E 7iI°'a+' [^] (^) + ^ • ^" [^] (^) (205) 

i=l 

jn-l,„o„.e, ^^^^^^ 

< 1^ • lp„ [W] (So) + S E 77^A+' [^] (^) (206) 

1=0 

with 

( 1 ifn^k-l 
l = l{n):^ < 2 ifn^k- 2 (207) 
[3 ifn<k-2 

The constants B, Bx depend only on the mass M and lower order energies of 
the Ricci coefficients (i.e. D*^ {tq)) and curvature which are e-small by the 
ultimately Schwarzschildean assumption. 

In view of the powerful techniques of [8] , where estimates for the deformation 
tensor and the Ricci-coefHcients are obtained from appropriate estimates on the 
Weyl tensor, it is reasonable to state 

Conjecture 5.6. The spacetime of Theorem 5.5 is weakly ultimately Schwarz- 
schildean to order k + 2. 

Here by weakly ultimately Schwarzschildean we understand that the space- 
time is ultimately Schwarzschildean in the sense of Definition 2.1, except that 
the highest integrated decay bound on the Ricci components in (42) degen- 
erates at the photon sphere. This is of course a consequence of the fact that 
we can only prove a degenerate spacetime estimate for k derivatives of curvature. 

We close the section with a brief discussion of the relevance of these results 
to the fully non-linear problem, i.e. the problem of improving the assumptions 
made on the Ricci coefficients. For this it is essential to prove in addition that the 
spacetime under consideration is ultimately Schwarzschildean to order fc — 1 with 



41 



all constants and decay rates improved. As remarked after the definition (51), 
choosing t = Tf sufficiently large, the energy (t/) is e ■ A-small. Hence 

the estimate (203) implies that I^^ ^■'^""■'^'^a j^^j ^^^^ < e • A, as the first 
term on the right hand side of (203) is arbitrarily small up to time Tf by Cauchy 
stability. Moreover, the estimate (206) tehs us that '' \y\l\ (A^ (r, oo)j < 

eT~^. At least in principle - i.e. modulo the estimates for the Ricci-coefficients 
from curvature - this is sufficient to improve the assumptions at that order, in 
particular assumption (49) and the assumption l'" [*H] ^^M. (r, oo)^ < e on the 
Ricci-coefficients. 

There is one problem, however: We will not be able to improve the decay 
rate for the lowest order energy. In fact, due to the presence of the D" term 
in (205), the decay rate of the curvature can never be better than what we are 
assuming on the the lowest order energy of the Ricci-coefficients. The existence 
of this term originates from an error-term in the energy estimates, which is not 
cubic. As discussed in the introduction, we will show that this term is lower 
order in differentiability (allowing one to prove the Theorems above) but we will 
not be able to eliminate its non-cubic nature entirely. 

As a consequence, the results of the paper do not quite allow one to move 
immediately to the fully non-linear problem: There is still a missing ingredient 
regarding the nature of this error-term. This problem can in fact be studied in 
the context of a linearization of the equations. This is work in progress by the 
author and collaborators [14]. 



6 The T-energy 

In this section, we derive the energy of the n-commuted Weyl field HVfW as 
arising from the ultimately Killing field T. In Schwarzschild, this energy involves 
all components of curvature with good signs but with a degeneration of some 
components at the horizon. To stabilize the estimate under perturbations, we 
will have to invoke the redshift. Those estimates are carried out in section 6.3. 
In the following section, we simply allow for an e-small negative contribution 
near the horizon in our estimates (cf. Lemma 6.1). 

The uncommuted Weyl field (n = in the expression above) requires special 
attention. This is because the natural Bel-Robinson energy for W will involve 
the term J which does not decay. To remedy this we derive a renormalized 
energy directly from the (renormalized) Bianchi equations in section 6.4. 



6.1 The boundary term 

Let W denote any Weyl-field^^ with null-components a, /3, p, ct, /?, a. 

We choose A" = 3^ = Z = T in (5). A computation using the "Schwarz- 
schildean" normal h = ^^/k^e3 + \ /. ^4 (which is close to by the 



^^For prospective applications, the reader can think of W = CX,W for i > 1. 
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ultimately Schwarzschildean property) shows 



„3 



2„3 ^ 



Note that in the Schwarzschild case p vanishes on the horizon while q equals 1, 
leading to a degeneration of all components except a and /?. Since |p — fc~ | < e 

k+ . 

and !<? — < e by the ultimately Schwarzschildean assumption we conclude: 
Lemma 6.1. 

/ Q (T, T, T, ns) dfis^ > b [ \Wfr^drduj 

JSx Js^n{r>rY-^(rY-2M)} 

\W\^r^drduj (208) 



S^n{r<ry} 



and i/ie analogous identity for the slices T,r 



16 



The spacctime term Kj^^ \W] is discussed in section 6.2, while we postpone 
the analysis of the much more intricate error term A'J-^-^ [W] to section 7. 

6.2 The error Kj^^ [W] 
We have to understand 

m = 3Q"^^^ [W] ((^)7r„^) T^T, (209) 

for a Weyl-tcnsor W (cf. footnote 15). 

Proposition 6.2. T/ie error term Kj^^ [W] satisfies 

[ <e(Ti)-^ sup ^[W](e,) 

: / r^~^\\afdt*drdLj. (210) 



''"1 J M{Ti,T2)n{r>R} 

Moreover, if (ti,T2) is replaced by M (ti,T2) (^anc? Et- &j/ S^-J, t/ie Zast term 
can be dropped. 

Proof. In the region M. {ti , T2) f^ {r < R] and for (ti,T2) the estimate is 
easily obtained without the last term in view of the uniform {t*) ''-decay of all 
components of the deformation tensor, cf. (55). In the region r > R we have 
have to take into account that a docs not appear on the characteristic slices. 



^Note that the a component will not appear on Nout (^S^ ^j^. 
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which gives rise to the last term. In particular, using the formulae of section 
(3.8) and the null-components of the deformation tensor we can estimate the 
most difficult terms which arise from the contraction (209). 



p (W) a (W) • (^)i y^dt*drduj < 
f j du\ J dvduj^p{wf + e J \\aiW)\mt*y^ ^dt*drdLu 



<—= sup E[w,tr]+—=\ r^-^\\afdt*drdLj (211) 

and using that the component n of the deformation tensor decays like i a , 



j ||a(W)f |(^)n| ^dt^drdiu 



< -L / r^-^\\afdt*drduj. 



'^l J M(Ti,T2)n{r>R} 

□ 

6.3 A non-degenerate energy 

Note that in Schwarzschild the Ricci coefficients take the following values close 
to the horizon: 



trH^VLKi -— K -— and trH^n^O (212) 



4 — 7.2 47\// 

Define 



7(r) + — fc-j and y {r) ^ ^ (r) —k_ , (213) 

where ^ an interpolating function equal to 1 in r < ry and equal to zero for 
r > ry + and Cred = 2Sfc a small redshift constant (fc being the number 

of derivatives in the ultimately Schwarzschildean assumption). We will see the 
reason for that later.) For r < ry sufficiently close to the horizon we have in 
Schwarzschild 

-037 =~[dt-dr]l= ^ > 



|047l = |[fc+5t + fc-a,]7|< 



Cred r"^ 4,MCred 
1 



lOM 



These inequalities are stable under perturbation as we are working in a regular 
coordinate system. We summarize this as 

Lemma 6.3. Choose c^ed — 25^- can find ry > 2M such that in the region 
r < ry the following inequalities hold 

\lPil\ + m\ + \trH\<^^ , ^<-AI-trH<2, 
-037 > , - < -4Am < 2 , < i— ^ < — . 

4:MCred 2 Cred 10 
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Defining now the vectorfield 

r = + 7764 , (214) 

we observe using the formulae (66)- (70) in conjunction witlr the bounds of 
Lemma 6.3 tliat in r < ry the bounds 



8 A/ ' - ISMCred ' ~ ISAfc^erf ' 

|(^V^^^ab|<^ , r^-^'\\ + r'^^'\\<e (215) 
hold. Finally, we set 

N = T + Y, (216) 
which is future directed timelike everywhere, and state 
Proposition 6.4. The N -boundary-term satisfies 

Q m (TV, N, N, O > h l^lap + |«p + + |/3p + (p^ + j (217) 

everywhere on the black hole exterior. The quantity 

j^NNN ^ 3Q ((T+Y)^^ T + Y,Y + y) (218) 



^-TVATAr [IV] > 5 . Q [y^;] (jv, AT, AT, n^) /or r < ry 
\j^NNN [>v] I < s . Q [>v] (T, T, T, ns) /or ry < r < ry + 



2^ 



Proof. Since = (f + ^s+d + ^) ^4 = xe3+?;e4 is timelike, the bound (217) 
is immediate. For the bulk term, only the statement in r < ry is not obvious. 
For this we need to explicitly compute the contraction Qabcd [VV] ^^'^'it°-^N'^N'^. 

Qa,c, m (^K'^'^N'^N^ = Q3333(^^^''x' + ^3334 {2<^''^n^'x' + 2xi^^n'') 
+Q3434 + 4(^)7r34xy + (^)^44^2^ ^ Q^^^^ {2<^y)n'^y^ + 2^y\^^xy 

+ (34444^^^'^'^^/^ + (^'tt-^^ [QAB^-jx'^ + 2QABMXy + QAB44y 

Q3Acd m ^""K^^N^N^ + Q4Acd [W] (^V^^iV^TV^ . 

Note that the terms in the last line are e-small in view of (215). As x and y 
are uniformly bounded above and below and in view of the properties of '■^■'tt 
summarized in (215), we find (using the formulae for the components of the 
Bel-Robinson tensor collected in section 3.8) that the wrong-signed underlined 
term can be absorbed by the other, positive, terms. This yields the desired 
bound in r < ry. □ 

We remark that the above proposition should be thought of as the analogue 
of the redshift for the wave equation, cf. [1, 7]. 



45 



6.4 A (decaying) T-energy at the lowest order 



As mentioned at the beginning of this section, the T-energy at the lowest order 
involves the energy of p itself, which does not decay. We will now derive a 
renormalized T-energy at the lowest order, using the renormalized null Bianchi 
equations directly: 

Proposition 6.5. At the lowest order, we have the estimate 

I \\W\\^r^drduj < WWfr^drduj 

1 



+£ sup 



ri,T2 JS^n{r>ri'} 



\W\'r^drduj + B\\r'+'' {m - m 



ss) 



l+S 



\w\\ 



-B\\p' 



|^-5Hss|l 



where we recall that W contains the renormalized component p = p + . The 
same identity holds for the slices T,r and regions M (ti,T2). 

Proof. As in the previous section we will use the vectorficld T for computations 
and then argue by stability of the estimate in view of the closeness of T to 
T. Write the Bianchi equations as 



T{\\m 



kHrH + 2-^trH 



|/3|P = fc^ {Tp\{-p,a)+E^{P))p 



±- [diva + Ei (/?)] p (219) 



k^ 

k^trH + -^trH 



= -k^ 
k\ 



div^+Esip) 
divP + Ei [p) 



(220) 



kHrH+-^trH 



4 



c^rW + Es (cr) 
CT/rll3 + Ei (cr) 



(221) 



J^X 

2k^trH + -^trH 



11^11' = +fc^ [-diva + E3{p)]p 



k} 



(222) 



Ei (/?) ) /3 

"^x 

Note that the terms proportional to p^ and already decay in view of the 
assumptions on the Ricci coefficients. 

We now multiply the four equations by appropriate weights (depending only 
on r), add them up and integrate over the spacetime region A4 (ti,T2). The 
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weights can be read of from expression (208) and ensure that after integration 
by parts all spacetime terms which appear decay: Multiplying (219) by wi = 

, (220) and (221) by W2 = r"^ (||) and finally (222) by = 

(fc^)^ and integrating with respect to dt* dr duj yields the Proposition. 
As an example we look at the a spacetime term arising from equation (219) 
after integration by parts: 




The first term in the second line is zero in Schwarzschild and decays in the 
perturbed case, in view of the assumptions on the Ricci-coefficients. The second 
term is easily estimated by the terms appearing on the right hand side of the 
Proposition. □ 



7 The error 



In order for the identity (60) to be useful for us, we need to understand the terms 
j^xyz jyyj K^yz [W] on the right hand side. The two terms are of a differ- 
ent nature. If we commute with the ultimately Killing field T, K^'^'^ 



the error arising from the commutation and hence expected to be small for any 



X, y, Z with bounded coefficients. The term K^'^^ 



on the other hand, 



is only expected to be small if A", y, Z are ultimately Killing, cf. Proposition 
6.2. This entire section is devoted to controlling the error as 
arising from commutation with T . Remark: As this part of the argument does 
not depend on the assumption that tr^'^\ ~ (nor the gauge F = 0), we will 
derive all the formulae including this term. 



7.1 General null decomposition of J {T,W) 

Let W be any Weyl field which itself satisfies the inhomogeneous Bianchi equa- 
tion (W)^^^^ = .//37(5- We consider the commuted Bianchi equation (cf. (74)) 

D" (CtW) = Jp^s (T, W) + CtJp^s . (223) 

Assuming that J also arose from commutation with such a vectorfield, it be- 
comes clear that we have to understand the structure of the term J (T, W) and 
Lie-derivatives thereof. We define, following [8], 

(^)p, = i5"(^)5f„,, (224) 

'^^qa/i7 = - D^'-'^^^Pc - \ (^^^p75o/3 - '-^hpgo.^) , (225) 
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and decompose J (T, W) as 



where 



J (T, W^) = i ( Ji (T, VK) + J2 (T, M^)) 



(226) 
(227) 

(228) 



Note that q^^^ is a Wcyl current, i.e. it satisfies 

q/375 = -q/357 and = g'^'^q/s^s ^ -^q34S - ^q43S + S'^^qABS ■ (229) 

We recall that the error-term K2^^ [W] in the energy identity for the Bel- 
Robinson tensor arises from (T, W) + ,P (T, W) via formula (63). Hence we 
have to estimate the integral of 



D{T,W) iy,z,u)^Ki 



CtW 



(230) 



for appropriate vectorficlds y,Z,U. We proceed with a null decomposition of 
this term: 



Proposition 7.1. We have 

D{T, W^)333 = 4a (CtW ) • 6 (T, W) 



CtW] •S(r,M^) 



(231) 



DIT, M^)334 = 8p ( CtW ) • A (T, VF) - 8(t ( /ItM^ j • K (T, VK) 

-KT, VK) 



(232) 



DiT, W)iA-i = 8p ( CtW ) • A (T, VF) + 8ct (^^tV^J K (T, W^) 

+8/3 f £tW^) / (T, W^) 



D(T, M^)444 = 4a (£tM^) • 6 (T, M^) - 8/3 (^£tT^) • S (T, W^) 
where we have used the null decomposition of J {T, W) : 



A(J) = - J434 

K{J) = \e^^JAAB 

1 

^\J)a = -^JiiA 



A(J) = -J343 



K{J) 



AB 



J3AB 



S(J) = -J33^ 



Qab — 2 (^^4B + Jb4A — {5'"^ Jcad) 5ab) 

^AB ~ 2 ("^^33 + Jb3A — (S^^JcSd) Sab) 



(233) 
(234) 



(235) 
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Proof. See 8.1.7c in [8] 



□ 



Clearly, from the point of view of decay in the interior, the worst terms of 
J (T, W) are the terms in (T, W) which are proportional to the non-decaying 
component p.^^ We collect the null decomposition of both terms in the following 
Lemma. 

Lemma 7.2. We have the following formulae for the null-decomposition of 
(T, W) 



A(ji)^p 



--trH (^)5f34 _ 3^^^ (T)^33 _ 3^^^ / ^5 
2 2 — 8 V 



3 



-trH 



4 — 2 
K f Ji) =0 = 6 fji) 



(236) 

(237) 

(238) 
(239) 



where = denotes equality up to terms of the form ^'^'^tt-IJ) [W hut not p), i.e. de- 
caying quadratically as T is ultimately Killing. 

Proof. Direct computation. For instance, for the component A we find 

4A = ^^^DiWMM + ^^^£'4^^3434 + 5f'^^L'4l^A434 
+^3-^ £1^1^3434 + ^f^'^i:»3W^A434 + ^'^^£^^^5434 

= ^^^^A3434 + A^'^^Dip + 4^44^3p - er^^^^p - 6*yA^^^p 



(240) 



The other components are computed similarly, using the formulae in the ap- 
pendix. □ 

Lemma 7.3. We have the following formulae for the null-decomposition of 
J2 (T, W) 



A J' 



1 



D4 {trir) + D^iJma + '2pp 



2q(iA + i:'3^A4 + -^DAtriT 



-p [2ppA + Dii^_ 



Bab ( J^) = -^P [DAi'iB + Dsi'AA - Sab (£'^V'4a) - 2qaAB] 

3 

K ( J^) = --p- CAB [DAi'iB - pcre^s] 



(241) 
(242) 
(243) 



^^Note that (T, W) only has derivatives of p and is hence "a derivative better" than the 
terms we mentioned. 
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A 



1 



Bab {J^) = [DAip3B + DbiPsa - Sab {D^^i'sA) - Zpa^s] 

3 

K (J^) = -TP • <^AB [Da'^PsB + q(J£AB] 



(244) 

(245) 
(246) 
(247) 



where = denotes equality up to terms which do not involve the curvature com- 
ponent p and are hence of the form "{W but not p) ■ D'-^^^tt", i.e. decaying 
quadratically. 

Proof. Collecting the terms proportional to p for (T, W), we find^^ 

1 



A (J) -J434 p 



1 3^ 

2^4 + ■^Oba<^A4B 



1 3 

K (J) = -e J4AB = -^P<]A4B<^AB 



I {J) 



^iJ)A 
1 , 



1 



-J. 



2 

1 



34A 



44A 
1 

2 



:P^4A4 



® — 2 ("^-445 + Jb4A — [P'^^ Jc4d) SaB^ 

3 

= (qA4B + qB4A (qC4D<5'^^) (^As) • 



(248) 

(249) 
(250) 
(251) 

(252) 



with the analogous formulae for the bared quantities obtained form interchang- 
ing 3's and 4's in these formulae. We define the antisymmetric tensor 



(253) 



and compute, using the definition of the deformation tensor and commuting 
covariant derivatives 



(\A4B = DiTTBA - DBTr4A + gp4'5AB 

^ i {D4DbTa + D4DATB - DbD4Ta - DbDaT^} - ^D^trTrSAB 
1 



:P4Sab 



2 yR4BA3T^ + R4BA4T'^ + R4AB{a)T'-°''' + DaD-iTb — RbA43T^ — DaDbT^ 

^DitrnSAB + \p45ab 

-^SabDa {trn) + Da^Jab + ^ i-2qaAB + '^PP^ab - 2paeAB) + ^PaSab 



-"^^Note the typo regarding the ©-equation in [8], where the factor of i is missing. 
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qiAi = Datth — DiT^Ai 

= i [2DADm - DiDATi - DWaTa] 

= i [D4DAT4 - DiDiTA + 2Ram3T^] = 2pI3a + Da^jaa 



^ ^ 1 12 

q34A = £'47rA3 - £'A7r43 - 2pAg34 = -2(7/3a - D'iijjAi - -DAtrn + -pA ■ 

as well as 

qA3S = -'^SabDs {trn) + i:'3V'3S + ^ (-2paAB + 2qpSAB + 2qGeAB) + \p3^AB 



q3A3 = -2g/3, + £'3V'A3 



q43A = 2p^^ - i:>4V'A3 - -^DAtrn + -pA . 



from which the formulae of the lemma are easily derived. 



□ 



7.2 The structure of J (T, W) 

We combine Lemmas 7.2 and 7.3 expressing the derivatives in terms of the 
slashed derivatives and the null component Va ■= i'AA and Qa 4>A3- 

Lemma 7.4. Let {tt, dec. RRC} denote any term which is of the form "con- 
stant times a component of the deformation tensor" or "constant times a de- 
caying Ricci rotation coefficient" . We have the following formulae for the null- 
decomposition of (T, W) 



-04(tr^)-y Va + 2p[p 



I{J)a^IP 



{tt, dec. RRC} 



'-^AtrTT + IPsVa + 2qPA + {tt, dec. RRC} 



(254) 
(255) 



and 
A (J) 



S ( J)^ = --p [1P4VA + 2ppA + dec. RRC}] 
Sab {J) = Ip \^t>lVA - 2qaAB + dec. RRC} 



K{J) 



:P ■ ^AB [-^ a'Pb - pc^AB + {tt, dec. RRC}] 



/a {trir) - 'JZ^^Q^ + 2fy ( p + ^ ) + {tt, dec. RRC} 



L{J)a^IP 



-llAtrTT ^IP^Qa- 2p§_^ + {tt, dec. RRC} 



(256) 
(257) 
(258) 

(259) 
(260) 
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S(J)^ = --p [iP:,Qa - 2q§_^ + dec. RRC\\ (261) 
Qab {J) = \p [2^2 Qa - 2paAB + {tt, dec. RRC}] (262) 
K{J) = -^p- CAB [-^aQb + qaeAB + dec. RRC}] (263) 



Proof. Note that 



= f^ipAi + \trH (203P -4:Qp + 7734) 



(264) 



and 



Dz^Ai = ^z'ipAi - ZAipai - 2ZBi>AB - 2[hpA4 

DilpAi = lp4i>A4 - ^aV'34 - 2YbV'AB + 2nipA4 



(265) 
(266) 



DAtpiB + DsipiA - Sab {D^ip4A) = -2t>l-tp4A - HabiPs^ - 2HcbiPac 

= 2t>llpA4- HABi^SA (267) 

and similarly for exchanging 3 and 4 indices. Hence for A we find 



A (J') 



^ " 34 - -trHtrn - -trHirn - -trHS^^TTAB + -trHtrn 
-^ 4 2—2 -^-^ ' 4 



-i04 (trn) - i/^VA + 2pp - ^trH (2^3^ - 4Qp + 7134) 



which upon adding produces 

A ( J) = [ - ^04 (trn) - I/^Va + 2pp - UrH {21/ ^p - A^lp) 

-trxT^M - ^trH_n44 - ^trHS'^^iTAB 

Since up to decaying Ricci coefficients 2pp~ \trH (21/ ^p — 4f2p) = 2p (p + ^3^) 
the result follows. For later purposes we also collect the explicit form of /. It is 



L{J)a^iP 



iy^irTT + IP{Pa + 2ql3A - Za^im - 20Pa - trHji^'^ - 2trH^*'^ 



□ 



52 



7.3 K^y^ 



Using the previous Lemma we can now look at the spacetime integrands D{T, W). 

Proposition 7.5. The terms D{T, W) arising from (T, W) + ,P (T, W) sat- 
isfy 



D (T, = 604 (poP (CtW) ■ r) + Qdh ip^a [CtW] ■ V 



-'SCi 



qp\a{W)\^-2pp\PiW)\ 



(268) 



D (T, W)443 = 603 (^poV ■ 13 (CtW)) - 304 (po • tr^^'^^^P (CtW)) 
-6|>i (-Va Pa P (£tM^) ,rA Pq<t (^tW^)) 

+6Ct{pp(p+'^^ + ppoa^ + qpa\l3\'') 



(269) 



D (T, W^)333 = -603 (po^ (^CtW^ ■ S) + 6div (poa [CtW^ ■ V 

-^CtUpMW)]'' ^2qp\p{W)\'' 



(270) 



D (T, W^)334 = -604 [poQ-§_ [CtW)) - 303 (po • tr^^^^^P [CtW 
+&Tt)i (-Qa Po P (CtW^ ,-Qa Po (CtW 

+6CT{qp(^P+'^^ + qpo^^ + PPoWf) (271) 

where = denotes equality up to terms of the following form 

(^CrW'j -{poTT , padec. RRC} (272) 

and lower order terms of the form 

(w-Wyin, dec. RRC} (273) 

Proof. First note that aU terms arising from (T, W) are of the form (272), 
(273) or even lower order. Denote by po = p (W) the non-decaying p component 
of the uncommuted Weyl tensor we compute, modulo cubic terms in which all 

three components decay and terms of the form ^£^1^^ ■ {7r,dec. RRC} • po: 

D (T, VF)444 = 4« (^CtW) ■ ^po {PIVa - qaAs) 
-8/3 (ZtW) ■ (-^pQ {2pPa + IP^Va) 
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Using Lemma 3.8 we can write 

D(T, 1^)444 = 6gpo 

6(00 



Va- — {1P4P0) 13 [CtWJ ■ Va 



diva [CtW^ - 1^4^(3 (^CtW 

+604 (^pol3 (£tVk) • Va) + Gd^v (po" (^CtW) -V^-Ga (^CtW) fpoV 
and hence 



D [T, = 6po 



{ 2trH + 2n + -trH ] P [CtW] - a (CtW ) fpo - Jaaa 



Va 



PoP (^CtW) ■ v) + Qdlv {p^a [CtW) ■ v) - SCt [qpoa^ ~ Spp/?^] 

Note that the terms in the last hne are pure derivatives. Here we have used in 
particular the Bianchi equation (86) and the fact that the adjoint of ^2 div. 
Similarly, again modulo cubic terms in which all three components decay 

and terms of the form [CtW] ■ {7r,dec. RRC} • pq: 



D (T, W) 



443 



CtW] ■ -po 



2M\ 
—J 



+8a [CtW^ (^--po (-2pa - c4rl Va) 
+8/3 ( CtW) ■ ^po UqPA + 03^A + ^'fAtr^ 



12ppo 



, 2M\ f 2M\ ( ^ (2M\ 1 



\ 7" 



+12ppo 
^^PqVa 



Uqpo 



1 [-p[ CtW), a [CtW] - 03/3 £tVI^ 



+'ipotr'^'^\ 



04P [CtW] - divl3 [CtW 



+603 (po^ • /3 (CtW)) - ilPi (po • t^^^^P (CtW)) 

-61^1 (^~Va Po P [CtW) , +Va Po o (CtW)) . (274) 

Using the Bianchi equations (88) and (91) yields the statement in the lemma. 
The computation for D (T, ^1^)333 and D {T, W^)334 proceeds analogously. □ 

Note that since these terms are going to be integrated, we have essentially 
gained a derivative: The worst error-terms involving po are not of the form 

Po ■ u (jItW) ■ Dtt but only po • u (^CtW) • tt. It is not hard to see that the 

structure we revealed at the first commutation with T, survives to higher orders. 
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7.4 K^y^ 



Consider the n + 1 times T-commuted Bianchi equation: 

n-1 

D" (Cp^w) = C^Jp^s (T, W)+y^ CtrJpyS (t, 



i=0 



W] = J' 



jhard | jeasy 



We use the index "easy" for the second (summed) term, in view of the fact 
that a T-derivative has aheady fahen on W in this expression. Hence this term 
already decays quadratically. On the other hand, the most difficult term in 
jhard ^^iq one where all derivatives fall on the deformation tensor, in view of 
the non-decaying component p of W. We set 



easy j^XyZ 



hardj^XyZ 



(275) 



with the terms arising from the different terms on the right hand side, J*^"*" 
and j'^°-^<^ via formula (63). In view of the quadratically decaying structure of 
jeasy easily see 

Proposition 7.6. Let (7?., g) he ultimately Schwarzschildean to order k + 1 
(k > 7). For 1 < n < k — 1 and any A > 0, 



^ ^ easy j^e^e^Ca 



< A 



sup E"*^' [W] (s,) 



a,b,e€{3A} 

1 7) / ~ 

+ Bxe{Ti)'^ -supE [m] [^r 
dt*drdujr^ (||I?"+iaii2^ 



(276) 



/'''I JMn{r>FI} 

Moreover, as in Proposition 6.2, if Ai (ti,T2) is replaced by M. (ti,T2) (and S,- 

CtW 



by Tir), the last term can be dropped. For n = 0, easy p^eae^e^ 
Proof. Observe that schematically 



= 0. 



n— 1 

J,asy = [{d^T^) (4"'^) + (^T^) (dO^T^w)] , (277) 

where the = ignores all lower order terms which arise from the commutation 
of D with the Lie derivative^^ and all combinatorial (integer) factors for these 
terms. We then estimate, first in r < R: 



E 



easy j^Ca eb 



< I dt*\\Cp^W\\L2 



Mirr,r,)n{r<R} „_f„c6{3,4} 
s=Q 

n-1 



'these terms will not only be of lower order but also introduce additional decay 
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to which we apply Cauchy's inequahty and Sobolev embedding: 



< A • sup e"^^ [W] ( 



-Sa - Tif ■ sup E'+-^ {tr) ■ sup E"+l-' [W] (tr 



i=0 

n-1 



+Bx -nf- sup E"-'+3 [VF] (S,) • sup E'+i (e,) . 

One easily checks that for = fc — 1 = 7 (the worst case) the ultimately 
Schwarzschildean assumption ensures that one can absorb the (t2 — ti) -term 
using the decay of one of the swp-terms. 

We can do the same estimate in r > R, the only thing we have to be careful 
about is that not all null-components appear on characteristic slices. This is 
easily accounted for by adding the spacetime term appearing on the right hand 
side of (276). Note in this context that potentially divergent terms like 

J dt*drdLu r2 a (l^+^W^ ■ (iP^iC^^^ (^i)) • Cra (278) 

(i.e. where all terms of the cubic expression only decay like -) cannot appear 
in K2 [W] because of the null-structure of the error-terms (cf. the signature 
considerations in [8]). At least one component in the cubic error-term has 
improved decay in r. □ 

For hardj^xyz ^^^^ 

Proposition 7.7. Let {TZ,g) be ultimately Schwarzschildean to order k + 1 
(k > 7) and X = p^-es + (?A'e4, y — Vye-i + qye^, 2 = vz&z + ^lz&A: be vectorfields 
for bounded spacetime functions pi and qi satisfying \T {pi) \ + |T (g^) | < e • t~3 
in the interior r < t* . For < n < k — 1 and any positive A2 we have 



f hardj^XyZ 
JM 



< (As + e (ri)~^) sup E'"*"^ [W] ( 

^j.+ l,nonde, ^^^^ ^ ^ -„ ^^^^^^^ 



^A±\\c^wr+\\wA 

/ dt*drdwr5 (||I?"+ia||2) +A2Sup|k| / WC'^+^Wf . 

J Mr\{r>R} H,i Jn 



(279) 



As in Proposition 7.6, if M (ti,T2) is replaced by M (ti,T2) (and E,- by Yjt), 
the penultimate term can be dropped. 

Proof. Observe that the null-decomposition oi J = J (T, W) almost commutes 
with Ct'. 

k{CtJ) =£y(A(J))-|-(decayingRRC)-(A(J),....,£(J)) , (280) 



56 



where the error is of lower oder and introduces additional decay. Similarly for 
the other null-components. Now the £^-derivative of each null-component of J 
is of one of the following forms 



ttLtDp + Dtt ■ CtP 



(1) [t:D {W but not p) + Dtt ■ {W but not p)] or 

i.e. a product of two (or more, arising from commutation) decaying com- 
ponents 

(2) £5;7r Dp 

(3) C^Dt: p (these terms were essentially computed in Lemma 7.4) 

In the expression for A'2 arising from these terms, terms of type (1) can be esti- 
mated as in ^°-''yK2-^'^ We only mention the term involving the highest derivative 
of the weakly decaying component ^i. In view of signature considerations this 
term is of the form^^ 



dt*drduj r a I Cj, 



r>R 

\r'^[W but not a,p)\\L 

-n+l^nondeg 



n+1 



W 



^a^^ (^(^ij j • (VF but not a,p) 
dfdrdijri (||I?"+V||2 + IIP' 



e 



I 



Mn{r>Ft.} 

m{M (ri,r2) 



dt^drdur^WV+^al 



Mn{r>R} 



For the terms of type (2), we first note that using the Bianchi equations 
for p, they can be transformed into terms of type (1) and terms of the form 
C^TT {trH, trH_) p. Realizing that the decay in r at infinity is not an issue for 
the latter term in view of p decaying like we integrate it by parts in T to 
obtain, modulo lower oder boundary terms which are already present on the 
right hand side of (279): 



C^TT po] 



Po 



(281) 



for n > 1, while the right hand side equal to /_^;((^^ y-'T''^ Poj W, in case that 
n — 0. An application of Cauchy's inequality yields the terms found on the 
right hand side of Proposition 7.7. 

Finally, for the terms of type (3), we can redo the computation of Proposition 
7.5 and transform this term into terms of type (1) and (2): 

Lemma 7.8. We have the following generalization of Proposition 7.5 



-3Ci 



qp\a |2 - 2pp\p (^C^W) 



■^'^Note again that, in view of signature considerations, terms involving twice the curvature 
component a in conjunction with the (weakly r-decaying) null-component cannot appear. 

^^Note that this term is easily estimated if Ai is replaced by A1, in view of the uniform 
decay of the deformation tensor (55). The problem here is that a does not appear on the 

characteristic hypcrsurfaccs Nout [S^ j^j . 
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D (T, VK)443 = 603 (^poC^V ■ /? (4+1 W^)) - 304 (po • Cl}.tA^K p (^CtW)) 

-0^ (^-c^Va po p ,c^Va po <y [cf^'w)) 

-ppoa (^C'^Wy + qpo\\C'^f3\\^) 



-6Ct{pp 



2M\ 



and similarly for the other two components. Here = denotes equality up to lower 
order terms which are of the form arising from (1) or (2). 

Proof. Wc show the statement only for the first component, as the others are 
treated completely analogously. 

D (t, C^W^ = 4a (^C^+^w'^ ■ 6 (T, C^W) - 8/3 • S (T, £^1^) 

= 4a V) • £?;e (T, W) - 8(3 V) • £JS (T, W) 

Now we can use Lemma 7.4, where we computed the part of the null-components 
and S, which is proportional to p. Since commuting with the Lie-T-derivative 
only introduces lower order terms we obtain 

D (t, d^W^^^ EE 4a {t^+'^W^ ■ ^po (tlC^VA - qa (fjM/^ 

-8/3 (4+1 ly) • (^-^po (2p/3 (4W^) + 044^a)) , (282) 

from which point on we can follow the computations of Proposition 7.5 to pro- 
duce the result. □ 

We now integrate the main terms collected in the previous Lemma. Since 
they arise from 

D (T,C'r}.W^ {pxe3 + qxei,pye3 + qye4,pze3 + qzei) , (283) 

we observe that for vcctorfields satisfying pi = {i = X, y, Z) on the horizon, 
the boundary term on the horizon will vanish. In general, the strength of the 
horizon term will be of the size of the pi's on the horizon. Hence, for the terms 
appearing in the Lemma 



/ D (T,Cr^w) (pies 9164,^263 + (7264,^363 + 9364 

JM{ti,T2) ^ ' 

<i?-sup|p,| / U^wZv^+^wf ^^WA 



+A2 sup E"+' [W'] (E,) + Sa, sup E"[fH](I]0 

r6(ri,T2) T6(Ti,r2) 

-|-spacetime terms of the form (281) (284) 

□ 
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7.5 Summary of the error estimates 



Combining Propositions 7.6 and 7.7 wc can summarize our estimate of the er- 
rorterm 7^2- We note in passing that the curvature terms in Proposition 7.7 can 
be controlled by the energies of the Ricci-coefficicnts one order higher, simply 
by the definition of curvature in terms of derivatives of the Ricci-coefficients 
(cf. the null structure equations in section 3.9). Hence: 

Proposition 7.9. Let (72,, g) be ultimately Schwarzschildean to order k + 1 
(k > 7) and X = Pxe^ + qxe^, y = Pye^ + qyCi, Z = pzea + lz^i be vectorfields 
for bounded spacetime functions pi and qi satisfying \T (pi) \ + \T (qi) \ < e • * . 
For any positive A and < n < k ~ 1 we have 

I K^^^\C1}+^W] < (X + Bxeinr^) ■ sup E"^\w] (tr) 

+B ■ [% (n, T2) + Ba • m (n, T2) 

+e(ri)-^ • [M/] (^(Ti,r2)) +Asupb,| / \\Cp-'wf 

^ ' H.i JU 

(285) 



For X = y = Z = T the last term vanishes and no highest order curvature-term 
appears on the horizon. 

Finally, if (ti, T2) is replaced by M (ti, T2), St by and the I}-energies 
by the C- energies, the estimate holds without the penultimate term on the right 
hand side. 



8 The redshift (from the null-Bianchi equations) 

8.1 The main result 

The main result of this section is 

Proposition 8.1 (Estimates close to the horizon). We have for < n < k — 1: 

Kil. m (StJ + 1"+' [W] {n {t,,t,)) + Xtlf [W] {M (ri, r2)) < 

2 . 1::;^ [W] (StJ + B ■ i;.t';';<i^ m {M (n, r2)) + Err^l^ (n, r2) 

(286) 

with the error 

Errlt' m (n,r2) = e Xttl% I^l (-^ (^i'^^)) 
+B ■ lr<4+M m {M {n,T2))+B- E,<,^+M (St,,I]t,,H) (287) 

where e arises from a pointwise estimate for W and the constants B depend only 

^ r 7 . /^^^\ , —max{2,n\,deq , , , , ^^ 

M and, m (286), the I^<^^+m [% [M {ti,T2))- energy. 

Proof. This will follow from the sequence of estimates proven in the remainder 
of this section. □ 

Note that we are estimating n + 1 derivatives of curvature, requiring only 
an e of n + 1-derivatives of the Ricci rotation coefficients. 
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8.2 Estimates for the uncommuted equations 

Recall that we defined 

1 



Cred 



(288) 



in (213). Multiply (87) by -fa and integrate the resulting equation over the 
spacetime region y = M. (ti,T2) C\{r < ry}- Integrate the angular /3 term by 
parts (producing diva) yields upon inserting (88) the identity 



7 - 7 ( 4r2 - -trH 



03 (711^11') + [-037 + 7(4fl + 4tr^)] M?)] = ^dt* drdtu [a,^ 

with the error 

e-y [a, P]=j [E^ (a) - 417a] a + 27 [£^3 {0) + 2^ §_ , 

for which we ignore the entirely harmless cubic term 2afi'^ ^7-^=^ arising from 

the integration by parts (this term can always be estimated as e (||a|P + 
using the ultimately Schwarzschildcan assumption). Inserting the relations 

^3 (7^') = (7^') = Da {{e^T 10^) - 10^ [-2^1 + trH] 

04 (7"^) = Di (7^2) ^ Da {{cif 7a2) _ [-2Vl + trH] (289) 

back into the above identity we can write 

j^^dt*drdLo(^-Da{{e,ri\\af)+ -1047^7(30) \\af + 
Daiie^iTlWPf) + [-IPzl + im+'^trH)] \\§\\^)] = f ^dt*drdue^ [a,§} . 



(290) 



By Lemma 6.3, in the region under consideration both 

- ^047 - Sjfl >b and - ^37 + 7 (6fl + 2trH_) > b 



hold and we obtain the estimate^^ 



< B 



Sx2n{r<ry} 
ST^n{r<ry} 



B 



A^n{|r-ri^|<- 



n{Ti,T2) 



+B 



y 



[a,P} (291) 



^^We have not written the measure exphcitly here as no confusion can arise from weights 
in r in this region. 
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The same procedure applied to the Bianchi equations (90), (93) and (94) 
produces the identity 



^ y/gdt*drdLO I hi 



+ 2^3 (7/5' +7'^') + 



i^47 + 7 (2tri7 - 2f]) 



^037+^7*^^ 



yfgdt* drduj (p, a) 



with the error 

\l, (p, a)] = (p) P + 7^ (a) a + 7 [ii;4 (^) - 2f]^ ^ 



(292) 



(293) 



(again we ignore a harmless cubic error-term of the form /3 {p,a)'^ ^^7-^=^, 
which is readily estimated). Upon integration over y this leads to the estimate 



J'S^^n{r<rY} 



/3| 



(ll^l 



< B 



^rin{r<rY} 



im 



H{ti, 
2 



-B 



m' + \\P^^\\')+B / e-,[/3,(p,a)] . (294) 

Mn{\r-rY\< ''^ 2 } 



where we have taken into account that both inequalities 



1 



-7:041 + 1 \-^ + TitrH] >b and - $37 + 1 [ +n + -trH] > b 



hold in r < ry 

The other two sets of Bianchi equations are a little more subtle. We can 
not achieve a good sign for both of the spacetime terms at the same time. For 
instance, using the renormalized Bianchi equation (95) and (92) we derive 



J ^dt*drdLoi^0^ (7p'+7<7' 



1 



03 {iP^ 



1 



hi - 7 (2il - trH) 



13']= ^dedrdwe^[{p,a),p] 



y 



with the error 

[(p, a) , /3] = ih (p) p + (fx) <T + 7 [E3 W) - 2m /? (295) 

After integration, the /3 spacetime term will admit a good sign, while the {p, a)- 
term will have the wrong sign as 



i 047 + 7 (tri/ + r!) 



(296) 
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is negative close to the horizon. However, since (/5, a) is aheady controlled by 
the previous estimate, we can add a multiple of the estimate (294) to also control 
the spacetime term of /3. Similarly, for the final set of equations: 



1 



1 



-T037 + 7 -:trHa - 2na 



P^\= j V9dt*drduje^[/3,a] (297) 



with error 



[/3, a] = -7 [E^ (a) - Ana] + 7 [S4 (/3) + 2^/3] /? 



(298) 



Hence while the a-spacctime term is positive, the /3-spacctimc term has the 
wrong sign (which will only get strenghtened after the contribution from the 
integration by parts). However, the /3-term has been controlled by the previous 
step and we can add a multiple of the latter to control this term. Adding the 
estimates up we can summarize our findings in the following 



Proposition 8.2. We have the estimate 



T,r2n{r<rY} 

B 



\W\\ 



\Wf < B 



S,in{r<rv} 



Mr\{\r--rY\<- 



-} 



\w\\ 



-B 



y 



67 (p, o-)] 



B 



B 



B 



For the error-terms we have 
Lemma 8.3. 



67 [A <t)] 



[/3,< 



y 



< \\{9\-y\ssf\\L^ 



\wr+B\\p' 



\m-mss\ 



(299) 



(300) 



Proof. This follows from inspecting the error-terms individually. Most terms 
are of the form "decaying Ricci rotation coefficient" • "decaying curvature com- 
ponent" and are estimated by the first term on the right hand side of the lemma. 
The only other term arises when the curvature component is p. This is accounted 
for by the last term. □ 



8.3 The higher derivative redshift estimate 

The estimate of Proposition 8.2 can be derived for all higher derivatives as well, 
this being essentially the statement of Proposition 8.1. We will now prove the 
latter. 
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Estimating ag^^ and 

From the equations for 13^^^ and a43„j , derived from (116) and (122) we obtain, 
modulo terms which will vanish after integration and a cubic term involving the 
angular derivative of the measure 



^4117^3,1,1 



^3ll7^3,. 



2||/3, 



^7 + ^^ {^,^)tr{H)r 



-047 + ^+ (a3nj itrH] = [fi3«, ,^3„J (301) 



where 



74'"' (^i)^i3. 



27^;^"' 



27C33 /3„ /3 



—ini 



For the latter expression we prove 
Lemma 8.4. For any A > 0, 



1 ' :i3n. 



< 



[6 + 2 {number of 3 's in ni)] 7 ri|la 



^2(nwm6ero/4'5mnO]7f^ll^„JI' + A I (ll«3,JI' + ll^3„J 



■3ri; 



B\\ (5K- 



5K, 



■l~\-l,deg 



SS) \\L^ ■ ^r<ry 

■l-\-l,deg 
|r-rv|< 

— /+l,(iet? 



Vt^] (X (Ti,r2)) 
-2M [W^] (A^ (Ti,r2)) 



-B\\W\ 



-l.deg 



[in,M^] (X (ri,r2))+E!:J,t [$K,W-] (S 



[fH](Al (Ti,r2)) 



Remark 8.5. The first term on the right hand side is the redshift term, which 
has a good sign when brought to the left hand side?^ The contribution for (3^^ 

has the wrong sign. However, since we have the large term —037 already avail- 
able on the left hand side of ( 301 ), we will absorb the negative contribution. The 
terms in lines two and three will eventually be absorbed by the left hand side 
(once we added the estimates for all the components). The remaining terms 
involving Ricci- coefficients are contained in the expression for (287). The re- 
maining lower order curvature terms will eventually be controlled by reiterating 
the redshift estimate, while the term in the fourth line is present in (286). 



Proof RccaU that E^f'' (a) = eI'''' (a)-C34 [«„,] 
the commutation term 



Li£43 



£34 



. Starting with 

ni o 



-C34 [a„J a 



S,ni 



[a 



,34n, 



£3n, 



we see that we need to push through the 4-derivativc. Every time it hits a 
3, Lemma 3.5 applies and introduces a lower oder term. Since terms which 
contain I derivatives are lower order (and accounted for by the term in the 



^^As is familiar from the wave equation, its strength improves, the more 3-derivatives are 
taken. 
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penultimate line of Lemma 8.4), we only have to focus on the highest order 
term. Furthermore, as both J7 and Z_ decay, we only need to follow the term 
which is proportional to in the expression -F34Ai...Afc of Lemma 3.5. It is easy 
to see that this term will be 

-C34 [a„J "an, = - ("34n, " a3«,4) ag™, ~ (number of 3's in n,) 2^a^^^a^^^ 

The C33 -term is dealt with analogously and gives a negative contribution. 

Inspecting the term iJg"' (/?) and how it arises inductively from derivatives 
of £"3 (/3) we observe that it consists of the following type of terms (note that 
£1 decays) 

• products of curvature components with Ricci coefficients with both com- 
ponents decaying 

• derivatives falling on the expression —"SYp in £3 which are problem- 
atic in view of the fact that p does not decay. 

The terms of the first type are unproblematic if sufficiently many derivatives are 
taken and are easily seen to be estimated by the terms of the right hand side in 
the lemma. For the terms of the second type we make the following observation: 
The worst term is clearly the one where all derivatives fall on y . This is because 
if a derivative falls on p it can be replaced by p and terms which decay using 
the Bianchi equation, hence gaining one derivative. It follows that except for 
the term where all derivatives fall on y, we can estimate p pointwise and fc- 
derivatives of Y_ using the spacetime energy for the Ricci coefficients appearing 
on the right hand side of the lemma. We still need to estimate the term 

^7P[0"'03i:]^3„, (302) 

For this we will be able to exploit a cancellation with the other a-terms. We 
first note from equation (131) and Lemma 3.5 that 

- [043 - ^34]„, = +2^^a3n, - 2fla4„^ + terms prop, to (303) 

Furthermore, inspecting how £4"' (a) arises inductively from derivatives of 
£"4 (a) , we observe that it consists of terms of the following type 

1. products of curvature components with Ricci coefficients with both com- 
ponents decaying 

2. derivatives falling on the expression 4fia in £4 (a), which in view of the 
fact that r2 itself does not decay are problematic. 

3. derivatives falling on the expression —3Hp in £4 (a), which in view of the 
fact that p does not decay are problematic. 

Again, the terms of the first type are easily estimated. For the terms of the 
second type we observe 

Derivatives all falling on the same factor. If all derivatives fall on a we find a 
contribution of 4f2a3,i, , which adds to the 2f2a3„j in (303) to produce the factor 
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of 6 in the Lemma. If alH + 1 derivatives fall on we can (since the I + 1- 
derivatives have at least one 3 derivative in it) repeatedly use the commutation 
formula 

03041) ^ JpilPsfl - 21703^ + 2fl04fi + {Zb - Zb) fs^ (304) 

and the structure equation (175) to create terms of the first type and a term 
where I derivatives fall on p (recall that J7 decays) . For the latter we can again 
use the Bianchi equation for or to gain a derivative or transform it into a 
cubic term. 

Derivatives splitting up. If the derivatives split, we can estimate one of the 
terms pointwise and the other in L^. Writing fl = fl + ^■ — ^■ one sees that 



l-s+1 

a 



< A ||a3,jP + f 



t'^' [w] (3;) + supEt<,^ miJ^r) 



where B also depends on the D' [fH]-cnergy which is bounded by the ultimately 
Schwarzschildean assumption. 

For the terms of the third type we argue similarly as for the /3-term, 
which leaves us to control the term where all derivatives fall on _ff . Combining 
it with the (3 term we need to estimate 

1 = j^ip[^ [2030"' y] ^3,,^^ + [0"' 03l] a3„,^ (305) 

Note that we can choose the ordering of derivatives as commuting only intro- 
duces lower order terms which have already been estimated. The idea is to 
integrate the first term by parts moving a three derivative on . Modulo 
a lower order term we will obtain [i^^ ^ for which the Bianchi equation yields 

—diva^m ■ Integrating again by parts to move the angular derivative back to the 
Y_ term and commuting it through, we obtain cancellation of the highest order 
terms after plugging in the structure equation (158). This leaves, ignoring the 
boundary terms for the moment the highest order term 

7Pa„,a3„, (306) 



y 

which after another integration by parts is seen to be of lower order. It is easy to 
check that in this process only boundary and lower order terms appear, which 
are collected on the right hand side of the Lemma. □ 

With the lemma at hand we can integrate the identity (301) in y and obtain 
favorable spacetime- and boundary terms if the first term on the right hand side 
of the Lemma 8.4 is taken into account. 



^'*In rewriting these terms like this one picks up lower order terms, which involve (i mul- 
tiplying i-derivatives of a. These are accounted for in the lower order energies on the right 
hand side of the Lemma. 
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Estimating ^g^^^ and (p,cr)3„, 

Using the equations for p^zm and cr33„, (derived from (106) and (111)) and for 
^Mn ('is'^i'^^d from (116)) yields, modulo a total divergence on Sf^) 



2 



'^3 [l{pln.+4n^)]+\lt)Al\\P^j' 



.2 I ^2 > 
■71; ' 3n/ > 



2 

11^3 Jl^ 



37 + {P3.n,)ltrH 



(307) 



(/^'^)3n, '^3„J = ^3 ' (P) • Pan, + ^3 ' (^) ' '^S" 

^3„, ( ^4"' (^) + [204^^^ - ll)^trH] + F43 (^) 



-7C: 



34 



3tI;. 



(308) 



We have the analogue of Lemma 8.4: 
Lemma 8.6. 



(P.cr) 



3ni '^3n, 



< 



7f2 [4 + 2 [number of 3 's m n;) 



— 3ni 



-2 / 7f2 (num&er o/ 4 's m n;) II (p, ct)3jj || +A 
Jy ' Jy 



1^3 J 



(P''^)3n, II 



3n, 



+ error-terms as in Lemma 8.4 



Proof. The proof is completely analogous to the previous one. Hence we will 
only reveal the structure regarding the term which is proportional to p. The 
most difficult term is 



PI 



3n, 



(309) 



First use the structure equation (173), which modulo terms of lower order pro- 
duces a familiar /S B -term which can be integrated by parts as in the previous 

lemma and a term which is like —0 ^^aXJS^^ . Moving the 4-derivative onto the 

/? -term, inserting the Bianchi equations /? = —'^Psni — *'i^cr?,ni +l.o.t. and 

O Tl I O 41: 'Tl' I 

moving the angular derive back onto the Y-term yields a term -1-30"' d^vY ■ p^m 
(note that the cr3„j-term is of lower order in view of the structure equation 
for ci/rlY). For the ^0 '^trH -term wc similarly move the angular deriva- 
tive onto B , insert the Bianchi equation p^ms = ^df^P^ + 1-o.t. and move 

0711 OTLl 

one derivative from P3„,3 back to the trff-term. This leaves the highest order 
term — | j/)"' 0^trH p^^n^ . The sum of the two highest order terms is however of 
lower order by virtue of the structure equation (159). Again due to the good 
main-term in the first line of Lemma 8.6 we obtain non-negative boundary and 
spacetime terms upon integration of (307) in y. □ 
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Estimating jS^m 

Using (101) in conjunction witli (107) and (112) wc obtain 



1 



03 b\\M']+7;^i [7(pL,+0] 



-037 + ^9^ W3n,)ltrH ||/33„,|1 



[P^ni , (p, cr) 



an; J 



(310) 



[/33n, , {p, cr) 



a.3n. [ (a) + 



3n, 



-7C34 [(p, Cr)„J (p, (t)3„^ - 7C33 [/3n,] ^3n, 



+ ^^43 (P) 
- F43 (ff) 



As before we have 
Lemma 8.7. 



[^3n,,(p,cr)3„J < / -fn[2 + 2{number of'S's mni)]\\{p,a)^,^^\\ 



-2 J^-fQ {number o/4's m m) ||/?3„,||2 + A {Wkn.W^ + II (p,f^)3„, IP) 

+A / (||/33„JP + ||(P,^)3,JP) 

+ error-terms as in Lemma 8.4 ■ 

Proof. We only present the control of the p-integral. The crucial term is, after 
using the structure equations to establish that 



^Ip^trH - ^03trH 



+ Fi3 (p) = 



p [Sdiv {Z_ — Z) + 2ilp3 — 2f2p4] + quadr. dec. terms , 



(311) 



1= I IP 
ly 



iip'^'div^z- z) 



P3r, 



-U"'ftrH + -3r'lp,Z]p,^, 



(312) 



Starting with the last term we move the 3-derivative onto to obtain yp3„, + 
*o'3ni +l-o.t using Bianchi. Moving the V derivative back to the left, the cr-term 
is seen to be of lower order in view of the structure equation for ct/rlZ. The 
P3n, term which arises cancels to highest order with the Z-pait of the first term. 
Hence we need to establish a cancellation between the Z-part of the first term 
and the third term. This is achieved by moving the ^ derivative onto iSsm , using 
the Bianchi equation p43n, = div[izni + l.o.t, and moving back the 4-derivative 
onto trH_. The structure equation (162) finally establishes cancellation to the 
highest order. □ 
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Estimating a^m 

Using (99) and (102) 



03 [7ll«3nJP] +04[7ll/33„J|- 

^37 + 79" {a3ni)'-ftrH ||a3„J 



^47 + 2^+ (/33»,)7i'^if 



(313) 



67 [aani , /33riJ = -27C34 [/?„,] /33„, " 7C33 [q;„J a3„, + 7i?3"' (a) ■ asn, 
+27/33„, f (/3) + [/3 [^4irH - 203tri7] + F43 (/3) 



As before we have 
Lemma 8.8. 

67 [a3„,, /33„J < / 27^1 [-2 + 2 + 2 (numfoer 0/3 's m nj)] ||/?3„, IP 

-2 / 7ri (num&er o/4's m 71;) ||a3„J^ + A / (||a3>i, |P + ||/?3n, |P) 
Jy Jy 

+A / (lk3„JP + ||/33„,||') 

+ error-terms as in Lemma 8.4 

Note that the tern 27;33„, (£:4"' (/3)) now makes a negative contribution 
because of the minus sign in E4 (/3) = — 2f2/3 + while both 27/33„,_0"'i^43 (/3) 
and — 27C34 [/?rii]/33ni still contribute a positive il-term. The critical p-term 



PI 



2-30"'03r -/Jan, -30"'03i7a 



3iii 



(314) 



is seen to be of lower order by integrating by first using the structure equation 

03r-04Z = ^ + l.o.t. , (315) 

and then integrating by parts twice (moving 04 onto /33„, and returning the 
angular derivative arising from the Bianchi equations back to Z) and finally 
using the structure equation 



,Z + l.o.t. 



(316) 



The previous for redshift estimates will allow us to control all 3n;-derivatives. 
However, to control all derivatives via the Bianchi equation we need control over 
at least two 47i;-derivatives as well: 
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Estimating /34„, 



^03 bWin,\?]+\lPi [7l|a3n,|| 



1 



^47 + 2z9+ {a3n,)ltrH 



(317) 



^7 [P4n; 1 ^3ni \ 



-7C43 [/3„,] /?4„, - -7C34 [«„,] a3^ 



For the obligatory Lemma regarding tire error-term we can be more naive, since 
we already have a large (proportional to ■^-^) spacetime term containing all 
three-derivatives available: 

Lemma 8.9. 

e7[/?4n,,«3n,] < {2 + 21) \\2^n\\ [ ( || /^a™, f + || /34„ + !1 «3„, f ) 

Jy 

{\\(X3ni\\'^ + 11/^4)1; 11^) + error-temis as in Lemma 8.4 (318) 



(319) 



+A 

Proof. Everything goes through as previously, except that the terms 

-7C43 [/?„,] Pin, + iPin, ■ £^1"' (/?) 

will introduce a highest order mixed term of the form 

7 [-2 - 2 • (number of 4's in n/)] r2/3„,3/3n,4 



(320) 



For this term we use Cauchy's inequality. We leave the considerations for the 
term proportional to p to the reader. □ 

Remark 8.10. The highest order terms on the right hand side of (318) will he 
absorbed by terms on the left once we add the estimates for all the quantities, 
since we control these derivatives already with a largeness factor of ^7^-^ from 
the previous steps. 

Estimating a^m 



^03 [7||a4„JP] +04 [7ll/?4„JP] + [-^037 + ^- {(^in,)ltrH 



+2 



-i047 + 2i?+ {Pin,)ltrH 



(321) 



67 [aim ' f^ini] = -27C43 Q!4„, - 7C44 [/?„,] /34„, 4- 7i?3"' (a) ■ 



which goes with 



llPitrH - ^Ip^trH 



F43 (a) 
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Lemma 8.11. 

67 [Pint ^ 4n/ II 1 1*^471/ I 

y 



+\ / (l|Q^4n, 11^ + j|/34ri, 11^) + crror-terms as in Lemma 8.4 

Once more, we remark that the main terms wiU be absorbed by the left 
hand side, once the derivatives for aU quantities are added. Using the Bianchi 
equations we now prove that the control over the derivatives we have considered 
is sufficient to control all derivatives: 



Lemma 8.12. We have the following estimates on the spheres S^* ^ in the 



¥V C ILUUC Lite JULlULUTUy COl tt/ / ttt tC^ UlL Lf iC ^jJliClt^^y 

region r < 5M : 

K+i ^ K 



K+1 . K I. 

l=Q ni •'S^ 1=0 m •'S^ 

+ 11 (P> '^)3n, f + ll/33n, ll' + W^Zn, f + Win, f + Win, \? 

2 , II oi|2 , II / ~ _\ |2 , II o||2 , II „ ||2 



+B 



|a||^ + ||/3||^ + l|(/5,a)|^ + ||/J||^ + ||a| 



K 

^EE / ll^?"'(5K-5Hss)||' (322) 



1=0 n, 

Proof. Use the commuted Bianchi equations to estimate ||yQi„ |P from ||/3 |p, 

i — ni o 

liV/?^ IP and ||a4„JP from j| {p,cr)^j^^ |p, etc. cf. Lemma 10.6. Finally, use the 

fact that W34 = y^u + l.o.t. holds for the curvature components to estimate the 
remaining angular derivatives from the 34-derivatives. Note that r-weights are 
irrelevant in this region. □ 

We can now complete the Proof of Proposition 8.1. For fixed I add up 
and integrate the equations (301), (307), (310), (313), (317), (321) using the 
Lemmata for the error-terms. Do this for all I and sum over all permutations 
ni to establishe Proposition 8.1, except for additional lower order terms (l- 
derivatives of curvature) appearing on the right hand side (cf. the last two lines 
of Lemma 8.4). Iterate the estimate for these terms, until one finally needs an 
estimate for zero derivatives, at which point we insert the estimate of Proposition 
8.2. 



9 Elliptic Estimates in the interior region 

In the interior region, we can estimate all derivatives provided we control T- 
derivatives: 

Proposition 9.1. Let {TZ,g) be ultimately Schwarzschildean to order k+l with 
k > 7. For n > 1 we have 

\V"W\'^r'^drduj < bV] [ \CrrW\^r'^drduj + B [ \W\'^r^drduj 



En 



+ B-K" '[^](Sr) 
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Moreover, 



'^TY- '^ <r<R+M 



[W] (tW (ri,T2)) < B ■ D"-i m (ti,T2) + 
B f l^Iy (i-^sA' \C\.W\^ + \W\^}r^dt*drdLu 



-n.nondeq 

I , 

r-s 



Here the constants B depend on the mass M and lower order energies of the 
Ricci- coefficients and Weyl- curvature, which are bounded (or decaying) by the 
ultimately Schwarzschildean assumption. 

Proof. (Sketch.) Write the Bianchi equation as a div-curl-sysicm for the electric 
and magnetic part of the Weyl-tensor, the latter defined with respect to the 
timelike normal to E, n-^. Apply standard estimates from [8]. The (quadratic) 
inhomogeneity is easily estimated by lower order energies of the Ricci-coefficients 
and Weyl-curvature (using Sobolev embedding). □ 



10 The region near infinity 

In section 8 we revealed a remarkable hierarchy in the Bianchi equations, which 
allowed us to obtain estimates for all derivatives of curvature. This hierarchy is 
also present near infinity and somewhat even more remarkable. This is because 
close to the horizon all curvature components decay at the same rate in r (there 
is merely a difference in the strength of the redshift factor), while at infinity 
the curvature components and their derivatives each have a characteristic decay 
in r, which the estimates need to reveal. It turns out that using appropriate 
rP-weighted multipliers (for some positive p) at the level of the nuU-Bianchi 
equations is sufficient to establish the correct asymptotics. In principle, just 
as near the horizon, it suffices to commute the Bianchi equations with 3— and 
4--dcrivatives (and obtain angular derivatives from the wave equations satisfied 
by the components). However, for the optimal r- weights (as exhibited in the 
energies of section 4), this commutation seems unavoidable. Since every set of 
Bianchi equations will admit a different r- weighted estimates, depending on how 
many derivatives have been taken, we have introduced the concept of boundary 
admissible tuples and matrices in section 4.2, which should be understood as 
encoding the r^ weighted decay for the individual components at each order. 
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10.1 The main result 

Theorem 10.1. For any boundary admissible matrix P we have for < n < k 
Ep [W] {Nout (t2, i?)) + [W] {V (ri, T2)) < Ip [W] {N^ut {ti,R)) 

+B ■ r/_Z^r<R+M m [M (n, r^)) 

+B 



sup E [m]{Nout{T,R))+l [5H](2?(ti,T2)) 

Te(Ti,T2) 



(323) 



The constant B in the second line depends on M and on lower order curvature 
energies which are bounded by the ultimately Schwarzschildean assumption. 

10.2 Estimates for the uncommuted equations 

Multiplying (85) by rP2-2 _ ^ j-^j. positive integers p2, q produces 

1 



^3 [rP^-^l - ti)' \\a\\ 



P2 



—tr (H)r qr 

2 2^ 



1 „P2-3 



= _2rP^-2 (1 - ;i)9 ^1/3-0 + rP^--^ (1 - /i)" {E3 (a) - Ana) 

Using that div is the adjoint of 'P2 ^^'^ inserting the Bianchi equation (86), we 
can write this as 



1 



Da ((63)" rP^-^ (1 - m)' M') + Da [(64)"^ rP^-^ (1 - m)' II/3|| 



+ \\a\\\P'-^l - fi)" 

+ ||/?||V^-3(l-/i)'? 



^^2-2 1 03(1-Ai) 

— SUr — —rq 



2 " - 2 

3r trH + dQr - {p2 - 2) $^r - qr 



04 (1 - M) 



(1 - ^)« ( (£;3 (a) - ma) ■ a + 2 {E^ (p) + 2VLf3) • ^ + 2 [div [ap)] ) 

We integrate over the region taking into account the following observations: 

• the terms in the last line are error-terms. For the last term we integrate 
the angular derivative by parts to produce a cubic error-term in view 
of the fact that both r and the measure approach something spherically 
symmetric. 

• The first square bracket is equal to \{jp2 — 2) + O 

• The second square bracket is equal to [(8 — P2) + O (i)] . However, choos- 
ing g to be a sufficiently large negative number [q = —4 is good enough), 
we can ensure that the O (■i)-term has a positive sign. 
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Using that 1 — yu w 1 in this region and that the measure brings in another 
power of 2 in r, it follows that for sufficiently large R we have the estimate 



dvduj2 rP^\\af + 



X-2 



-B ■ 



^O.deg 

^R-M<r<R+M 



[W] {M (ri,r2)) 



dt* drdij < B 



dudu}2 rP^\\l3\\^ 
dt^ dr duj 
ep2 f3] 



dvduj2rP^\\af , (324) 



with the error 

ep, [a, /3] = rP^-~^ (1 - ^if ( [E^i (a) - 411a) ■a + 2{Ei (/3) + 2fl/3) • /? + 2 [d/w (a/3)] 

It follows that for 2 < p2 < 8 we obtain good terms on the left hand side, 
provided R is chosen sufficiently large. However, there is a non-linear restriction 
from the error-terms in the penultimate line as we will see in the next section. 
Let us define the renormalized quantity p = p + ^ . Multiplying (89) by 
(1 — p.)"^ P and integrating by parts as before to insert (91) and (92), yields 
for sufficiently large R the estimate 



dvduj2rP^\l3\'^ + _ duduj2rP' {p^ + a^) 



T2-R 



dt*drduj\B\'^rP^-^ 



(P3 - 4) + O I i 



/ dfdrduj (f + (T^) rP^-i [(6 ~p3)+p]<B ( 



-B 



Tf{].deg 

^R-M<r<R+M 



[W] {M (Tl,T2))+4 



[/3, (p,cr)] 
dvduj2 r^Hp? 



(325) 



with the error 



^P3 



[/?, (p, a)] - _ _ 2^1/3) ■P + E^(p)p + Ei (a) a + div a) /3] 



hence (325) produces a good estimate for 4 < p3 < 6. However, note that 
P3 < 4 is also admissible as long as < P27 a-s we can estimate the wrong 
signed /3-term by adding a multiple of the previous estimate (324). 
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We now consider the underlined quantities. We start with equations (93), 
(94) and (90), from which we derive 



dv du2 r^"" {p^ + cr^ 
dt*drduj (p^ +cr2) rP"-! 
+ / dfdrdu r'P"-^ [(4 - p4) +tA< 
[W] {M (Ti,r2))+4 



-° ' ^R-M<r<R+M 



duduj2rP^W+ 
. - 6) + O 

dvdL02r^^ [p^ +Sr'^) (326) 



with the error 

[(P, ,^ = (1 - m)' [ (^4 (^) - .^+E3{p)p + Ei (a) a - div [/3 • (p, a) 

Obviously, the first spacctime term demands p4 > 6 while the /3-tcrm re- 
quires p4 < 4 in order to be positive. However, since we already control the 
spacctime integral of (p^ + tr^) from the previous estimate, the restriction is 
Pi < min(p3,4) only. 

Finally, we multiply (88) by r^^ /3 (1 — /i)° and integrate in the region to 
arrive at the estimate 



dvduj2 rP^^ I^P + 
dt*drduj IBl^rP''-^ 



T2-R 



du duj2 rP'-" 



(P5 - 8) + O - 



-B 



+ 

M,deg 

^R-M<r<R+M 



dt*drdLo\a\^rP'-^ 



(2 - ps) + M 



< 



[M^] (7W (ti,T2)) + 4 



dwdw2r^=|/3p (327) 



with the error 

a] = IrP^'^ (1 - m)' ( (^^3 {§) + 2n§) • ^ + ^ (^4 («) - 417a) • a - 2dh {a ■ §) 

Wc conclude that this estimate requires p < 2 and p > 8 to make both volume 
terms positive. However, by the usual argument of adding multiples of the 
previous estimates, only p < min (2,^3) is relevant. Wc summarize this as 

Proposition 10.2. For a boundary Q- admissible tuple (0,p2,P3,P4,P5j 0, 0) we 
have, for sufficiently large R, 



+1 



%P2,P3.P4,P5,0,0) [W] {Nout {t2,R)) 

m {V{n,T2)) 



(0,P2-1,P2-1,P3-1,P4-1, 

<E(o,p,,p3.p4,P5,o,o) [W] {Nout {ri,R)) 
+B ■ fRtZ<r<R+M m {M (n, T2)) + B ■ i?rr(o,p„p3.p„p„o,o) (328) 
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whe 



Err 



(0,P2,P3,P4,P5,0,0) 



dt*drdu! [a, /3] 



V' 



dt*drdu! [/?, (p, cr)] 



V' 



dt*drdu! Bp^ [{p, cr) , 0j 



dt*drdujep^ [^,a^ 



Proof. Add up the estimates (324), (327), (325) and (326). 



□ 



For this proposition to be useful, we need to estimate the error-term in 
terms of the bulk-term on the left hand side. Clearly, this may impose further 
constraints on the admissible p's. We now show it doesn't: 

Proposition 10.3. For a boundary 0-admissible tuple (0,p2,?'3,P4,?'5, 0, 0) we 
have 

£'^^(0,P2,P3,P4,P5,0,0) < 



Wo 



Jl(0,p2-l,p2-l,P3-l,P4-l 



^ ■ Vp2-i,p2-i,p3-i,P4-i,p5-i,o) [^] P(n,r2)) . 

Note that in view of the smallness assumptions on the Ricci rotation coeffi- 
cients given in Definition 2.2, the estimate of Proposition 10.2 closes. 

Proof. For the a-part of the error of estimate (324) we have 

J rP' {E3 (a) - iQa) a 



< J j rP^~^\a\^dt*drdujJ J rP^+ij {E3 (a) - AQ) \^dt*drduj . 



(329) 



We note that 



rf^+i| {E3 (a) ~ An) \'dt*drduj 



< C 



\\Hf{p'+-') + {\\Vf + \\Z\f) 



rP^+^dt*drdu; 



<C||pV||l- / \\H\\\P'-^dt*drdoj + C\\rm\V\\^ + \\Zf)\\L^ rP- 



-3|«|2 



dt*drduj 



from which we obtain the constraint p2 < 7 — d (otherwise the spacetime-term 
involving is not small). 

For the /3-part, recalling that in our frame Y = 0, 



rP^+^\ {E^ (/3) + 2np) fdt*drduj <C {\\Vf + \\Zf) \a\^rP''+^dt''drdLO 



<C\\r^ {V^ + \\Zf) 



rP^-^W^dfdrduj 



which does not impose further constraints on p2. 
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Turning to [/3, (p, a)] we have 



rP^H ■ P dt*drduj < y J rP^-mi3pdt*drduj^ WH^t^Wl'^ / r2P3-P2-3||^||2^i-*.^^^^^ 
from which we read off the constraint 2p3 — P2 ^ 3 < P4 ^ 1, 



y rP^Y-a-pdt*drduj rP^--mi3\\^dt*drduj^ WV^v'^Wl^ J r^P^^-P^^^WaW^dt^drdu 

from which we read off the constraint ps < P2 + 2, 



rP'^pZ ■ (3dt*drduj < W / rP^-i||/3||2di*drda;J / r^P-^-P-^-'^\\ZW^dt*drduj 



from which we read off the constraint 2p3 < p2 + 7 — 6 in order for the ||Z|p 
spacctime integral to produce a smallncss factor. A similar estimate is obtained 
for the term *Za. Turning to £'4 (p) we estimate 



J rP^'H- apdt*drdu} < ]J J rP3-'^\\p\\'^dt*drduj^ \\H^r'^\\L^ J rP3-'^\\a\\'^dt*drdu 

which imposes p3 < p2- The remaining terms do not create new constraints and 
we can turn to e^^ [{p, a) , 0j, for which the most difficult terms are 



J rP^H ■ apdt*drdLO < ^ J rP^-^\p\\^dt*drduj^^\\mr^L-=- j r^P-^-P^-^aW^dt^drduj 
imposing 2p4 — < + 2 and 



rP^pZ- ^dt*drduj < ^ J rP^-mi3\\^dt*drduj^ Wp^r^WL-- / rPi~^\\Z\\^dt*drduj 
is fine since p4 < 4 anyway. Finally, for e^^ a] we note 



rP'- {V,Z)-apdt*drduj rP^-^\P\\^dt*drduj ^^\\ {V,Zfr^\\L^ J r^P^-P^-3\\aPdt*drduj 

revealing p^ < 2 + p^, which is already satisfied as < 2. Inspecting all the 
constraints on the p's obtained, we observe that they are all implied by the 
assumption that the p-tuple is boundary admissible. □ 

This establishes Theorem 10.1 for n — 0. 

10.3 Estimates for the commuted equations 

To prove Theorem 10.1 we will first prove the higher derivative version of Propo- 
sition 10.2. The focus here is on the main-terms and their r-weights, while all 
errors are simply collected on the right hand side. The latter will be estimated 
separately in Proposition 10.9 in the next section. Combining the two Proposi- 
tion will imply Theorem 10.1. 
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Proposition 10.4. For P a boundary admissible matrix with second row being 
the tuple p ~ (pi,p2, •■•,P6j 0) we have the estimate 

fT;''' [W] {Nout (t2, i?)) + r^'^' [W] (V (n, T2)) < lp+' [W] {Nout (ri, i?)) 



-° ' %-A/<r<_R+A/ 



[W] (tW (ri,r2) 



B ■ Errp^^ + B ■ Low'^+^ 



(330) 



Here the error is given by (all integrals are taken over 'DJ^^ ) 



k j 



j=0 j=0 1=0 rii 



ep3+2i2 



'Pinn^i ' (P4n, , Cr4n, ) + / ^^2+2(2 ^1 '/34n,,^i '"Sn, 



+ 



ep4+2i2 



6p6+2i2 



Low™ = Low^p^ ^ / dvdur'^^^LoWnl {P2,P3,P4,P5,P6) 

j=0 j=0 1=0 n, •'^out 

+ X! X] X! / ^'"^^ LoWni iPl,P2,P3,P4,P5) 

j=0 1=0 n, "'■^n 

where ni denotes a l-tuple of 3 's ami 4 's wzt/i 6em(; the number of 3 's and I2 
being the number of 4 's. T/ie integral is short hand for J = J ^/gdt*drduj and 
the integration region is always • The precise definitions of the [..., ...] are 

given in the following section, while LoWnl (•■•) '■^ defined in Lemma 10.6. 
10.3.1 Derivatives in the 4-direction 

We start with the Bianchi equations for a34„, (derived from (99)) and /344,i, 
(derived from (102)). First we relate a34„, to a4„,3 using the commutation 
formulae. Since the resulting equation is of the same form as the original Bianchi 
equation, we proceed as for the non-commuted equations, i.e. we multiply the 
equation for q;4„,3 by r^"-^ (1 — oiAni and integrate by parts. This yields the 
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estimate 



dvdw2 r^||a4ni|| 



dudui2 r^llAn, I 



+ / dt*drdLj\\ainif rP-^ 



PI? 



j03r + 7?- (a4„Jtr {H)r + o{- 



+ / dt*drduj ||^4n, iPr-P-i 2^9+ (^4„,) r-trH- plp^r + fi 



< B 



Sp [aim , P4n,\ 



+ S-l'+iit|,.<fl+M [W] {M {n,T2)) 

+2 [ dvdiJ2rP\\a4n,\\^ (331) 

with the error-term 

ep [a4„,,/34n,] = rP'^ (1 - m)" (4^"' (a) «4n, + ^E^' (/3) • /34„, + 2div [a4„, • /?4n,] 
One reads off that the left hand side is positive for r > i? sufficiently large if 

2(4 + ;-s(a4„,)) <p<2(4 + / + s(/?4„,)) (332) 
holds, which is equivalent to 

2 + 4 (number of 3's in ni) < p < 12 + 4 (number of 4's in n;) . (333) 

From the Bianchi equations for /334„, (derived from (101)) and (p, o')44„j (derived 
from (107), (112)) we obtain, after relating /334„, to /34„,3: 



1 



dv du2 rP 1 1 Pin: 11^+2 / X, - '^'^2 {Plni + O'li, ) 



+ / dt*drdLjrP-'\\/34n,f -^lp3r + ^^ (MrtrH + O 



i / ^ dt'^drduj rP-' {pI^^ + <, ) [-p04'^ + 2^+ (P4„, ) r iri/ + /z] 



< 



e-p [Paui, (P4n,, tT4„,)] 



s - r 



k^^M^c^ii+MlW^] iM (ri,r,)) 

d«dc^2r''||/34„J|' (334) 

with the error-term 

ep [/34„,, (P4„,, fT4„,)] = rP-^ (1 - //)'' (4""' (/3) • /34„, + 4"' (p) ' P4„, 

+ d/z; [( 

The left hand side of (334) is positive for r > i? sufficiently large if 

4 + 4 (number of 3's in ni) < p < 10 + 4 (muirber of 4's in n;) . (335) 



78 



Note that only the upper bound is relevant, as we already control the a4„, 
spacetime term from the previous step. 

Thirdly, combining the Bianchi equations for (p, cr)^^^^ and P^^^ 

out \ T2 . R J Ti— R. 



+ / dt*drdurP-^ (pL, +^L, 
1 



2,^2 



+ - dt*drdLorP-^\\B IP [p04r + 2i9+ \rtrH + , 



< 



dvduJ2 (pL, +cr|„,) 



(336) 



4ni 



with error 

and a good left hand side for 

6 + 4 (number of 3's in ri/) < p < 8 + 4 (number of 4's in ri/) . (337) 
Finally, using ^^^^^ and 044^, 



2 It^2-R 



du dido r^WctA 



^Jt^drdu^WP^J^r 



-J 
2' 



+ \ j^^^ dedrdurP-^\a,^^ f [-p$^r + 2^+ {ciin,)rtrH + p] < 

+2 



2,^2 



dvd^2rP\\P^J' (338) 



with error 



rP-^ (1 - p)'^ ( i^s^"' {§) ■ + i^^"' {a) ■ «4„, - 2d^ 



and a good left hand side for 

8 + 4 (number of 3's in n;) < p < 6 + 4 (number of 4's in ni) . (339) 
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10.3.2 The missing derivatives in the 3-direction 

Controlling all 4-derivativcs is not quite sufficient to control all derivatives from 
the Bianchi equations: We will also need and a-i^ . The first can be 

obtained by starting from the Bianchi equation for 043^^ and then using 13 , 
which yields 



( dvdu2rP\\a^^f + ( du(hj2r-P\\l3 



■D ^ 



p 



2'^ ' ~ \—3ni 

dt*drdLO r^'^Wi^^^W^ [2^9+ (/^^__ ) r trH - plp^r + fi 



< 



V' 



.3n, 



dvdu:2rP\\a,J\' (340) 



with error 

and positive left hand side for 

10 + 4 (number of 3's in ni) < p < 4 + 4 (number of 4's in ni) . (341) 

Note however that only the upper bound is going to be relevant in view of the 
previous estimate (338) for the spacetime term of ||a4„Jp. Finally, from the 
Bianchi equations for P^^^ and a^:}ni 



-l-R 



+ dt*drdio\\l3\\^rP-' -Ljp^r + d-{l3\tr{H)r + 0(- 



— 3ni 

2i?+ ("an, )rtrH - p$^r + ^^ 



< 



rP\\P^J'dvdL02 (342) 



ep [^3„, , a3„J = (1 - m)' (^4'"' (^) ' ^3n, + 2^'^' ' ^3™, " div 
with the left hand side being positive for r > i? sufficiently large if 

12 + 4 (number of 3's in n;) < p < 2 + 4 (number of 4's in n;) . (343) 
As usual, in view of (340) only the upper bound is going to be relevant. 
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10.3.3 Angular momentum operators 



The derivatives obtained so far are in principle sufficient to control all derivatives 
using the Bianchi equations only, just as we saw in section 8. However, for the 
optimal 7'-wcights, we need to commute with angular momentum operators as 
well. 

We first note that we can derive the identical 4+2 estimates of sections 10.3.1 
and 10.3.2 for the fii-commuted quantities. We will see the same constraints on 
the admissible p's (as commutation with the fi^ doesn't change the signature of 
the curvature components, cf. section 3.3), while the errorterms are replaced by 

ep [^]^«4„,,^^^'/34„,] = rP(Ef''^^ {a) ■ nja^n, + 2Ef''^^ (/?) • r!^'/34„, 

+2div Vtlaini -^iPini ) 

and analogously for the other components. We also note 

Lemma 10.5. Given a covariant tensor f tangent to the two-spheres S'ft, „ on 
TZ there exists a constant cq close to 1 such that 



1 

Co Js^ 



f\^dfi^<co / {\f\' + r^\ff\')df,. 



(344) 



where we denote |/si/P = J2i l-^f2i/P- Moreover, for f a scalar we have 

-[ r^\ff\^dfi^< [ \l^f\^d^^, < Co [ (r2|y/|2)d/,^ 
Co Js^ Js^ Js'^ 



(345) 



while for f a one-form or a symmetric two-covariant traceless tensor tangent to 
the surfaces S'^ 

(346) 



"0 JS2 JS2 



Co Js 

This lemma allows us to gain improved decay in r for the ^-derivatives as 
exhibited by the E [l^]-energies (cf. (182)). 



10.3.4 Retrieving the remaining derivatives 

Controlling the derivatives appearing in the estimates of sections 10.3.1 and 
10.3.2 suffices to control all remaining derivatives from the Bianchi equation 
(Cf. Lemma 8.12): 

Lemma 10.6. For j > 

B I (r'MI/a.«4„,f + r^=||/'n./34„,f + r'^||/^,. {p.a),^^ f + r'^'WCn.P^J? 

+r'''Uk^4nf +r''''Uk^,j' +r'''U'^^^^^^ +B f^Low^j {q2....,q,) 

•J S" 



81 



with 

LouS! fe, ...,96) = [r'1|£;r"' (/3) r + r'^' (pr""* («) f + llC' (p) f ) 

(|1£;3"'^^- (/3) iP + (^)ll 
(lli^-"' (p,a) IP + ^Ef'' (a) f ) + (^) | 

Proof. This follows by integrating the commuted Bianchi equations over and 
integrating by parts. For instance, in the case j = 0, from the bound on p4„, : 

52 JS2 



</ \\V^|3„,r< |l(p,a)4,J|^+ / ||i?r'(p)||^ + ||£;r(a)||^ (347) 

JS2 Js2 Js2 

where K is the Gauss curvature of the S^. This in turn controls 
Ik3„,f <8 / + 2 / |li?3'(«)lP 

< 8 / PnP2V*,Pn, + 2 / 1|£;3"' («) II' + 2 / ||i?^'' (a) f 



<-4/ /3„,4^/?„, -4K||/3„J|2 + 2 / ||£;^M")ir 

<4 / II^AJP + 2 / \\E^^{a)r (348) 

In other words, 

/ ||a3„J|2 + ||y/?„J|2< 

8( /, II (p.-^)*,,, ir + 11^^?' (p) ir + m (^) ii' + iii?,-' w ii') (349) 

5'"' 

Similarly, we estimate the pair y [pnn o'm) and /33„, from , the pair '^(3 and 
(/3n,3,o-ri,3) from o^^^ , and finally 

JS2 JS2 

</ \\'D2aJ\'<[ \\^J'+ f jE'S'iPjW (350) 



□ 



10.3.5 The summed asymptotic estimate 



The claim is that we can add up the asymptotic estimates to estimate all deriva- 
tives at order k + 1 with the weights appearing in the E-energy. Before we do 
that let us note the following lemma, which establishes that permutations of the 
ni are lower order. 



82 



Lemma 10.7. Let u S {a, /3, p, cr, /3, a} he a curvature component and Up a 
p-tuple of 3 's and 4 's. We have 

np-2 



where the sum is over all p — 2 tuples containing one less 3 and one less 4 then 
the original Up. In addition 




Proof. The second statement is immediate from the definition. The first state- 
ment follows by induction using Lemma 3.5 and the pointwise decay of the 
Ricci-coefficients. □ 

Proposition 10.4 now follows easily from the next Proposition, which reduces 
the problem at order ^ + 1 to the problem at order I. 

Proposition 10.8. Let p ~ (pi, ...,p6i 0) be a boundary admissible tuple with 
associated bulk admissible tuple p ~ (pi, ...,p-j). For I > we have 

[W] {N^ut (Sl^n)) + 4^''"^' m {Vll) < B ■ E^+i [Vt^] (iV„„, (5^^.^)) 
+B ■ Ir1mZ.<r+m m {M (n, r2)) + B ■ Err'+' + B ■ Low'+' 

+B ^ (e; [w] {N^ut {sI^r)) + 1;'"^^^ [w] {v;i)) . 

i=l 

We remark that the last line estimates lower order terms arising from com- 
mutation of the ordering of derivatives. It is absent for / = 0. 

Proof. We start with the tuple of length I (denoted ni, with li being the number 
of 3's, I2 the number of 4's) consisting of all 4's (i.e. = and I2 = I) and apply 
the estimates as follows: 

• (331) with p ~ pi + 2I2 

• (334) with p = p2 + 2I2 

• (336)with p ~ ps + 2I2 and - in case that pa — S and Z2 = - in addition 
with p ~ 8 — 6 

• (338) with p ~ p4^ + 2I2 and - in case that p^ = 6 and ^2 = - in addition 
with p ^ 6 — 6 

• (340) with p = P5 + 2I2 and ~ in case that = 4 and ^2 = - in addition 
with p — A — 6 

• (342) applied with p = pe + 2^2 and - in case that pq — 2 and ^2 = - in 
addition with p ^ 2 — 6 
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Remark: The reason for the additional apphcation with a (5-loss in the last 
four estimates is caused by the fact that, as we saw, in this "extreme" case, 
the spacetime term loses an additional power compared to the boundary term, 
i.e. instead of the usual weight pa — I for the spacetime term, it will only ad- 
mit p4 — 2. Applying the estimate also with a delta loss ensures that we get 
the spacetime term with an improved weight of p4 — 1 — 5. (We also obtain 
weaker r^~''-weighted positive boundary terms in this process, which we simply 
discard.) 

We next add the estimates above, so that all spacetime and all boundary 
terms on the left hand side are positive: This is already automatic for the first 
three estimates in view of the constraints (333), (335), (337) being satisfied. 
For estimate (338) we observe that the possibly negative (if ni is the zero tuple) 
signed spacetime term containing IP can be absorbed by adding enough of 
the previous estimate (336) which controls this term. This uses that the tuple 
is admissible: The spacetime term of the second quantity of the i*'' estimate has 
weight 7-Pi+2'2-i jj-^ ^jjQ extreme case) while the spacetime-term of the 

first quantity in the {i + lf^ estimate has weight 7.Pi+i+2'2-i_ Since pi > pi+i 
holds for an admissible tuple {pi > pi+i + 2 in the extremal case), one can always 
control the first spacetime term of the (i + 1)*'' estimate from the second of the 

Next, an application of Lemma 10.6 allows us to control ||Dm„JP where u 
is any curvature component, D is either a 3-, 4- or an angular derivative and ni 
is the tuple consisting of only 4-derivatives. In particular, the spacetime term 
of ||a34...4|P is controlled with the same weight as ||p44...4|P, namely 7-^2+^'"^. 
Modulo lower order terms, which are covered by Lemma 10.7, this is equivalent 
to controlling ||a4fi, |p for tuples hi which have I2 = I ~ 1 and li = 1. Hence we 
can iterate the procedure and start with the estimate (331) again. Now only 
the upper bound in (333) is relevant, as we already control the spacetime-term 
of |ja4ft, p from the previous step. Since for an admissible tuple pi -\- 2I2 — 1 < 
P2+2 {I2 + 1) — 1, we can apply (331) withp = P1+2I2 and then step down: (334) 
is applied with p = p2 + 2I2, (336) applied with p ^ p^ + 2I2, (338) applied with 
p = P4 + 2I2, (340) applied withp = P5 + 2I2 and (342) applied with p = pQ + 2l2- 
Moreover, we can add these estimates so that all spacetime and boundary terms 
are positive. Using once again Lemma 10.6, we now control all derivatives of 
the form Du„^ , where now ni is a tuple consisting of cither all 4's, or all but one 
entries being 4's. In particular, we control the spacetime term of ||Q!334...4p, 
which is, modulo lower order terms, equivalent to controlling ||Q;4ftJp for tuples 
hi satisfying I2 = I — 2, li = 2. We hence reiterate our estimates.... 

Using this stepping down procedure we will eventually control all derivatives 
of the form Dum where n; is any tuple of 3's and 4's. At this step one can 
bring in the wave character of the Bianchi equations (M34 — = l.o.t holds 
for any curvature component u) to estimate all derivatives. However, there 
is a drawback: This procedure will only improve the r-weight in the energy 
by a power of one per angular derivative, since a 34 derivative pair gains only 
j,o+2 terms of weights. To prove the full decay, one has to bring in angular 
momentum operators: We apply the estimates of section 10.3.1 and 10.3.2 to 
the rii commuted equations (for which they are valid by the remarks of section 
10.3.3). More precisely, to any arbitrary fixed tuple of length j we apply l—j fli- 
derivatives and redo the algorithm we outlined above. This produces precisely 
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the error-terms collected in Err^'p^ . □ 
10.4 Controlling the error-terms 

We now prove the analogue of the error-estimate of Proposition 10.3: 

Proposition 10.9. For a boundary- admissible matrix P with second line being 
the tuple p — (pi,p2 7 •■•,P6jO) we have, for any A > and m > 

W^T' < A • I^'^' [W] {V (n, r^)) + B ■ fn-M<r<R+M [W] {m (n, r^)' 
+Ba sup r"+' m (Nout (t, R)) + [m] {V (n, rs)) 

TG(Ti,r2) 



(351) 



and 



maxil^rti) 



+B ■ (1" m {Nout {Sin)) m (Kl)) , (352) 

where in case that m = 0, i?o = e_R a small constant (which gets smaller 
as R is chosen larger). Generally, the constants B depend on the mass and 
lower order energies of the Wey I- curvature, which are bounded by the ultimately 
Schwarzschildean property. 

Before we turn to the proof let us outline the main idea. Inspecting the 
structure of Err^^^ we see that the error-terms consist of derivatives of the 
inhomogeneities £'3^4 (a, ...,a) in the Bianchi equations. As is manifest in our 
energies used, taking a 4-derivative or an angular derivative improves the (point- 
wise) decay in r by one, while a 3-derivative does not change the r-weight. Hence 
the higher derivative estimates for the error would be a completely straightfor- 
ward generalization of Proposition 10.3, if it was not for the fact that we are 
actually claiming stronger decay in r for some higher derivatives than the naive 
improvement above provides. For instance, we applied the estimate (324) with 
P2 = 7 — S, while we applied (331) with p = 10 — 6 instead of the naive 9 — S 
which would follow from the above reasoning. In order for this to work, there 
has to be a cancellation of the slowest decaying terms, which we are going to 
unravel in the proof. Note that this special improvement only occurs from the 
0*^ to the 1^* order, as for the second order, we claim no additional improve- 
ment. This is immediate from the structure of the decay matrices: All tuples 
following the second line are all identical to the second. 

The cancellations are apparent by doing the following computation: 



Lemma 10.10. In the standard null-frame for which Y ~ , we have the 
following expression for the error-terms 

Et (a) + 4 {$4^ a - AHp^ - A*H(Ji + HE^ [p] + *HEi {a) 

+2?(i^4 iP)) - 3/9 (-2nH - a) + (F + 4Z) i/34 

-2i'trH + ^trH {V + 4Z) + ^4 (y + 4Z) j §/3 (353) 
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El (a) = El (a) + F34 (a) + a trHtrH + ntrH - hVLtrH - Ap 

+2H ■ H + divZ - hdivZ + Z_- ^-bZ ■ z\ (354) 



El (/3) ^2H-(3_^ + 2nil^ + Ya4 + 3Zp4 + 3*Zct4 + 



hY - -trHY 
4_ 2 — 



-Ann - 2a 



)4*Z - -trH*Z 



a + *F3 (a) 



^^Z - -trHZ 



P + F3 {p) - ^ftrHp - Hfp - ^*ftrHa - *Hfa 



El {(3) = -2VL(3i - (2^4 f7 + 2ntrH) /3 + alp^ {2V + Z) + {2V + Z) 

+ ^trH {2V + Z)a- ^^trH -a- H ■fa + trF^ [a) (355) 



El{p) 



-Ha, A 



2nH ~ a) a + P ■ IP4 (y ~ 2Z) + 2P ■ IP4Y 



{V - 2Z) ■ ^ + 2Yj3i + -trH {V -2Z)- trHY/B 

+trHP ■ ^trH + §_ - trF^ (^) (356) 



Et (p) = -\Ha4 - ^IPiK +{V + 2Z) ■ /34 



04 {V + 2Z) + -{V + 2Z) 



P - 20trH -H-fl3 + trFi (/3) (357) 



El {0) = 2np^ + 2HJ34 + 2j3U)iK + 2^ {ip^^ + trHn) 
-3Zp4 - 3*Zcr4 - ip$iZ_ - 3(704* Z + itriJ {-iZp - 3*Z_a) 

+ ^pf (trH) - '^a*f (trH) + H ■ f p - * H ■ fa - F^ (p) - * F4 {a) (358) 

El (a) = 4^04 - 3^p4 - THoi + {V - AZ) - iplJ)Ji - 2,(t$4*H_ 
+a {A^^Vl + AtrHn) + [04 {V - AZ)] - ^trH {V - AZ) ®P_ 

Hdiv§_ - *Hci/rl§_ + {"ftrH) - 2s (F4 {§) ) (359) 



El (a) = Aila^ - SHps - 3" Has + {V - AZ) i/3 + Aalps^ - 2s {F3 {/3)) 



Ip3*K + trH*H_ 



hK + tr HH 



h [y - 42;) + -trH {V - AZ) 



p + H_div§_ - *Hc4rl§_ + 2 l^trg) ®§_ 

(360) 



86 



Proof. A lengthy but straightforward computation. □ 

Remark 10.11. The remarkable thing is that inserting the structure equation 
for H , (171 ), which in our gauge (Y = reads 

l^iHAB + trHHAB = -2VlHab - a as , (361) 

always cancels the terms with the worst decay. For instance, E3 (p) only decays 
like r~^ pointwise, while (p) decays like r~^ , while the naive improvement 
would only suggest r^^ . Similarly, while E^{a) decays only like r^^ (the pH- 
term), E^ (a) decays like r^~ . 

Proof of Proposition 10.9. We turn to Proposition 10.4 to estimate the expres- 
sion for Err^^^ . For j < m and some < I < j , the first term we have to 
estimate is (recall that I2 denotes the number of 4's in a tuple of lenght /) 



We conclude that it suffices to understand the case j = 0,1 = 0. This is a 
consequence of the structure of the errorterms E, E, E and the fact that we only 
have to worry about the number of derivatives involved and about the amount 
of decay in r. Since the E are quadratic in Ricci-coefficients and curvature, the 
number of derivatives is not a problem. Moreover, since taking a 4-derivative 
improves the r-dccay of both Ricci coefficients and curvature components by a 
power of 2 in the energy (i.e. wherever the derivative falls on), while taking a 
three or an fJj-derivative does not change the decay, the additional weight factor 
of r'^ in (362) for each additional 4-derivative is naturally incorporated. 

We hence turn to the case j — 0,1 ~ 0, starting with the first of 7 terms, 
Cpi [a4nfc , /34Tifc ] , which is in turn is written out explicitly below (331). For a 
boundary admissible matrix P with second row p = (pi,p2, ■•••,P6j 0) and as- 
sociated bulk admissible matrix P with second row p ^ {pi ~ l,pi — l,p2 — 
1, ....,P5 - l,pe - 1): 

\\Et (a) f + \\F3i (a) f ) rfi+i < \r^if\ f dt^drduo ||a4|| ^p^-^ 
V Jv 

+ \\r' (h, V, ...)' ^ dt*drdu:[{m\^ + | {p,a)^ ^ r^^-^ 



dt*drdu\\HPrP'~^ 



V 



+ (e, [pf + Ei {<jf) + IIp^oVI^IUo 

+ \\p^r^\\L^ / dt*drduj\\afrP^-'^+ 
Jv 

+ lk^l/3nU~ / dt*drdujrP^-^ {Wl/trHW" + \\ll)i{V + AZ)f) 
Jv 

< ER ■ fp^ [W] {V) + e- 1 [% (D)+B- 1" [% (D) 

the last step following as long as pi < 10 — 5. Remarkably, the last term in the 
expression 

' Et [a) ■ a^rP^ = f {E^ [a] + F34 (a) + a [...]) • a^r^^ , (363) 
V Jv 
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where [...] denotes the square bracket in expression (354), actuaUy generates a 
good lower order term. The point is that the leading order term in [...] is —-^ 
at infinity (arising from the trHtrH_-teim), while all the others decay like 
and are easily estimated as before. We write 

a [...] • a, = \U), {r^^ Waf [...]) + Waf ([...] r^^) + ^vH [...] r 

and find, using the structure equations and the pointwise estimates for the 
Ricci-coefficients that a negative boundary and spacetime-term is generated, as 
long as pi < 12 (which consequently have a good sign when brought to the left 
hand side). Besides an improved weight for a on I, we gain a spacetime term 
/ dt* drdujr'P^~^\\a\\'^ , which in view of pi < 10 — 5 considerably improves the 
estimate from the uncommuted equations (where we obtained only control over 
/ dt*drdu)r^~^\\a\\'^)?^ We summarize 

Et {a) ■ a^rP' < -b f dt*drdujrP'-^af + B ■ fnt^KrKR+M [W] {M (n, T2)) 



en ■ li [W] (P) + e • I [5H] {V) + B ■ 1° [fH] (V) 
Similarly (note Ej (/?) = Ef (/?), in view of C44 [f3] = 0), 

/ Et (/?) • fSirP' <en-ll [W] (V) + e ■ t [Dl] (V) + B ■ l" [Dl] {V) . 

We remark that it is crucial here that F = and hence F4 does not introduce 
3-derivatives (cf. Lemma 3.5). The other error-terms are dealt with analogously 
with no further difficulties involved. To give one more (interesting) example, we 
consider the worst part of the error-term 

j dt*drduj rP-^ El (p) • p4 + r-P^ EI (cr) • + r^^ eI (^) ^ , (364) 

which arises from the ^(p, (7)4 , -pair of estimates (cf. (336)). The leading 
order contribution from the first term is estimated (provided is not equal to 
6) 

j dt*drduj rP^Ha^pi < Bx j dedrdi^WHW^r^Wa^W^rP^"^ 
+A [ dt-'drdLo\pi\^r^P''-P^-^ < {Bxe + A) • li [W] {V) , 



with the last step following provided 2p^ — p4 — 3 < P2 — 1, which holds for 
a boundary admissible tuple. If P4 = 6, a similar computation leads to the 
condition {p3~ P2) +P3 < 8 — S, which is also valid. For the ^-term in (364), a 
typical term is estimated (say p^ 7^ 8) 



J dt*drduj rP^Hfii§^ j '^^*drduj\\Hfr'^\\^ 

+X j dt*drdoj\l3i\'^rP^''^ < [Bxt + A) • ^ [1^] {V) 



^^An even stronger improvement follows by a different argument. Namely, since we control 
the component |/34p in the main term of this estimate (331), we also control llZ^op with 
the same r- weight from the Bianchi equation (cf. Lemma 10.6). The Poincaro inequality 
Jg2 < /g2 then generates the improved decay of a. 
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requiring < pi, which is true for a boundary admissible tuple. For p^ = 8 we 
obtain 8 + 6 < pi- 

The estimate for (352) is straightforward. Again one can convince oneself 
that it suffices to understand j = and I = 0. In this case, for the first term 
(cf. Lemma 10.6) 

/ (E4{pf)rP^dvduj <\An^r'^\ [ \/3\^rP^-'^dvduj 

+ 1 {V,Z) f ^ , \a\\P^-Uvduj < en ■ E^^ [W] {N^ut {Sj^ «)) , (365) 

and similarly for the spacctime term. The other 7 terms arc dealt with similarly. 

□ 



11 A spacetime estimate for the Weyl-tensor 

Recall the energy (196). It turns out that we have control of this spacetime 
energy in the interior, provided we can control the p and the a component of 
the Weyl-tensor under consideration. To establish this result (Proposition 11.2) 
it is convenient to introduce the concept of spin-reduction. 

11.1 Spin reduction 

Consider the Maxwell pseudo-tensor 

Tap = r ■ WaP34 ■ (366) 

It is clearly antisymmetric. Its null components are 

TA3^2rP^ , TA4 = 2r(3A, (367) 

7^34 = 4rp , Tab = 2raeAB ■ (368) 

Note that in the Schwarzschild case, T corresponds to the electromagnetic field 
created by a point charge of strength proportional to M . 

As is readily computed using the formulae in the appendix, if V^Wabcd = 
J^bcd-, then T satisfies the inhomogcncous Maxwell's equations 

y'Tap = J-^'"' := rjpsi + rjp (369) 

with J7^34 a null-component of the right hand side of the Bianchi equation and 

J3 = P {trH - 203r) + a-H - 2Yj3 + 2ZP , 
j4 = p {trH - 204r) + a-H + 2Y§_ - 2Z/3 , 

JB = 1^ (trH - 2U)ir) + i/3 [trH - 2J])^r) ~ eAB*P_pHAF , 

-eAB*pFlLAF -Ya + Ya- p{Z -Z)+a{*Z + *Z) . (370) 

If the background in exactly Schwarzschild, the right hand side of (369) van- 
ishes, which is the origin of the name "spin reduction" . Note the "linear" (non- 
quadratically decaying) error-terms that arise in j in 73 and J4 that occur at the 
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lowest order (when p = p {W)) but that all terms are decaying quadratically in 
case that W = CipW, i > 1. 

The associated pseudo energy momentum tensor 

Qap = Ta-fF^-' + CT)^^ C^J-)/ = 2J-„^J-/ - ^go^pT^ (371) 
has null components 



Qab = (^SabP ■ P " - 13b §_^) + ^Sabt^ (p' + a') 



Q 



iA 



Br^ (p/3 - a* 13) 



(373) 



Note that Q is traceless, j^^Qab = Q34- Finally, we have the spinl-encrgy 
identity arising from a vectorfield X, 



(374) 



[J"] - Qa/J^^V"^ and [T] = (j}'^"'^ J"/ + {J^p"^^*^)^ (*J-)^^) 

These formulae should be compared with the spin2-identity (60). For the error 
one observes that 



J5z,T,' + {JlssA (*-^).') (^3)" = r^Kr^^^ m (375) 



:ild 



= ( + (^*),34 CJ"),') i^^r = r'Kr^'^ [W] (376) 



(377) 
(378) 



and recalls that the first two terms have been estimated before (however, without 
the additional weight) in section 7. The additional r^-weight does not alter 
the treatment of the intricate "linear" error-terms, in view of the strong decay in 
r of the component p ^ j^- However, the additional weight at infinity requires 
that an e of an r-weighted spacetime term needs to be added to Proposition 7.9. 
For this we recall (195) and state 

Proposition 11.1. The estimate of Proposition 7.9 holds for the expression 

(379) 

(with X = pxe3 +Qxei as in Proposition 7.9) as well, provided the spacetime- 
term e ■ Ip^ 



1 




IM 


\ J 0034: 



w 



M ) is added to the right hand side. 
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Proof. As mentioned, this is clear for the terms proportional to the non-decaying 
component p. It remains to look at the behavior at infinity. For the terms in 
j-spini g^j.jgj]^g fj-oni the original Jd34, one performs a null-decomposition as in 
[8] and inspects each term individually. For instance, 

/ dfdrduj rHrx^'^'> i-a- p{CtW) 

Jr>R 

< e [ dt*drdu}r"\a\^ + e [ dt*drdujr^\p {CtW) p < e • Tp^ 

Jr>R Jr>R " 

For the terms arising from j we have to look at the formulae (370). The esti- 
mates for these terms are straightforward and analogous to the above with one 
important exception: It is crucial that in our gauge y = 0, as otherwise the 
term r'^Y ■ a could give rise to an divergence). □ 



W 



11.2 The main result 

Proposition 11.2. Let W be a Weyl-field satisfying the inhomogeneous Bianchi 
eguations, T the associated pseudo Maxwell-tensor and recall (62). We have the 
integrated decay estimates 

sup Esp^nl m {tr) + E [W] {% (n, Ta)) + /^"^^ [W] (m (n, Ta)) < B ■ Esp„a [W] (s,,) 



+B r f dt*drr^ d J'' + lap] + B ■ ErrMora [W] 

Jro J ^ 

(380) 



sup £;,p,;„i [W] (tA + E [W] {H (ri, T2)) + r"^"'"' [W] [M (ri, T2) ) < 

(Ti,r2) ^ ^ ^ ^ 

B ■ Esp^ni m (Sri) + B j j dt*drr^ doj [\p\^ + \a\^] + B ■ ErrMora [W] 

(381) 

with 

ErrMora m= J2 If ^^2^'"'' 
a,b,c=3A ■'M(Tur2) 

Remark 11.3. In applications W will he the n-times Lie-T commuted Weyl- 
tensor, whose K2 [W] error-term was estimated in section 7, while the error- 
term arising from K2 [T] was estimated in Proposition 11.1. 

Proof of Proposition 11.2. We will apply the energy identity for the Bel- Robinson 
tensor (60) with the vectorfields 

X = 2X0 + 2C-r = /(r) (-3963 + 964) + 2C-r and 2r = peg -^(764 . (383) 

Here / is a bounded function behaving like /' ~ 7^-1-'' near infinity for a 
small constant > 0, while C = 2sup|/| is a constant chosen to make X 
timclike everywhere in TZ except for an e-small region close to the horizon. 



r'^K^, 



[F] (382) 
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For the following computations, we will use the "Schwarzschildean" normal 
2ris = \fk^&z H — 7r=^4:j which is close to ns, and then argue by stability of 



the estimates (cf. (64)). We compute 



2\a\ 



2„3 / 



-4 (p2 + a^) -p^qf. 



It is clear that by adding C • T to Xq we obtain the bounds 



pq Nk. 



Q{X,T,T,nf^) >b-Q{T,T,T,n'^ 



>{N,N,N,n'^) and 



(384) 



Q {X, T, T,n^)<B-Q (T, T, T, n^) + e • Q {N, N, N, n'^) . 



For the bulk term Ki [W] wc first compute the components of the deforma- 
tion in case the metric is Schwarzschild 



(2^o)_34 
^SS 



4 (1 - m) - ^/ 



(2Xo) 33 _ 1 r 



(2Xo)^44 ^ ( f (2Xo)AB / ,,\xAB 

ss=2y—j J,r , 't^ss ^ - - P)o 

and then decompose 

K^^^ [W] = A'f/^ [W] + A'i^^^ [W] + K^J,:^,, [W] , 

where 



1 



\ r 2 r 



k+ 

X 



(385) 



(386) 



(387) 



m = (1 - M) / (' - 3^) [p' + - ^ (1 - A*)' [P' + 

are the expressions one would obtain if the components of the deformation tensor 
of Xq took on exactly their Schwarzschildean values (and if T was replaced by 

k+ 

the Schwarzschildean T = -^k^e^ + 2^^^) arid 

K^Jrlor m = (^^'- - Vss) Q [T, T) + 2 ((^).) Q (X, T) 

+2 (^^'tt) Q (^T - f , + (^(-^'tt) Q (t -f,T -Tj , (388) 
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is the error introduced. One easily sees that there is no way to make Kia and 
Kif, globahy positive at the same time, in view of the components p and a 
having the opposite sign compared with the other curvature components. 

To avoid degeneration of control near the horizon we will add to the energy 
identity of XTT a little bit of NNN (we drop the [Wj-argument for the moment 
to simplify the notation) producing 



M 



XTT 
b 



+ / J, 



XTTu 



At 



XTT 

error 



K. 



XTT 



eJ. 



NNNu 



NNN 



H 



J, 



jXTT u 



+ eJ, 



NNN^fj. 



(389) 



for e > a small constant. We choose e so small that in particular 



1 



1 



ry - 2M 



holds. Observe that K^^^^ has a positive sign close to the horizon and near 
infinity so that we only have to estimate this term between two fixed constant 
r hypersurfaces. Note also that the boundary terms in the second line of (389) 
have a sign and control the non-degenerate energy. 

To obtain the estimate (380), we would like to choose / (3M) = and f,r to 
be a smooth function which is always non-negative and vanishes quadratically 
at r = iM . This almost gives estimate (380), except that the /3-tcrm now 
degenerates like (r — 3M)''. Moreover, the same problem occurs when we try to 
prove (381). In this case we would like to choose / to change sign at 3M and /.^ 
everywhere positive. Clearly, the /3-tcrni would still degenerate quadratically 
at 3M. We resolve this problem using the X-cstimate for the Maxwell-pseudo 
tensor: Applying (374) with the same vectorficld X , wc find 

/ [F] + + I [F] n^^ + [ [F] ,4 = f [F] . 

JM -IT. JH JT.Q 

(390) 

For the boundary term J^nf^^ — Q {X, n„j,) we have the analogue of (384): 

b-Q{T,nnJ-e-Q{N,nn^) < Q (X, n„ J < B • Q (T, n„ J + e • Q (^, n„ J . 

(391) 

As in the spin2 case, we decompose the bulk term as 



(392) 



where now 



m\ 



(393) 



16 (p' 
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and 

Ktrror [-^] = (^""^^ " ^""^^Ss) 'Q + C- ^^V • Q . (394) 

We remark that adding to X a little bit of the vectorfield, one can eliminate 
the degenerating weights close to the horizon in the standard fashion (it is 
unnecessary at this point though, as we already dealt with this problem in the 
spin2 estimate). 

Finally, to obtain the estimate (380), we choose / (3A/) = and / to be a 
smooth function which is always positive and vanishes quadratically at 3-A/. In 
this case the desired spacetime integral is controlled by 

Id^9 [w] (n, r^)) KBjdt* drr^ J dJ-^—^^ [\p\' + \a\'] 

+B f [K^r m + K^,^^ m + ex^^ m + ki m + k^, . (395) 

JM 

The right hand side is now estimated by inserting the identities (389) and (390): 
The resulting boundary terms have a good sign on the horizon and are otherwise 
estimated by (384) and (391). For the error-term we have 



JM 



^^1, error ^ / „ -^1, error [^^] 

JM(ti .To) 



< 



I M(ti,T2) J M(ti, 

e • /"^^^ [W] (m (ti, T2)) + e • sup E,p„,i [W] (tr) , (396) 

as an consequence of the ultimately Schwarzschildean assumption. Hence the 
estimate (380) follows after absorbing the terms on the right hand side of (396) 
by the main terms on the left. 

Similarly, to obtain (381), we choose / to vanish linearly at 3il/ and f^r 
everywhere positive. Note that since f^r no longer degenerates at 3Af. we now 
need non-degenerate control of p and a on the right hand side. The estimate 
(395) now holds without the degenerating weight on the right hand side and the 
non-degenerate energy on the left. The second line is estimated as previously, 
using the energy identities (389) and (390). □ 

11.3 The estimate at the lowest order 

For W ~ W the original Weyl-field, the integrated decay estimate of the previous 
section is not useful, as the right hand side is expected to grow in t in view of the 
component p appearing. However, using multipliers directly on the renormalized 
nuU-Bianchi equations, just as we did in section 6.4, we obtain 

Corollary 11.4. For the original Weyl field W we have the following estimate 
for the renormalized components: For any h > 

sup / \\w\\'+[ \\w\\'+f -^\\wr<BH[ \\wr 

{Ti,T2)Jt^ J-H(Ti,r2) JM{tuT2)'^ "'S^j 

+B j dt* r drr^ f doj[\p\^ + \a\^]+Bh\\p^\\L^ f \\D\-D\ss\\^ 

J Jra J J M{ti,T2) 
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12 The components p and a 



As we have seen in the previous section, the components p and a play a special 
role in the decay. In this section, we derive the wave equations for p and a and 
their T-commuted higher order analogues, define appropriate energies for them 
and prove multiplier estimates for the latter. 

As mentioned previously, the natural variables to look at are the renormal- 
ized quantities 



^3 



and = r^<T (C^W) , (397) 



for which we defined the energies E\V(f>n] [^-rj ^^.d Ideg [D4'n\ (^) (and sim- 
ilarly with ■0„) in (198) and (199). Note in this context that the relation 
E \D(j)n] {pT^ < B ■ Ep^^ \W] i^T^ holds for the boundary admissible decay 
matrix defined in (195). 

The quantities (397) will be seen to satisfy an inhomogeneous Regge- Wheeler- 
type equation (section 12.2). In section 12.4, we derive an integrated decay 
estimate for the homogeneous part of this equation. Such type of estimates 
are known [15, 16], in the context of the Regge- Wheeler equation on an exact 
Schwarzschild background. The crucial point here is that the lower order terms, 
which are typically the source of difficulty in deriving such an estimate, are a- 
priori controlled from the ultimately Schwarzschildean assumption. The second 
crucial ingredient consists in controlling the inhomogeneity which is present in 
our case. This is particularly intricate in view of the additional r'^-weight intro- 
duced by the renormalization (397). The analysis is the content of section 12.5. 
Finally, we remark that all error-estimates are carried out for the p-equation 
only, since the a-analogues are easily inferred (section 12.3). 

12.1 The main result 

We start by stating the main result of this section, an integrated decay estimate 
for the quantities (/>„ and ipn , which follows from combining the decay estimates 
for the Regge- Wheeler operator with that for the inhomogeneity. 

Proposition 12.1. Recall the boundary admissible decay matrix Pp (cf. (195)) 
and its bulk associated matrix Pp, as well as (397). For 1 < n < k — 1 



sup E [D<j>n] (f^r) + E [V^n] {U (ti, Ta)) + Ideg ['D(j)n] (m (ti, T2) 
(ri,r2) ^ ^ ^ 

<BE [Dcjjn] (Sn) + e Tp^''''"' [W] (n, T2)) + supE"+' [W] (s.) 

+A • ©"+1 [m] (ri ,T2) + Bx- ©" [m] (ri , T2 

• Ideg [2?0n-l] (M {tu T2)) + B-E [Vc^n-l] 



for any A > 0. For n = 0, the same estimate holds with the last two lines 
replaced by the term B ■ [SH] . Moreover, rhe same estimates hold replacing 0„ 
by i>n- 
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Proof. Proposition 12.1 will follow from applying Proposition 12.13 and esti- 
mating the error-term via Propositions 12.14 and 12.20. See section 12.6. □ 

Remark 12.2. Note that the energy I^i^g [Vcjjn] provides at the same time a 
non-degenerate integrated decay estimate for (t>n~i- 

The above result should be understood as follows: It is possible to control a 
spacetime integral for the T-commuted renormalized p and a components from 
the L^-energy at the expense of 

• borrowing an e from the spacetime integral of a weighted energy (decay 
matrix Pp) of all the other curvature components (second line) 

• borrowing a small amount (A) of the highest energy on the Ricci-coefhcients 
(third line) 

• and, provided the analogous estimate is available one order lower (last 
line). 

Remark 12.3. It is crucial to note that the analogous degenerate estimate for 
all curvature components requires the full D" [$H] instead of D" [9^] . The point 
here is that in proving a non- degenerate estimate for all curvature components 
one only loses derivatives for the quantities p and a , i.e. one only needs a higher 
order estimate for p and a, not for the other curvature components (cf. Proposi- 
tion 11.2). For these scalar equations, however, one can exploit the fact that one 
gains non-degenerate integrated decay for a particular radial derivative, which 
allows the gain of a small A when controlling the error-terms. 

Clearly, one can eliminate the last line of the estimate of Proposition 12.1 
by iterating the estimate: 

Corollary 12.4. Forl<n<k-l 

sup E [V(j)n] (Sr) + E [2?0„] {n (ri, T2)) + Ideg pC^n] (M (n, T2)) 
(Ti,r2) 

n 
i=0 

+A • D"+i [9^] (n, r2) + Bx ■ B" M (n, rs) 

for any A > 0. For n = the same estimate holds with the last line replaced by 
the term B ■ The same estimates hold replacing 0„ hy ipn- 

12.2 The scalar wave equation for p 

We want to derive the equation for the p and a components of the commuted 
Bianchi equation 

D-^ (C^W) = Z],^s : (398) 
where for the inhomogeneous term we have by (75) 

n-2 

Tp^s = C^^'Jp^s (T, W) + J2 ^tJm& {t, IVf^-^ZrW) . (399) 

i=0 



I^;^^'^^ {W\ (m (n, r^)) + supr+^ \W\ (S.) 
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Recall the notation for the null-decomposition of the inhomogeneity introduced 
in Proposition 7.1, 4A [-3"] = 3434 etc. and also the notation for the inhomo- 

geneities in the null-decomposed Bianchi equations, £'4 (p (C^W) ), cf. (91). 



Proposition 12.5. The quantity p yClj^W j for ?i > satisfies 

2 



+p (^C^W^ (itrHtrH + Up (W) + {dec. RRC)^^ ^^^i + S ^- > 
where the right hand side has the following structure 



-H ■ YI3 [Cr^W^j + trFs i^P [C^Wjj - 2nEi (^p 
+2y^fl^ (C'^w) + div (2H ■ p (C^w) + y • a (c^W 



+'ic'ilrlZa {pi^W^ + iZ-^^AP (^t^") + 3 (^*)^ ^ (^JW^) (400) 

F!^ = -2 (A [3"])3 + m{k [r]) - Idiv {I [T]) + 3po • p (C^w) (401) 
and F^l being the conjugate quantities. 

Remark 12.6. The expression on the left hand side is essentially (see below) 
the Regge Wheeler operator of the metric g acting on p [C^,W^ . If the right 
hand side exhibited only quadratic error-terms, proving decay would be straight- 
forward. However, this is only manifestly true for and Ffl, while for T2 
and there are in fact terms proportional to the non-decaying component p. 
The structure of the error-terms will be discussed in section 12.5. 

The last term in (401) will be seen to cancel by contributions from the other 
terms in case n > 1. For n = 0, however, 5" = and + ^ = 3p^ + 3p^ and 
hence 

Corollary 12.7. For n = we obtain 

{lt)^Jt)^ 0304) P - '^'f^P + \ {^trH ~ Ij)^p 
+ \ (^2trH - 04P + p (strHtrH + 6p + {dec. RRCf^ = + Z? • 
Note that the right hand side decays quadratically.^^ 



^®In particular, the equation for p itself decouples modulo quadratically decaying error- 
terms. However, p does not satisfy a Regge- Wheeler equation and only the renormalized 
quantity p -|- can be seen to satisfy an equation with the Regge- Wheeler operator on the 
right hand side. Unfortunately, this renormalization also introduces "linear" error-terms of 
the form "p X "D (IH — 5Hg5)", in view of the fact that p itself does not decay. A more careful 
choice of the renormalizing function may remedy this problem, cf. the remarks at the end of 
this section. 



97 



For the following proof we will, to avoid overloading the notation, denote by 
a, a, p, a the null-components of the n-commuted Weyl tensor CJ^W, which 
satisfies the n times T-commutcd inhomogeneous Bianchi equations. Moreover, 
we denote by po the (non-decaying) p-componcnt of the original Weyl curvature 
field, i.e. po = p (W). 



Proof of Proposition 12.5. Start from the scalar Bianchi equation for p = p yC^W 

Pi = Ipip + '^tr [H) p = divP + Ei (p) - ^^434 . (402) 
Taking a 3-derivative we obtain 

P43 - = div {Es (/?) + ^3^4) + ^3^4 (p) + ^trHEi (p) 

-P ■ 1/trH - 1 • y/3 + trFg (/3) - ^03^434 - \tr (H) 5434 , (403) 

where 

Es (13) = 2H ■ P + 2ni3 + Ya + 3{Zp + *Za) . (404) 
We first deal with the left hand side. Using the null structure equation 

IPstrH = -^trHtrH - 2divZ + 2ntrH -H -H 

+2{Y-Y + Z-Z) + 2po (405) 

with Po denoting the p component of the original Weyl-tensor, we derive 

3 3 

LHS = P43 - f% = IPzIPiP + 2 it^H) $^p + - {trH) l^iP + f^P 
+ ^P (trHtrH + 2divZ + 2ntrH - H ■ H + 2 {Y- Y + Z ■ Z) + 2po) • (406) 

Turning to the right hand side of (403) we decompose 

1 13 
RHSi = - - (^434)3 + div {Z3A4) = -2-030434 - -^tr {K)Z43A + div (^3^4) , 

RHS2 = IP^Ei (p) + divEs (/3) + ^trHEi {p) - P ■ ftrH ~H-fp + trF^ (/?) . 
Using again the Bianchi equations we compute 

RHS2 = IP^Ei (p) + ^trHEi {p)~p- 1/trH -H-i^p + trF^ (/3) 

2n (^]pip + hrH p~Ei (p) + ia434^ + 2f^n/3A 

+div (2H ■P + Y-a^+ idivZp + iairlZa + 3^^^^/? + 3 {Z*)^ ^^cr (407) 

and observe that all terms containing p without a derivative are cancelled by 
terms in LHS. The term Q^lpip is moved to the left hand side. Repeat the 
same computation starting from p^ and applying a 4-derivativc. This yields the 
conjugate terms. □ 
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To obtain the energy naturally associated with the equation, we will have to 
rescale the equations of Lemma 12.5 by r^. Recall that (/)„ = r^p(^CJ^W^ for 
n > 1. Define he vectorfield 

L = -i(2-p)e3 + i(?e4, (408) 

which equals dr in Schwarzschild with respect to the regular (t* , r)-coordinates. 
RecaU 2T = pes + qe^ and 2T-^ = -pes + qa and hence {I - p)T + pL = . 

Corollary 12.8. For n > 1, we can write 

(2 - p) (0„) - L {pLcj)„) + {L (p)) T {(bn) - 2 (1 - p) Ti0„ 
-qf^(t)n ~ g ^^^3^^ 0« = I ^^-7^" + (^decaying RRC , ^^V) • I?0„ (409) 

where \E\ < e is a small error (C"^ in the metric, independent of n) arising from 
the fact that the metric is close to Schwarzschild, and the last term denotes 
a collection of terms of the form "any derivative of multiplying any null- 
component of the deformation tensor of T , or an expression d\ — " 
The renormalized quantity (po also satisfies (409), except that terms of the form 
Pq ■ (^decaying RRC , (^■'tt) need to be added to the right hand side. 

Proof. For n > 1, we first write the equation of Proposition 12.5 as 

(0403 + 0304) 0„ - 2y - 21]030„ - 2il040„ + 0„ {6p (W)) - J"" , 
which holds modulo the error specified in the Corollary. Secondly, we write 

(0403 + 0304) - 2r!03(/)„ - 2il04(/.„ = 

- [(2 - p) {4>,,) - L {pL4>^) + {L ip)) T (0„) - 2 (1 - p) TL0„] , (410) 
q 

which again holds modulo the error specified. For n = 0, we start from Corollary 
12.7 directly and do the same computation. □ 

Clearly, the exact form of the last term on the right hand side of (409) will 
not matter, since any such term will be easily estimated in view of the strong 
pointwise decay available for the expression in brackets (cf. Definition 2.2). 

We finally remark that, in case that the metric is exactly Schwarzschild and 
Edington-Finkelstein coordinates are used, the operator on the left of (409) 
reads 

-5„a^0n + ^-^A520„ + (l-/i)^(/.„, (411) 

which is the Regge- Wheeler operator in its perhaps most familiar form. 
12.3 The scalar wave equation for a 

The computations for a arc similar. In particular, we have the analogue of 
Proposition 12.5 (this can either be computed directly or inferred from duality): 
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Proposition 12.9. The quantity a \ CVj,W^ satisfies, for n>Q 

(0403 + 0304) ^ {CtW) - 2f^a 
+ ^ (^2trH - lp,a (^W^) + ^ (^2trH - Ip^a (^C^w) 

2 

+(T (strHtrH + 12p (W) + {dec. RRC)^^ = ^ + , 

where the right hand side has the following structure 

Gl' = 035^4 (cr (^tW")) + ^trHEi (a (^t^^)) + *P • ftrH 

+H ■ f*l3 (fjw) - trFs (*/3 (^JW^)) - 2il£;4 (a (^t^^)) 
+2^*11/3 (rjH^) - c^rl (2H ■ §_(c^W^ +^-a (^JVk) ) 

^icdrlZp {C'^W^ + SZ^cr - i* Zf p [c'i^W^ , (412) 

G^ = -2 {K [rDa + 4il (K [r]) + 2c4rl {I [5"]) + 3po • ^ (^t^^) , (413) 
and 02 being the conjugate quantities. 

For the n = case, one can use the nuU-structure equations ci/rlZ = —^H A 
i£ + A y + cr plus its conjugate to conclude 

Corollary 12.10. For n ^ 0, one obtains 
{ip^lPs + lp-:,lpi) a - 2f\ + ^ (^trH - Ul^ lj)^a + ^ (2trH - 04a 
+cr (strHtrH + 12p (W) + (dec. RRCf^ =Gi+G^i+HAH + 2YAY. 
Note that the right hand side decays quadratically. 

12.4 The integrated decay estimate 

We apply general multiplier identities to equation (409) to derive energy esti- 
mates. As previously mentioned, decay for the (homogeneous) Regge- Wheeler 
equation on Schwarzschild has been shown ([15, 16]). In this section we present 
an alternative multiplier approach based on ideas of [1]. The main result is 
Proposition 12.13 which together with the error-estimate of Proposition 12.14 
will finally imply Proposition 12.1. 
Multiplying (409) by the multiplier 

X =[a{r)T+b{r)L + c{r)](l)k (414) 

with o, b,c being bounded functions, one derives (let us momentarily drop 
the subscript n) 
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Lemma 12.11. We have the following basic multiplier identity 



+ [V(j)] ^ Err [V<j)] + [r^ J"^'] (X^k) 



where 



X 



[Vcjy] = Vt 



ia (2 - p) |r0p + Qap - b (1 - p) ) 



If 6M + £ 



aq + c L{p)]^^ + b{2-p)T{^)L{^) 



+c(2-p) T(0).0-2c(l-p)L(0)-0 



-a (1 - p) - ib (2 - p)^ |r0|2 - ibp|L0p + ibgiy^p 
p • L (c) - bg^^^^ ) 02 _ (0) r (0) _ PC 



ogy./) • T<f) - bqfcf) ■ L<j)~ cq (fcj)) 



(415) 



(416) 



[V<j)] = 



{l-p)L{a) + -L{b{2-p))-c{2-p) 



^pL{b)-^bL{p)+pc 



-^L{bq) + — + cq 
2 r 



-T (0) • L (0) [pL (a) + bL (p) + 2c (1 - p)] 



1 / 6M + S 

2^r'^^^ 



-L {pL it)) -qc ^3 



and i/ie error-term is cubically decaying, schematically of the form 

Err [Dcj)] = (^decaying RRC, ("^^tt^ • V(j) ■ {V(j) or </)) . (417) 

Multiplying this identity by and then integrating over M (ri,T2) with 
respect to the measure dt^drdOdcf), we realize that is a boundary term 

modulo terms of the form Err [Dc/)] . 

With the Lemma at hand, we can prove an integrated decay estimate for the 
Regge- Wheeler-type equation under consideration. In view of the trapping, the 
best one can hope for is to generate a spacetime-integral term which provides 
non-degenerate decay for a particuar radial derivative and degenerate (at the 
photon sphere) integrated decay for the remaining derivatives of (j). In the 
context of the wave equation, obtaining positivity for the zeroth order terms in 
(j) is intricate and typically involves the use of a Poincare inequality exploiting 
the positivity of the angular-derivative term. In our case, the situation is more 
favorable, as we are assuming convergence to Schwarzschild. In particular, the 
assumption of being ultimately Schwarzschildean to order A: -I- 1 implies that the 
integrated; non-degenerate spacetime energy of A: — 1 derivatives of curvature is 
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bounded. This means that in the putative degenerate integrated decay estimate 
for derivatives of ipk-i, the lowest order terms (pk-i are a-priori controlled from 
the ultimately Schwarzschildean assumption. Moreover, we can exploit the fact 
that the degenerate integrated decay estimate for derivatives of (j^k-i will provide 
control over a particular radial derivative of fpk-i near the photon sphere (cf. the 
energy (199)). Indeed, this allows to turn non-dcgcncratc decay into degenerate 
decay plus a small amount of the non-degenerate derivative term, as is spelled 
out in the following Lemma: 

Lemma 12.12. For any A3 > and n > 1 

/ dt*drdu2^TT<l)l < A3 • Ideg [T>(j)n] (M (n, T2) 

+ ^-Ideg[D(bn-l] [M {n , T2)^ + € SUp E [Vcj),,] (s,) (418) 

Proof. Decompose 

= + c^r^y = x{r)- + (1 - X ir)) cj^k (419) 

where x ('') is a function which is equal to 1 in [3A/ — ^M, 3M + and equal 
to zero outside of [3M - ^M, 3M + \M] . For 0^'"''^, the estimate (418) holds 
trivially requiring only the second term on the right hand side. For 0^''°*°" one 
has the estimate (recall that T-^ = ^9** + (1 — ^) dr in Schwarzschild) 

X'<P' < 6x'<P' {T^ {r - 3M)) 
= er-L (r - ZM)) - 12x^ (r - 3A/) {T^(j)) - l2xT^ (x) {r - MI) 0^ 

< 6T^ {x^<f' {r - 3M)) + ^ (r - 3Mf <f>^ + A3 (T^0)' , 

A3 

since T-^ (x) = in [3 A/ - |A/, 3 A/ + |A/] . Integrating the estimate with re- 
spect to J j^dt* drdoj gives the estimate of the Lemma. The contribution from 

the boundary term on S,- is e-small in view of the ultimately Schwarzschildean 
assumption (g (T-^,dt*) ^ e). □ 

Proposition 12.13. We have for n > 1 

E [P(^„] (S,,) + E [2?<^„] [n (ri, T2)) + heg [2?<^„] {m (n, T2)) < 

B-E[V(bn] (f^r,) +B 

B ■ Ideg [Vcj)n-i] [m (ti,T2)) +B-E [2?(/)„_i] [^r,n+j . (420) 

Proof. Consider the following multiplier: 

X = B1X1+ &2^2 + X3 + 64X4 (421) 

for Bi a large and 62,64 small constants depending only on A/ (to be chosen) 
and 

Xi=T , X2 = -7(r)L , X3 = il-p)fir)T + fir)pL + ^pL{fir)) 



dt*drdujr^T'' (X(/)„) 
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where we recall the redshift multiplier 7 (r) defined in (213) and 



Finally, 



(r - 3My ~ 



(422) 



(423) 



where ^ (r) is a cut-off function supported for r > ry 



> 2M and equal 



to 1 for r > ry- From Lemma 12.11 we compute, the integration being over 

^(Tl,T2)), 



J K^^ [V(t)] = J dt*drduj 



L{f){{l~p)T{<i>)+pL{<i>)f 



\ i^f ■ L ipL (pL if))) - \fr^L (^) ly^p + L 



MI + £ 
— ^3 — P^i 



Since / is monotonically increasing, the first term is already non-negative. More- 
over, as / changes sign at 3M, the angular term is seen to be non-negative, too 
(modulo a term of the form Err [Dcf)]] in Schwarzschild L(ff) = -^(l--^)). 
Finally, the lower order terms in (p itself are treated as error-terms and will 
eventually be controlled with an application of Lemma 12.12 and the ultimately 
Schwarzschildean assumption. Hence up to now we have shown the estimate 



B I ^dt*drduj 



K^"- [Vcj)] 



> b j dtdr*doj 



-4 ((1 ~P)T ic^) + pL + ~ 3^^' 



which controls all derivatives except the 9^* -derivative. The latter is retrieved 
by adding a little bit of X^. It's not hard to see that for our choice of (note 
that X4 is supported away from the horizon) we have 



B 



[ ^dt*drdijj + ( K^^ [2?0] +B [ K^" [D(t)] 



> b 



IM(ti,T2) ^ 

dt^drdw 



1 



2 (r - 3A/)^ 



' Mn{r>rY} 

Turning to the redshift vectorfield X2, we observe 



[ K^' [Vet)] +B [ 



>b dt*drduj 

' Mn{r<rY} 



dt*drdu) 



M>r>rY} 



101' + im' 
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Consequently, we can find small constants 62, &4 such that (reinstalling the 
subindex n) 



-B [ ^dt*drdij, (424) 



/A^(ri 



where /^"^ [P^n] (ri,T2)j denotes the energy (199) with (5 = 1. The last 
term is estimated by Lemma 12.12. 

Before we optimize the weight at infinity to the S = we specified, let us 
note that by choosing Bi large enough, we have for the boundary term 



I 



e 



^ r-^ ,] > 5 • S [V^n] [tr,'H+]-B-E [l?0„-i] Sr , ) . (425) 



This follows because X is timelike for sufficiently large Bi (noting that both 
X3 and Xi degenerate at the horizon) and that moreover, we can estimate the 
zeroth order term by the lower order energy. Finally, 

/ Err[V(j)n] <esupE[V(t)„](tr) . (426) 

Jm(ti,T2) r V / 

In summary, we have established Proposition 12.13, except that the spacetime- 
energy docs not have the correct r- weights at infinity yet. As this optimization 
is standard, it is omitted. □ 

12.5 Error-terms: The structure of the inhomogeneity 

It is apparent from Proposition 12.13 that we finally need to understand the 
error-term arising from the inhomogeneity. There are two important regions: 
The interior region, where we need to have sufficient decay in t and be careful 
about the degeneration of the estimates near r = 3M, and the exterior where 
we do not encounter the problem of degeneration but, on the other hand, need 
to worry about the weights in r for the error-terms that arise, especially in view 
of the fact that the renormalizcd quantities have introduced a large r'^-weight 
for these terms. 

Proposition 12.14. Let n > 1 and X he the multiplier used in (421). For any 
A > we have 

2 I- 

J2 / dt*drduJ2r^ ( J-f ) {X(j)n) 

< A • Ideg [V<j)n] (m (n, T2)) + A • sup S [!?(/)„] (s,) + A • £; [Dcj,^] {%) 
It;^'-' [W] [m (n, r2)) + supr+^ [W] (S.) 
+ supE"+' [$n] (s,) +r [5n] {m (Ti,r2) 

-B 



A • Ideg [VK-i] {m (ti, r2)) + A • [% {t^^t^) + Bx ■ O" [m] (n, r2) 



'Here one can simply redo the argument with an / in (422) which is bounded, monotonically 
increasing, changing sign at the photon sphere and decaying hke /,r ~ near infinity. 
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Remark 12.15. Note that the terms in the first line will he absorbed by the left 
hand side of (4-20) for sufficiently small A. 



Proof. We begin with the term J^" . It is obviously an expression which decays 
at least quadratically and hence improves upon taking derivatives. For r < R 
we can exploit the pointwise decay of the Ricci coefficients (and first derivatives 
thereof): 

dt*drduj2r^ (J"") iX(j>) 

M(Ti,T2),r<R 



< e • sup Je (I?0„,I],,,<h) • sup y/e"+' [W] (I],.,<i?) . (427) 

T ' T 

For r > R we can borrow a bit from the weighted energy for the Weyl field 
and estimate the X(/)„-term in a spacetimc integral, as we are far away from the 
trapped region: 

/ dt*drdLU2r^ (JT) (^0n) 

JM{Ti,T2},r>B. 

< Ai • Ideg [V<Pn] (m (ri,r2)) + e • Sa, ■Tp^^''"'' [W] (m (ti,T2)) (428) 

with the e arising again from the pointwise bounds on the Ricci-coefficients. 
The same estimate holds for the conjugate quantity 

Turning to J-2 we have to understand the structure of A = J434 and / = 
J3A4 and Lie-T-derivatives thereof, in particular the asymptotic decay in r on 
the one hand and the terms not decaying quadratically on the other. Indeed, it 
suffices to understand A and / as arising from one commutation in T, since the 
null-decomposition almost commutes with the T-derivative: 

n-2 

A [a] = [A (J (T, W))] + ^tA ( J (r, £J-i-*£tVk) ) + 3^4°' , (429) 

4=0 

where the last two terms on the right hand side decay quadratically. We first 
split 

-^2' =-^2" +-^2!., (430) 

where J-2a collects the terms in which at least two factors decay and J^jh ^^'C the 
potentially dangerous "linear" terms, i.e. those proportional to po- 
we decompose 



For 



•'2a 



• {Xcf) = r'T^^ ■ [X^) • x (r) + ■ {X<i>) • (1 - x (0) = 2( + 'B (431) 

where x (r) is an interpolating function which is equal to 1 for r <R~ M and 
for r > R. Note that 2t is supported for r < i? only, while *B is supported in 
r > R — M. For *B we have 



M{ti.T2) 



dt*drdtj r^Tl, ■ {XK) ■X{r)<\x- heg [m (n, ts)) 



1 

'^l J M{Ti,T2)n{r>R-M} 



r'^'x' i^^ar ■ (432) 



For the second term wc show 
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Lemma 12.16. 



JM(Ti,T2)n{r>R-M} ^ ' 

Proof. To estimate these quadratic error-terms in the far away region, we only 
have to be careful about the weights in r. To reveal the structure which is 
present for the error-terms (cancellation of worst terms, cf. section 10.4), we are 
going to think about the wave equation for p in yet another way. At the zeroth 
order we can write 

P34 + P43 - 2f^p = -divEi {§) + El (p) + dj[vE3 -I- El (p) . (433) 

We now imagine commuting this equation with the Lie-T-derivative directly to 
obtain the higher derivative version of the equation. Note that this commutation 
will not produce any additional difficult terms at infinity, but simply Lic-T- 
derivatives of the right hand side of (433).^^ However, the decay in r of the 
right hand side is not altered by application of T. It follows that it suffices to 
understand the decay in r of the right hand side of (433), more precisely, to 
establish that it decays better than r~'^~^ pointwise. We note in passing that 
the gauge condition y = is not needed here and that we will hence carry out 
the estimates including this term. 

g 

Let = denote equality up to terms which decay like r~2 or stronger. We 
observe that 

~divE4 {§) + divEs {(3) = div {Ya) (434) 

and, noting that in view of the null structure equation (174) Ip^X actually 
decays like (instead of the naive ^) also 

Et {p) + El {p) = Tt)lY-a + Y^$^p_ - 2YD3^ . (435) 

Adding up these terms one sees that the worst terms cancel by the Bianchi 
equation (88), producing terms decaying like r~^. □ 

We emphasize that it is absolutely crucial to exploit the cancellation of the 
worst error-terms as otherwise a r~'^-divergence would arise. 

For the term 2t in (431) we do not have to worry about r-wcights at all, only 
about the order of differentiability in view of the degeneration at the photon 
sphere. However, all the terms in the expression for J^2a decay quadratically. 
This makes the error-term cubic and the estimates are straightforward, putting 
the highest terms in and the lowest in L°°. Using Sobolev embedding, we 

^®It may produce terms proportional to p from commuting T past the derivatives on the 
left hand side. However, these decay strongly in r. 
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obtain 



M{Ti,T2).r<R 



{X(j)n) < {t2 - Tl) sup JE [P0„] 



i=0 



supWE'^ [^H] (SJ supWE"^ '[W]{tr 



n+l 



X! sup\/E 



■n+3— -i I 



[w] sup We [$h] 



< 



A • E [2?<^„] (e,) + (t2 - Ti)" supeL"^'J+' [9^] (s,) • supE"+' [W^] (e,) 



-B(r2 -Ti)^supE 



■ri+l 



VK] Er • supE [Sn] E^ 



and in view of the decay assumptions on the lower oder energies (recall the 
ultimately Schwarzschildean assumption and n < k — 1), the t- weights can be 
absorbed, producing the terms on the right hand side of the proposition. 

, i.e. the terms proportional to po- These terms 



It remains to estimate 



2b 



are particularly difficult near the photon sphere r = 3M (in fact, near infinity 
they are harmless, in view of the strong decay of po)- We first note that it 
suffices to understand the case = 1, since higher k are just T-Lie-derivatives 
of this expression (and commutation by Lt introduces additional decay factors). 
Hence wc want to compute the expression 



2b 



^4^0 + ItrxAo - 2QAa + dkh, - ip [W) ■ p (CtW) 



where 



Ao 



^04 [tru) - f^VA + 2pp - UrH {2f^p - 4flp) 



2qPA - ^3iZA - 2nrA - trxjT^ - 2trxTT\ + ^-^Va + 1^1^ a [tri:) 

(and similarly for the conjugate^^ quantities Aq, /g) were computed at the end 
of the proof of Lemma 7.4. 

In fact, we will only need to do evaluate J-'jf, modulo quadratic error-terms 
(as the latter immediately transfer to J'^a)- Wc make the a-priori convention 
that for the computations within the context of the errorterm 7^265 
means "equality up to quadratically decaying terms". 

^^The conjugate quantity of /9 is — /3, while all other quantities just receive an underline and 
p becomes q. 
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Lemma 12.17. // the vectorfield T is chosen such that its deformation tensor 
^"^•'tt is traceless (which it is according to Definition 2.1), then (n>l) 

r^^, = po ■ VC^-^ (tt, dec. RRC) + l.o.t. (436) 

Proof. Direct computation, using that 

2^^3/9 + 2g^4p - 2q1/^p = 2T (p) - 2qf^p 

= 2p (^CtW^ + ^trnp + 3qtrxp - 2qdivf3 (437) 



and 



3 r 

divio ~ —po 2qdiv(3A — i^ndivZA — 2QdivVA ~ trxdivir^^ 

1 9 

-2trxdiv7r^A + div^^^VA + -^f" (trn) 



□ 



Remark 12.18. Note that the Lemma once more exploits a cancellation of the 
worst terms: Naively, (4-36) only holds with n instead ofn—1. 

Remark 12.19. // the deformation tensor is not traceless, J^2b contains an 



trir. 



additional higher order term of the form PqC^ ^ (^Ip^Ip^ + Ipi^lp^ — 
However, the operator in brackets is the principal part of the Regge- Wheeler- 
operator. Hence via several integration by parts, moving the operators to 
in the spacetime- error-term arising from J^2b' ^^^^ term can also be seen to be 
of lower order. It is here where the technical assumption tr^"^^ = ^ ( discussed 
below Definition 2.1) helps us to reduce the number of error-terms. 

With the Lemma at hand, we pick a cut-off function x* (^)i which is equal 
to 1 for r < lOM and equal to for r > 11 A/ and decompose 

-^21, = ^2b ■ X* (r) + [1 - X. {r)] ■ (438) 
We can then estimate away from the photon sphere 

dt*drdw2r^J''2l [1 - X* (r)] ■ (Xcj^n) < A • Ideg [V<j)n] (m (Ti,r2)) + ^D" [m] , 

M \ / A 



while for the other term we integrate by parts, writing a shorthand / to denote 
^i,»ri dt*drdLu and noting that weights in r are bounded in the 
region under consideration 



26 ■ X* (r) ■ iN<j)n) - J r^pn {Cl^r^V {n, dec. RRC)) • (X^n) 

< +\2- sup E[V(j)n] (tr,n) +Bx,_ • SUp E" [$H] (Sr,^ 



X2 J (l?2/:j"^ (tt, dec. RRC))%^ J <f)l + B J (P^J^^ (tt, dec. RRC))^ 
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for any (small) number A2 > 0. Applying Lemma 12.12 with a A3 much smaller 
than A2 to the first term in the last line finally yields 

J r^J^b ■ X* ir) ■ iX<Pn) < A • heg P(t>n] [M (ti, T2)) 

+ (A + e) • sup E [V^n] (S,, n)+B^- I^eg [m (n, T2)) + B^'+i [m] , 

for any small A. Clearly, the first two terms on the right hand side will be 
absorbed by the terms available on the left for sufficiently small A (depending 
only on M). Note the importance of Lemma 12.12 to estimate this error-term, 
as it allows one to turn the non-degenerate into a degerate energy. 

Collecting our estimates for ^"2026' £1 2a 261 Proposition 12.14 is proven. 

□ 

Finally, we note that for the wave equation satisfied by a (and Lie-T- 
derivatives thereof), the terms Gi and G2 can be estimated in exactly the same 
fashion as J-" and J-"^ were for the p-equation. In fact, the estimates are much 
easier, as the "linear" error-terms on the right hand side all cancel. This fact 
is also easily inferred from commuting the scalar equation of Corollary 12.10 
directly with T and using that a decays. We conclude 

Proposition 12.20. The estimate of Proposition 12.14 also holds replacing T 
by Q everywhere. 



12.6 Proof of Proposition 12.1 

Combining Proposition 12.13 with the error-estimate 12.14 yields the estimate 
of Proposition 12.1 for the case n > 1. For the case n = wc run the argument 
of Proposition 12.13 again, since 0o satisfies the same equation by Corollary 
12.8. The only difference is that Lemma 12.12 is now not available. Instead, 
we simply control the zeroth order curvature term from the energies on the 
Ricci-coefficients: 

/ -^dt*drduj <B-n^ [% (n, T2) , (439) 

as a consequence of the null-structure equations. 



13 Proofs of the main theorems 
13.1 Uniform boundedness: Theorem 5.1 

Uniform boundedness for the top-order derivatives can be proven without prov- 
ing a (degenerate) integrated decay estimate at this order. This observation was 
first made in the context of the wave equation on black hole backgrounds in [4] 
(see also [7]) and is adapted here to the spin2-setting. We are going to use the 
slices E (no "ears" ) is this section. 

Proof of Theorem 5.1. To lighten the notation we write E, = E,-. , E^ = E^. n 
{r > ry - y-^^} and E" = E,, n {r < ry - 3^:^^}, = % (ti,T2), 
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= Ai (ti,T2). The identity (60) applied with three T-vectorfields yields, in 
view of Lemma 6.1 



W 



71 



w 



W 



3 = 1 



n 

(Si)+i?.^i?[4iy] (Si) 



W 



i=i 



j^TTT 



K. 



TTT 



Using the elliptic estimate in the interior, Proposition 9.1, we obtain 



[ly] (S2) < B • E" [ly] (Si) + e • E;<,^ [W^] (S2, Si) 



j = l "'^l? i=l,2 



T7-_ 1 

Bsupi; [«H] (S^) . (440) 



Adding the redshift-estimate of Theorem 8.1 and the expression 
Xvifn{r>rv} ^ [^] (^**) t° ^i'l^s yields 

1" ^\ (S2) + 1" [M^] (H?) + / d<*E" [W^] (St* )<B-T^\ (Si) 



^TTT 



K. 



TTT 



+BsupE"' ^[5n](S^) + B 



Mln{r>rY} 



B-Errl,, [5H](Ti,r2) 



Note that the left hand side controls in particular the time-integrated L^-based 
energy. For the errorterms Kf^^ and A' J^-^ we use the following Corollary of 
Proposition 7.9 (for the untilded region Ai {ti,T2)): 

Corollary 13.1. IfX=y = Z = Tin Proposition 7.9 we have, for any 
(small) A > and < n < k — I the estimate 



K 



TTT 



< A sup E"+' [W] (SO + Bx ■ C"+i m (n, T2) 



where the C"" energy was defined in Theorem 5.1. Moreover, in view of 



Proposition 6.2, the same estimate holds for K]^^^ 
the last term). 



( in fact, without 



Applying also the estimate for Err^^^ (ri,r2) given in Proposition 8.1, 
we obtain 

E" [W] (S2) + 1" [W] (Hi) + I dt*t" [W] (St.) < B ■ 1" [W] (Si) 



+B-C"[<n](ri,T2) + B / dt* 



r[TV] (St.) 
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We finally estimate the last term of the above estimate using the estimate 

(440) integrated in time.'^° 

Defining / (t) = E [W] (St), the resulting estimate can be written in the 
following form: For any A > and T2 > ti > tq 

f (t2) + r dtf (t) < Bx ■ (C" m (to) + / (to)) • max [1, (t2 - n)] 

+Amax[l,(T2-Ti)] sup / (t) + B • / (n) . (441) 

■re(Ti,T2) 

In particular, we can achieve that A < j^^- Boundedness of / (t) now follows 
from Lemma 13.2. □ 

Lemma 13.2. Suppose / (to) < oo and that f (t) satisfies the estimate (44V 
with A < j^gj < jq and B\ = B > \ fixed. Then supj-^^^ o^-j f (t) < C holds for 
a constant C = 2 • B • / (to) + 9 • 8 • • B (C" [^H] (to) + / (to)). 

Proof. We will write D ~ C" (to) + / (to). The proof is a bootstrap argu- 
ment. Consider the region A{t) = Ai (tq, t) for all t such that f (t) < C holds 
in A{t). This region is clearly non-empty and open as a subset of the black 
hole exterior J\4 {tq, oo). We now show it is also closed using the estimate (441). 
For this, let L = AB^ and define a sequence of slices t; — tq + i ■ L. Applying 
the estimate in a region between two slices Ti,Ti located in the region A(t) 
yields a "good" slice rf satisfying 

/(Tf)= inf fiT)<B-D + X-C+^-C. (442) 

Having found a good slice in each region M. {ti, Ti+i), we can apply the estimate 

(441) again, this time from the good slice to any T-slice in Ai (t^ ,Ti+2) (and, 
for the last good slice, in Ai {t^,t)). Note that the distance between the two 
slices considered is at most 2L. Hence the estimate yields 



sup f{T)<B-D-2L + X-2L-C + B-[B-D + X-C+^-C 

(^f.^.+2) 



L 



o ~ 1 , ~ 1 1 15 

<8B^-B-D + -- C + B^B-D + —C + -C < —C . 

2 16 4 16 

This uniform bound has been established for all slices in A (t) north of the first 
good slice, Tq . However, for the region up to Tq our estimate tells us that 

sup f {t) < B ■ f {to)+4:- B^ B ■ D + lc <^C . (443) 
/ 3^ 4 4 

It follows that A (t) is closed and that the uniform bound we established is valid 
for the entire black hole exterior. □ 



^•^It is very important that (440) does not have an error- term on the horizon for this argu- 
ment to go through. This is the underlying reason for the assumption that T is null on the 
horizon: It ensures both that the energy flux through the horizon has a sign and that highest 
order error- terms are absent on the horizon in Proposition 7.9 . 



Ill 



13.2 Integrated decay: Theorem 5.2 

Proof of Theorem 5.2. Apply Proposition 11.2 with Wi = C^W for 1 < i < fc 
and add the estimate of Corohary 11.4 at the lowest order. Apply Theorem 12.1 
to estimate the p and cr-term on the right hand side. This yields, for 1 < n < k 



n 



(t"1 ,T2 



(ti,T2) 



{TI.T2) 



+E [W,] {n (ri, T2)) + r'^'-' {m (ti, T2)) } + f [W] {^M (ri, T2)) 

n 



4=1 



W] (M (ti,T2) +supE" [W] (S, 



+ 5 -©5; [5n] (ti,T2) 

+bY, Err Mora im 

i=l 

where we wrote $j to indicate that we are adding both the estimate for (jji and 
ipi of Theorem 12.1. For the non-degenerate case we have 



n— i 

sup [W] (tr) + sup EspM [W.] (tr) + sup E [2?*,_i] (s 



(ti,T2) 



(ti,T2) 



(ti,T2) 



+S [W.] (ri, T2)) + /""""^-^^ [W,] (ri, T2)) } + f [T4^] (>i (n, T2 

n— 1 n 
i=l i=l 

T~^"' [W] (^M (ri, ra)) + supE" [W] (s. 



n-1 

+B J2 Err Mora im . 
1=1 

(444) 



It is important to note that we only need the n-derivative energy of the com- 
ponents p and (7 (and an e of the n-derivative energy of all components) to 
obtain a non-degenerate integrated decay estimate for n — 1-derivatives of all 
components of the Weyl-tensor. As a consequence, this estimate requires only 
Err Mora [VV„_i]. Note also that for n = 1 in (444) there is no ErrMora-ie^im. 

We next turn to Propositions 7.9 and 11.1 to estimate Err Mora [Wi]. We 
observe that for A sufficiently small, the error-term on the horizon (the last term 
in (285),) can be absorbed by the horizon-term available on the left of (444). 
(In particular, the assumption that T is null one the horizon is not necessary 
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for this argument to work.) We conclude, therefore, 



sup E° [W] (tr) + V I sup [m] itr) + SUp E 

- ^ ^ (ti,T2) ^ ^ (ri,r2) ^ 



('''l,'r2) 



+S [W.] (ri, ra)) + 7'^'=^' [W,] (n, T2)) } + f [W^] [m (n, T2)) 

n 



+ (A + e) 



-n,deg . 



Ip '[W]lMiTi,T2)]+SUpE [W]{^r 



+ B-©"[5n] (ti,T2) 



Note in particular that we need and not only D^' [5H]. On the contrary, 

because the non-degenerate estimate needs only Erri\jora [VVn-i], we have 



sup [W] (tr) + y I sup E,p,nl (±r) + SUp S[2?*,_i] (s, 

- ^ ~[ ^ (ruT2) ^ ' (ri,T2) ^ 



(Tl,-r2) 



+i? [W,] (ri, T2)) + /"""'^^s [W,] (ri, T2)) } + f [14^] (>J (n, T2) 

n— 1 n 

" ^ i=l 

[ly] [m (ti,T2)) +supE"[M/] (S, 



—n.deq , 

ip 



i3-D5; [5K](ri,T2) , 



which requires only B" [91] on the right hand side. 

These estimates provide control over the Lie-T derivatives of all components. 
In view of the elliptic estimates of section 9 this is sufficient to control all 
derivatives away from the horizon and from infinity (r < R + AI). In particular, 
defining'^^ 

ry - 2M -1 



ry - ^-^ <r<R + Mj, 



^T,int = Et n |ry 
Mint {ti,T2) = M (ti,T2) H | 

we have 

sup r [W] (tr,^nt) + l"'"'" [W\ (M^nt (ti,T2) 

(ti,T2) ^ ^ ^ 

n 



(445) 
(446) 



I^f ' [W] (ri,T2)) + supE" [W] (e, 



B-©"[5n](Ti,T2) (447) 



^^The reason that S^^int only reaches up to R is that S starts being null for r > i?, and 
standard elliptic estimates are no longer available. 
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and 



—n — l.nondeg . f .~ . , 
I [M^] (Xmt(Ti,r2) 



Tl-l 



j=l 



I^''''" [ly] ( (ri, T2) ) + supE" [VK] ( 



We now couple the degenerate estimate (447) with the redshift estimate of 
Proposition 8.1 and with the estimate near infinity, Theorem 10.1, to obtain 
Theorem 5.2. □ 



13.3 The decay iteration: Theorem 5.5 

Proposition 13.3. Recall the decay matrices Pj defined in (192)-(194)- For 
j = 0, 2, we have for any i = 0, k — 2 the estimates 



Ep+;^ [W] (S.) < I [sup Ep+' [W] + B - 



■'+2 r$H] (r, oo) 



(448) 



Ct^] (■^(r,oo)) 



< 



B 



m 



B 

T 



sup [W] (tr') + B ■ ©^+2 [9^] oo) 



and 



[W^] (r. 



oo < 



B ■ 



m] (r) 



sup E'^^ [W] [Y.r'j + B 



t'>t 



[9^](r,oo) 



(449) 



(450) 



Before we prove the above Proposition let us note how it implies Theorem 

5.5: 

Proof of Theorem 5.5. Theorem 5.5 follows from repeatedly (depending on the 
decay matrix) inserting the estimate (448) into (449) and (450). □ 

Proof of Proposition 13.3. At order i + 1, we apply the estimate (203) with 
decay matrix Pq, cf. (191), which has the associated bulk admissible matrix 

Q-5 Q-5 b-5 3-(5 1-5 
p^^l 9-(5 %- 6 8 - 5 7- 5 5-5 3-5 1-5 



From the resulting estimate we derive 



Ept' [W] (S.„^,) + ^ di*Ep+' [W] (S.) < B ■ [W] (S.„) 



-B 



[5n](T„,T„+i) , (451) 
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for the boundary admissible decay matrix Pi defined in (192), where we only 
used that Pi < Pq < Pq (recall the ordering relation defined in section 4.2). We 
conclude 



Pi 



E 



sup Ep+' [W] (tr) + B 



y+2 



[$n](r,c 



(452) 



which is (448) for j ~ 0. Since the square bracket is bounded (by applying 
(202) from data with Pq), we have established j decay of a weighted energy. In 
particular, since Pp < Pi we have shown this decay for the Pp-wcighted energy 
at order i + Applying (202) with Pi yields, after inserting (452), the estimate 
(449) for j = 0. In the same way, applying (203) produces (450) for j ~ 0. 

It is clear that this procedure can be iterated. We next apply our estimate 
(202) with the boundary admissible decay matrix Pi and its associated bulk 
admissible matrix 



Pi 



/ 5-S 5-(5 3-5 1-S 

8-S 8-5 8-S 7 - 5 5 - 5 3 - 5 1 -5 



Noting that P2 < Pi < Pi for P2 as in (193), we obtain, after adding the 
non-degenerate integrated decay estimate, the analogue of (451) 



[W] (s.„,,) 



dt* Ep+^ [W] ( 



< B • Ep^ [Vt^] + S • [m] (r„, T„+i) , 



(453) 



from which wc conclude 

Ep+' [W] (tr) < ^ [ sup rp+' [W] (tr') + B 



[9^](r,oo) 



(454) 



Another iteration, applying (202) with P2, will provide the estimate 



Ept' [W] {±r) < I [ sup Ept' [W] (tr') + B 



(455) 



and hence an r- weighted energy which decays stronger than -ij . This is enough 
to close the argument. We remark that, of course, there is no reason to stop the 
iteration at this point. To obtain the optimal result in terms of the decay in t, 
one can define additional admissible decay matrices and continue the iteration. 

□ 
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A Useful formulae 

Let (a) denote tetrad (=frame) indices, while Greek letters denote spacetime 
indiees. The Rieci eoefBcients are defined as 

7(c)(a)(b) = e(a)K;i.e^^e^^) (456) 

with the nuU-eomponents 

7(B)(3)(A) = Kab ' 7(s)(4)(A) = Hab , (457) 

7(yi)(3)(3) = 211^ , 7(A)(4)(4) = ZYa, (458) 
7(A)(3)(4) = , 7(^)(4)(3) = 2Za , (459) 

7(4)(3)(3) = 4il , 7(3)(4)(4) = 4fl , (460) 

7(3)(4)(A) - 2Va . (461) 

Due to the antisymmetry in the first two indices we can write e^.^^''.^ = ~7(q')''^- 
For the Weyl-tensor we note the relations 

W(am(c)(d) = e'(af'(b)e1c)eJd)W,,^<rT (462) 

and 

D(a}W(^b){c){d}if) = D{a) (W^(&)(c)(d)(/)) -7(ff)(b)(a)W^^''\c)(d)(/) 

Glossary 

-^3"' (q^) inhomogeneities in the (commuted) Bianchi 

equations, 25 

J {X, W) error-term arising in the context of commuta- 

tion with a vectorfield X, 22 

X^y^ [W], K^^^ [W] spacetime terms in the main-energy identity 

(60), 20 

L 2L = — (2 — p) 63+964 (L ~ dr in Schwarzschild 

in (<*, r)-eoords), 99 

N N = T + Y , everywhere timelike vectorfield, co- 

incides with T away from the horizon, 45 

Nout {S'^ r) outgoing null-hypersurface emanating from 

Sir. 14 

Pj boundary-admissible decay matrix, 37 

R r = R is a fixed timelike hypersurface close to 

null-infinity {R is large) , 14 
T 2T = P63 + ultimately Killing field, 16 
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T"*- 2T^ = —pes + 964, vcctorfield (non-dcgcncrate 

derivative in integrated decay estimate for p and 
cr), 38 

Y redshift vectorfield (and, unfortunately, a null 

component of the Ricci-rotation-coefScients, 
cf. the appendix. Since we work in a gauge 
where K = 0, no confusion can arise.), 45 

Et* slices of constant t*, 14 

ED^ (ti,T2) A-wcighted energy on the Ricci-coefficients, 19 

D" [5H] (ri, T2) norm on the Ricci-coefficients, 17 

SH, 9^™«™ collective notation for Ricci-Rotation coeffi- 

cients, cf. also the appendix, 16 

[W] (^^T^ weighted curvature energy, weight encoded in 

decay matrix P, 36 
Ip [W] iU) weighted spacetime curvature energy, weight en- 

coded in decay matrix P, 36 
Pi bulk-associated decay matrix, 37 

Et* slices of constant t* for r < P with null hyper- 

surfaces Nout {S'^ fl) attached, 15 
jCx modified Lie-derivative to make CxW traceless, 

21 

k^,k^, k~ functions of r, 13 

p 3-component of the ultimately Killing field 

T, vanishing on the horizon, equal to in 
Schwarzschild, 16 

q 4-component of the ultimately Killing field T, 

k+ 

equal to 1 on the horizon, equal to -7^ in 
Schwarzschild, 16 

Ty r = ry is a fixed timelike hypersurface close to 

the horizon (see Lemma 6.3) , 15 
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